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1.  INTRODUCTION 


Recent  years  have  seen  a  rapid  increase  in  research  on  quantum-well  physics, 
materials,  devices  and  theories.  The  emergence  of  quantum-well  electronics  has 
ushered  in  a  new  era  in  tailoring  materials  and  structures  for  specific  purposes,  and 
now  offers  the  possibility  to  advance  the  state  of  the  art  of  submillimeter  technology 
and  a  means  to  study  high  speed  processes  in  semiconductors.  This  activity  is  indica¬ 
tive  of  the  great  need  for  ultra-fast  active  components.  The  program  summarized  here 
has  developed  some  new  approaches  to  fill  this  need. 

At  Lincoln  Laboratory  we  have  placed  our  primary  effort  in  the  recent  ARO  pro¬ 
gram  on  the  area  of  resonant-tunneling  devices.  We  have  studied  the  physics  of 
resonant  tunneling  and  developed  oscillators,  multipliers,  and  mixers  into  the  submil¬ 
limeter  spectrum.  Preliminary  work  into  the  feasibility  of  resonant-tunneling  transistors 
has  also  been  undertaken.  Planar  antenna  arrays  have  been  investigated  initially  for 
application  to  detector  arrays,  but  they  may  also  have  application  for  oscillator  arrays. 
The  results  from  all  these  studies  have  been  so  encouraging  that  we  are  planning  to 
continue  development  of  these  devices  and  to  investigate  new  phenomena  involving 
coupled  quantum  wells. 
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2.  RESULTS 

Progress  during  the  contract  period  has  been  reported  in  the  semiannual  progress 
reports  to  ARO.  Below  are  listed  significant  highlights,  the  details  of  which  appear  in 
the  reprints  and  preprints  that  are  appended  to  this  report.  In  most  cases  the  work 
received  additional  support  from  other  sponsors.  An  overview  of  some  of  the  devices 
developed  under  this  program  can  also  be  found  in  Appendix  A. 

2.1.  Large  Negative  Differential  Conductance  at  Room  Temperature 

Early  resonant-tunneling  structures  required  cooling  to  achieve  negative 
differential  conductance.  However,  by  using  thin  AlAs  barriers  instead  of  the  AlGaAs 
alloy  and  by  reducing  the  impurity  concentration  outside  the  barriers,  we  have  obtained 
peak- to- valley  ratios  of  3.5:1  and  current  densities  of  4x10^  AJcvc?  at  room  tempera¬ 
ture.  This  development  has  greatly  simplified  the  fabrication  and  testing  of  a  variety 
of  devices,  in  particular  oscillators,  and  has  expedited  other  investigations  of  the  phys¬ 
ics  of  resonant  tunneling.  (See  Appendix  B  for  details.) 

2.2.  Room-Temperature  Oscillators 

We  have  achieved  fundamental  oscillation  frequencies  up  to  420  GHz  using  the 
negative  differential  conductance  region  of  the  resonant-tunneling  current-voltage  (I-V) 
curve.  (See  Appendix  C  for  results  up  to  200  GHz.  The  420  GHz  results  are  in 
preparation.)  This  oscillation  frequency  confirms  our  predictions  of  the  extrinsic  max¬ 
imum  oscillation  frequency  for  the  double-barrier  structure  used  in  these  experiments 
and  surpasses  any  tunnel  diode  result  reported  in  the  literature.  However,  the  power 
measured  at  200  GHz  is  only  0.2  |iW  for  this  material,  largely  because  this  frequency 
is  very  close  to  the  maximum  oscillation  frequency  of  270  GHz.  At  lower  frequencies 
(<  60  GHz)  we  have  measured  powers  between  50  and  200  pW.  Although  these 
powers  are  insufficient  for  use  as  a  local  oscillator  with  a  room-temperature  Schottky 
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diode  mixer,  they  are  more  than  enough  for  a  local  oscillator  with  superconductor- 
insulator- superconductor  (SIS)  mixers  operating  at  4  K.  We  have  also  measured  other 
properties  of  the  resonant-tunneling  oscillator  specific  to  its  application  as  a  local  oscil¬ 
lator.  Frequency  tunability  of  over  10%  has  been  demonstrated  near  100  GHz,  and 
injection  locking  is  possible  with  injected  power  20  dB  below  the  peak  oscillator 
power.  This  latter  fact  and  specific  noise  measurements  indicate  that  resonant¬ 
tunneling  oscillators  should  be  quite  low  in  noise  content. 

2.3.  Self-Oscillating  Harmonic  Mixer 

The  resonant-tunneling  diode  has  the  capability  to  act  as  an  efficient  mixer 
because  of  the  rapid  variation  of  the  dynamic  conductance  with  voltage  near  the 
negative-differential-resistance  region  of  the  I-V  curve.  For  a  diode  oscillating  with  a 
DC  bias  at  the  center  of  this  negative-differential-resistance  range,  the  time-varying 
dynamic  conductance  actually  has  a  predominant  Fourier  component  at  the  second  har¬ 
monic  of  the  oscillation  frequency.  This  implies  that  the  conversion  efficiency  should 
be  high  for  signals  at  twice  the  oscillation  frequency.  We  have  experimentally  studied 
the  properties  of  this  harmonic  conversion  with  a  resonant-tunneling  structure  that 
oscillated  at  50  GHz.  The  conversion  loss  of  a  100-GHz  signal  was  found  to  be  about 
12  dB,  which  is  near  the  state  of  the  art  for  harmonic  mixers.  This  suggests  that  self- 
oscillating  resonant-tunneling  mixers  could  be  very  useful  in  the  millimeter-wave  band, 
perhaps  replacing  conventional  Schottky  diode  mixers  in  applications  where  the  use  of 
a  separate  local  oscillator  is  cither  difficult  or  imdesirable.  (See  Appendix  A  for 
details.) 

2.4.  Multipliers 

We  have  used  double-barrier  resonant-tunneling  structures  as  microwave  and 
millimeter-wave  multipliers.  The  antisymmetry  of  the  I-V  curve  about  the  origin,  i.e.. 
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I(V)=-I(-V),  combined  with  the  presence  of  a  negative  differential  resistance  region, 
suggests  that  fifth-harmonic  generation  should  be  especially  efficient  with  a  DC  bias 
voltage  near  zero.  We  have  observed  fifth-harmonic  efficiencies  of  -23  dB  multiplying 
from  4  to  20  GHz  without  optimized  device  or  load  impedances.  Theoretical  predic¬ 
tions  suggest  that  efficiencies  slightly  greater  than  -10  dB  should  be  obtainable.  The 
maximum  output  frequency  should  be  on  the  order  of  the  maximum  oscillation  fre¬ 
quency,  as  discussed  in  Appendix  D. 

2.5.  Ultimate  Frequency  Limits  for  Resonant-Tunneling  Diodes 

We  have  calculated  the  times  required  for  charge  transport  through  various  parts 
of  the  double-barrier  structure,  including  the  depletion  region.  Using  the  experimen¬ 
tally  observed  oscillation  frequencies,  we  have  shown  that  the  WKB  approximation 
greatly  overestimates  the  storage  time  of  charge  in  the  well,  and  that  a  calculation 
method  based  on  the  energy  width  of  the  resonant  state  gives  results  consistent  with 
the  measurements.  Also,  the  transit  time  across  the  depletion  region  can  be  important 
at  higher  frequencies,  as  can  the  RC  time  determined  by  parasitic  series  resistance  and 
device  capacitance.  This  theoretical  foundation  provides  a  guide  to  design  of  higher 
frequency  structures  and  suggests  that  there  are  no  fundamental  or  practical  impedi¬ 
ments  to  development  of  resonant-tunneling  oscillators  up  to  several  hundred  gigahertz. 
(See  Appendix  E  for  details.) 

2.6.  Persistent  Photoconductivity 

We  have  made  the  first  observation  of  persistent  photoconductivity  in  resonant¬ 
tunneling  structures.  Illumination  at  low  temperatures  of  structures  with  Alo.3Gao7As 
barriers  results  in  a  persistent  shift  of  the  current-voltage  (I-V)  curve.  This  is 
explained  by  modifications  to  the  conduction  band  from  fixed  positive  charges  in  the 
barriers  and  equal  negative  charges  that  accumulate  outside  the  barriers.  The  fixed 
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positive  charge  arises  from  ionization  of  DX  centers  located  in  the  barriers.  Calculated 
I-V  curves  agree  with  our  measurements,  and  the  derived  DX  center  density  agrees 
with  estimates  from  material-growth  constants.  This  effect  has  been  useful  in  deter¬ 
mining  the  heterojunction  parameters  (as  described  in  Appendix  F)  and  in  estimating 
the  performance  of  the  structure  in  which  the  quantum  well  is  contacted  to  form  a 
resonant-tunneling  transistor.  In  the  latter  case,  the  change  of  the  well  potential  with 
respect  to  the  contacts  (caused  by  the  DX  center  ionization)  permitted  an  estimate  of 
transistor  transconductance  and  cutoff  frequency.  (Appendix  H  has  details  of  the 
transistor.) 

2.7.  Resonant-Tunneling  Transistor 

A  resonant-tunneling  transistor  has  been  proposed  and  fabrication  feasibility  stu¬ 
dies  have  begun.  The  persistent  photoconductive  effects  discussed  above  have  made  it 
possible  to  estimate  the  transconductance  and  unity-current-gain  frequency  fp  for  a 
transistor  based  on  the  double-barrier  resonant-tunneling  structure,  in  which  the  well 
serves  as  the  control  electrode.  With  conservative  assumptions  based  on  experimen¬ 
tally  obtained  current  densities,  the  projected  transconductance  gn,  for  a  l-p,m  emitter 
stripe  is  2000  mS/mm.  This  corresponds  to  a  fp  of  at  least  50-100  GHz  and  a  max¬ 
imum  frequency  of  oscillation  of  150-2(X)  GHz.  Moreover,  an  increase  in  current 
density  beyond  that  presently  attained  by  over  a  factor  of  10  is  predicted  theoretically, 
which  should  increase  the  fp  to  well  over  1(X)  GHz.  (See  Appendix  F  for  details.) 
We  have  developed  a  technique  for  growing  thicker  epitaxial  layers  that  are  attached  to 
the  thin  quantum  well  by  growing  on  the  sidewall  and  horizontal  surfaces  of  an 
etched-groove  structure.  Horizontal  barriers  and  wells  nearly  three  times  as  thick  as 
on  the  sidewalls  have  been  achieved.  This  should  simplify  the  fabrication  of  contacts 
to  the  thin  quanmm  wells. 
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2.8.  Picosecond  Optical  Measurements 

The  switching  time  of  one  of  our  resonant-tunneling  diodes  (with  oscillation 
cutoff  frequency  of  270  GHz)  has  been  measured,  using  optically  generated  picosecond 
electrical  pulses.  The  measurement  was  made  in  collaboration  with  Gerard  Mourou 
and  his  students  at  the  University  of  Rochester.  (See  Appendix  H  for  details.)  The 
observed  switching  time  of  2  ps  is  in  agreement  with  an  approximate  theory  of  switch¬ 
ing  time:  x  =  C(AV/AI),  where  AV(AI)  is  the  voltage  (current)  excursion  of  the  nega¬ 
tive  differential  resistance  region.  These  experiments  confirm  our  earlier  measure¬ 
ments  in  the  frequency  domain,  and  suggest  that  the  device  could  act  as  a  high-speed 
switch.  Promising  digital  applications  are  being  explored. 

2.9.  Planar  Dipole-Fed  Mixer  Arrays 

A  planar  double-dipole  antenna  containing  a  beam-lead  Schottky  diode  has  been 
developed  to  produce  the  best  sensitivity  of  any  planar  mixer  above  100  GHz 
[Tjec  =  49(X)  K  (DSB)]  at  140  GHz.  This  element  has  also  been  incorporated  into  a 
planar  two-dimensional  array  that  performed  well  for  element  spacing  of  about  two 
wavelengths.  (See  Appendix  I  for  details.)  With  the  Schottky  diode  replaced  by  a 
resonant-tunneling  diode,  the  double-dipole  could  act  as  a  resonator  for  oscillations  in 
the  submillimeter  region.  If  successful,  the  oscillating  elements  could  then  be  com¬ 
bined  in  an  array  to  greatly  increase  the  output  power. 

2.10.  Rocket-Plume  Spectroscopy  at  557  GHz  with  a  Compact  Solid-State 

Radiometer 

We  have  developed  a  compact,  solid-state  submillimeter-wave  heterodyne  radiom¬ 
eter  and  have  used  it  to  measure  spectral  characteristics  of  a  simulated  rocket  exhaust 
plume.  Features  of  the  557-GHz  water-vapor  line  profile  were  observed  in 
significandy  greater  detail  than  in  previous  experiments  through  increased  sensitivity 
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and  improved  frequency  resolution  (600  kHz).  A  prominent  self-reversal  feature  from 
a  warmer  central  core  of  the  plume  was  observed.  The  local  oscillator  of  the  radiome¬ 
ter  utilized  a  frequency  multiplication  chain  consisting  of  an  InP  Gunn  oscillator  at 
92.6  GHz  driving  a  frequency  tripler  and  second-harmonic  mixer  in  cascade.  The 
receiver  had  a  noise  temperature  of  4500  K  (DSB)  at  555  GHz,  consumed  3  W  and 
weighed  3  kg.  (This  project  was  carried  out  in  collaboration  with  N.R.  Erickson, 
University  of  Massachusetts,  Amherst.  See  Appendix  J  for  details.) 


APPENDIX  A 


T.C.L.G.  SOLLNER,  E.R.  BROWN,  and  H.Q.  LE 


Microwave  and  Millimeter-Wave 
Resonant-Tunneling  Devices 


Resonant-tunneling  devices,  which  may  be  capable  of  operation  at  terahertz  frequen¬ 
cies,  have  been  developed  and  tested.  Included  in  these  solid  state  microelectronic 
components  are  oscillators,  self-oscillating  mixers,  and  harmonic  multipliers.  A  char¬ 
acteristic  of  these  devices,  negative  differential  resistance  (NDR),  has  been  observed  at 
room  temperature.  Resonant-tunneling  transistors,  which  promise  operation  in  the 
terahertz  frequency  range,  are  also  proposed. 


PHYSICS  OF  RESONANT  TUNNELING 

A  resonant-tunneling  diode  can  respond  to 
electrical  impulses  in  picosecond  or  subpico¬ 
second  times.  Therefore,  these  devices  may 
provide  a  basis  for  developing  electronic  devi¬ 
ces  that  operate  at  terahertz  frequencies. 

The  essential  features  of  resonant  tunneling 
are  shown  in  Fig  1.  A  thin  layer  of  GaAs  (2  to  10 
nm)  is  sandwiched  between  two  thin  layers  of 
Al^Ga^.^As.  The  addition  of  aluminum  to  GaAs 
raises  the  band-gap  above  that  of  GctAs  so  the 
Al^Gaj.jfAs  regions  act  as  partially  transparent 
mirrors  to  electrons;  the  higher  energy  level  of 
these  barriers  reflects  the  electrons  back  to  the 
region  of  the  structure  from  which  they  came. 
The  charge  transport  across  the  structure  takes 
place  by  tunneling  through  the  thin  (1  to  5  nm) 
Al^Ga^.^As  barriers. 

This  structure  is  the  electron  analog  of  a 
Fabry-Perot  resonator.  As  shown  schematic¬ 
ally  at  the  bottom  of  Fig.  1,  the  resonator 
exhibits  peaks  in  the  electron  transmission 
(current)  as  the  incident  electron  energy  (vol¬ 
tage)  changes. 

The  physical  implementation  of  these  prin¬ 
ciples  is  summarized  in  Fig.  2.  The  la3rered 
material  is  grown  in  wafer  form  by  molecular 
beam  epitaxy.  Then  the  active  regions  are  de¬ 
fined  with  ohmic  contacts.  These  contacts  are 
used  as  a  mask  to  isolate  the  region  under  the 
contact,  either  by  etching  mesas  (as  shown  in 
Fig.  2),  or  by  proton  implantation,  which  makes 
the  surrounding  material  nonconductive.  Be¬ 
cause  the  contact  is  only  a  few  microns  in 
diameter,  electrical  connection  to  the  ohmic 
contact  is  made  with  a  pointed  wire  (whisker). 


In  some  cases,  the  whisker  acts  as  an  antenna, 
coupling  high-frequency  ac  fields  to  the  double- 
barrier  diode. 

Using  the  Fabry-Perot  analogy,  coherence  of 
the  electron-wave  function  is  required  across 
the  entire  double-barrier  region  to  maintain 
resonant  tunneling.  Any  scattering  that  occurs 
in  either  the  well  or  the  barriers  will  alter  the 
wave  function  phase  randomly,  destroying  its 
coherence  and  therefore,  the  conditions  re¬ 
quired  for  resonant  tunneling.  But  a  different 
picture  (see  Appendix.  “Resonant  Tunneling 
Theoiy"),  one  that  does  not  require  coherence 
between  the  parts  of  the  wave  function  outside 
and  inside  the  well,  can  produce  negative  resis¬ 
tance,  the  essential  characteristic  that  is  ex¬ 
ploited  in  the  devices  that  are  described  in  this 
article. 

In  1973,  Tsu  and  Esaki  (1)  derived  the  two- 
terminal  current-voltage  (I-V)  curves  for  flnite 
multiple-barrier  structures.  This  matching  tech¬ 
nique  has  been  remaikably  successful  in  explain¬ 
ing  experimental  results.  In  1974,  Chang  et  al 
(2)  were  the  first  to  observe  resonant  tunneling 
in  a  monocrystalline  semiconductor.  They  used 
a  two-barrier  structure  and  observed  the  reso¬ 
nances  in  the  current  by  measuring  the  I-V  ciuve 
of  the  structure.  The  voltages  at  the  current 
peaks  agreed  well  with  Chang's  calculations. 

A  decade  later,  interest  in  the  field  was 
renewed  when  Sollner  et  al  [3]  showed  that  the 
intrinsic  charge  transport  mechanism  of  a  two- 
barrier  diode  could  respond  to  voltage  changes 
in  less  than  0.1  ps  (>  1  THz).  More  recently, 
negative  differential  resistance  (NDR),  a  char¬ 
acteristic  of  resonant  tunneling,  has  been  mea¬ 
sured  at  room  temperature  [4]. 
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Fig.  11a)  —  The  energy  profile  of  an  electron  in  a  double¬ 
barrier  resonant-tunneling  structure  is  shown  here.  The 
top  illustration,  which  includes  the  energy  well  pro¬ 
duced  by  the  twoAlxGai.xAs  layers,  shows  the  unbiased 
energy  profile.  The  Alfia^.^As  layers  act  as  partially 
transparent  mirrors  to  the  electrons,  similar  to  a  Fabry- 
Perot  resonator,  (b)  In  this  plot  of  diode  current  as  a 
function  of  incident  electron  energy,  the  large  resonant 
current  at  point  A  corresponds  to  the  energy  profile  of  A; 
the  valley  in  the  current-profile  at  point  B  corresponds  to 
energy  profile  of  B.  Although  the  energy  of  incident 
electrons  is  higher  at  B,  the  absence  of  resonance  low¬ 
ers  the  current. 


SPEED  OF  RESPONSE 

Two  factors  indicate  that  the  response  time 
of  a  double-barrier  resonant-tunneling  struc¬ 
ture  should  be  as  little  as  0.1  ps:  the  resonant- 
state  lifetime  of  an  electron  and  the  time  re¬ 
quired  for  a  double-barrier  structiu^  to  reach 
equilibrium  after  an  electrical  impulse.  This 
predicted  response  time  is  the  time  required 
for  the  device’s  current  to  respond  to  a  sudden 
change  in  voltage.  The  comparison  of  the  mea¬ 


sured  and  predicted  times  not  only  gives  an 
indication  of  the  accuracy  of  the  models,  it  also 
indicates  the  fundamental  limits  of  operation 
imposed  by  the  structure  of  double-barrier  de¬ 
vices.  One  of  the  most  significant  limits  of  oper¬ 
ation  is  fn,ay.  the  maximum  useful  frequency  at 
which  the  device  exhibits  NDR.  This  frequency 
is  the  point  above  which  NDR  is  no  longer 
observable  at  the  terminals  of  the  device. 

If  the  voltage  across  a  double-barrier  struc¬ 
ture  is  instantaneously  changed,  the  current 
through  the  device  changes  to  a  different  steady- 
state  value.  As  in  a  Fabry-Perot  resonator,  the 
steady-state  charge  in  the  well  must  decrease 
when  the  current  decreases,  and  increase  when 
the  current  increases.  The  lifetime  of  any  reso¬ 
nant  state,  including  the  one  represented  by  an 
electron  initially  placed  in  the  well  between  the 
two  barriers  of  a  resonant-tunneling  structure, 
is  given  by  t  =  X/AE,  where  AE  is  the  enei^ 
half-width  of  the  transmission  probability  func¬ 
tion  throtigh  the  resonant  state.  It  takes  approx¬ 
imately  one  lifetime  to  fill  or  empty  the  well  to  a 
new  steady-state  value.  Since  the  carrier  trans¬ 
mission  probability  is  determined  by  the  ampli¬ 
tude  of  the  wave  function  inside  the  well,  the 
current  will  reach  its  new  steady-state  value  in 
approximately  the  lifetime  of  a  resonant  state. 

The  lifetime  of  an  electron  in  between  the 
barriers  has  been  calculated  for  three  repre¬ 
sentative  structures  by  SoUner  et  al  (5]  and 
ranges  from  4  ps  for  a  2.5-nm  AlAs  barrier  to 
0.16  ps  for  a  3.0-nm  Ga^  yAlo^jAs  barrier.  In 
Table  1,  this  time  is  found  from  r  =  l/(2n-f^E). 
The  transit  of  electrons  across  the  structure’s 
depletion  region  produces  an  additional  delay 
that  ranges  from  0.16  to  0.69  ps.  Total  delays, 
the  intrinsic  response  times,  therefore  range 
from  0.4  to  5.0  ps. 

For  signals  with  periods  much  shorter  than 
the  intrinsic  response  time,  the  current  lags 
the  applied  voltage  and  the  I-V  curve  for  the 
high-speed  signals  should  deviate  markedly 
from  the  dc  1-V  curve.  The  intrinsic  response 
time  was  first  measured  experimentally  by 
examining  the  difference  between  the  dc  I-V 
curve  measured  and  the  I-V  curve  inferred  from 
high-frequency  measiirements.  Because  the  re¬ 
sponse  times  were  in  the  picosecond  range. 
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Fig.  2  —  This  double-barrier  diode  has  1-  to  5-nm  barriers  and  a  2-  to  10-nm  well. 


operation  at  terahertz  frequencies  was  nec¬ 
essary. 

The  straightforward  approach  to  determin¬ 
ing  the  frequency  response  of  a  device  is  to 
sweep  through  the  I-V  curve  at  increasing  speed 
until  the  measured  I-V  curve  differs  from  the  dc 
curve.  Unfortunately,  it  is  difficult  to  measure 
current  at  frequencies  above  those  accessible 
to  stunpling  oscilloscopes  (100  GHz).  There¬ 
fore,  we  used  a  differential  measiirement  tech¬ 
nique.  We  swept  through  a  large,  low-frequency 
voltage-range  and  measured  the  diode  current 
with  relatively  simple  conventional  techniques. 
At  the  same  time  we  superposed  a  small  high- 
frequency  signal  across  the  diode.  The  low- 
frequency  current  revealed  changes  in  curva¬ 
ture  of  the  I-V  curve  that  were  proportional  to 
the  high-frequency  signal.  Besides  tractabil- 
ity,  this  method  of  measurement  closely  re¬ 
sembles  the  actual  operating  conditions  of  de¬ 
vices  which  utilize  NDR.  In  practical  operation, 
a  device  is  biased  into  the  negative  differential 
region  by  a  large  dc  bias  and  the  superposed 
voltage  doesn’t  extend  far  outside  the  NDR 
region. 

The  circuit  used  for  applying  ac  and  dc  ffelds 
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to  the  double-barrier  diode  is  shown  in  Fig.  3. 
The  comer-reflector  mount  was  originally 
developed  for  Schottky-diode  mixers  in  the  far 
infrared  [6].  Frequencies  between  dc  and  about 
20  GHz  can  be  applied  through  the  coaxial  con¬ 
nector  via  the  whisker  that  contacts  the  active 
area.  The  GaAs  substrate  is  soldered  to  the 
chip  stud,  which  is  at  the  same  ground  poten¬ 
tial  as  the  comer  reflector.  For  frequencies  of 
about  100  GHz  to  a  few  terahertz,  the  long 
whisker  acts  as  an  antenna,  and  the  conduct¬ 
ing  surfaces  cause  images  of  the  whisker 
to  produce  an  antenna  array.  The  array  will  re¬ 
ceive  energy  in  a  beam  that  has  a  frequency- 
dependent  direction,  but  is  approximately  45° 
from  all  siufaces  of  the  cube.  The  full  cone- 
angle  of  the  beam  is  about  20°. 

For  the  2.5-THz  measurements,  for  example, 
the  power  produced  by  an  optically  pumped 
methanol  laser  (about  100  mW)  was  matched  to 
the  antenna  pattern  with  a  lens,  thus  coupling 
about  50%  of  the  incident  power  onto  the 
antenna  and  producing  an  ac  signal  across  the 
double-barrier  structiue.  The  characterization 
of  the  comer-reflector  moimt  at  lower  frequen¬ 
cies  (using  measurements  performed  by  Fet- 
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terman  et  al)  led  to  an  estimate  of  50  0  for  the 
antenna  impedance  Z^.  This  parameter  was 
used  in  Eq.  ( 1 )  to  calculate  the  curves  in  Fig.  4. 

At  high  frequencies  the  current  responsivity 
R|,  which  is  the  change  in  dc  current  ( Al),  divided 
by  the  change  in  ac  power  applied,  (AFg^.),  is 
given  by 


Ri(V)  = 


.M(V) 

•^^ac 


21"(V)Za 

(1+Za/Rs)2 


(1) 


In  Eq.  1, 1"  is  the  second  derivative  of  the  I-V 
curve  at  the  ac  frequency  of  interest.  This 
expression  gives  the  I-V  ciuve  for  current  re¬ 
sponsivity  measured  over  the  voltage  range  of 
interest.  Measurements  were  made  at  1,  138, 
761,  and  2,500  GHz  on  a  double-barrier  diode 
with  3.0-nm  well  and  barriers.  The  results  at  dc. 


1  GHz  and  2.5  THz  are  shown  in  Fig.  4.  The 
measurements  for  138  and  761  GHz  were  essen¬ 
tially  identical  to  the  1-GHz  curve.  Since  the 
1-GHz  signal’s  period  is  long  compared  wiUi 
the  expected  intrinsic  response  time,  deviation 
in  the  I-V  curve  shouldn’t  occur  until  frequen¬ 
cies  that  are  orders  of  magnitude  higher  are 
reached.  The  reason  for  the  change  between  dc 
and  1  GHz  is  not  known,  but  it  is  probably  due 
to  slow  traps  in  the  material. 

By  2.5  THz,  the  diode’s  I-V  curve  looks  quite 
different  from  the  dc  curve.  In  one  direction, 
NDR  has  vanished,  but  it  remains  in  the  other 
direction.  Sollner  et  al  [3]  take  this  as  evidence 
that  the  intrinsic  response  time  for  this  device 
is  of  the  order  of  r  =  (2n-f)'^  =  6  X  10'^'*  s.  This 
result  agrees  approximately  with  theoretical 


Table  1  —  Measured  and  Calculated  Parameters  for 
Three  Different  Wafers  of  Double- Barrier  Diodes 


•Wafer 


1 

2 

3 

Material  Parameters 

Barrier  material 

AlAs 

Gao  7AI0  3AS 

AlAs 

Barrier  thickness  (nm) 

2.5 

3.0 

1.5 

Well  thickness  (nm) 

4.5 

4.5 

4.5 

Doping  outside  barriers  (cm'^) 

1  X10’® 

2X10^^ 

2X10^7 

Electrical  Parameters 

Peak-to-valley  ratio,  300  K 

1.7/1 

1.3/1 

3.5/1 

Peak  current  density  (XI 0^  A 

cm'^)  0.8 

1.2 

4.0 

Depletion  layer  at  bias  (nm) 

15 

30 

70 

Capacitance  (fF)® 

100 

50 

20 

Maximum  negative  conductance  (mS)°  5.0 

8.0 

13.0 

Series  resistance  (D)* 

10 

15 

15 

Oscillation  Characteristics 

DC  bias  Ig,  (mA,  V) 

0.7,  0.40 

2.7.  0.32 

3.0, 0.95 

fosc  (GHz)b 

20.7 

Theoretical 

43.7 

201 

Maximum  Oscillation  Frequency 

Ux  <GHz)= 

35 

70 

270 

Wr<GH2)‘* 

1,000 

500 

230 

f^E  (GHz)« 

40 

1,000 

400 

"Typical  values  for  a  circular  mesa  of  4-/im  diameter. 

^Maximum  observed  fundamental  oscillation  frequency. 

(27rC)-l(-G„„/R,-G2,„)’/2 

^From  depletion  layer  drift  time  assuming  a  drift  velocity  of  10^  cm/s. 

"From  calculation  of  energy  width  of  transmission  through  double-barrier  structure. 
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Intermediate 


Fig.  3  —  This  corner-reflector  mount  is  used  to  apply  signals  to  double-barrier  diodes.  For  frequencies  between  dc  and 
20  GHz.  the  signal  is  coupled  through  the  OSM  connector.  Above  100  GHz,  the  long  whisker  acts  as  an  antenna, 
coupling  the  signals  into  the  chip  from  a  full  cone-angle  of  about  20°. 


expectations  from  Table  1  for  3-nm  barriers  of 
AlojsGao  75As(f_;^E  =  1  X  10^®).  The  asymmetry 
of  both  the  dc  I-V  curve  and  the  response  may 
indicate  that  the  two  barriers  or  the  two  deple¬ 
tion  regions  are  not  identical. 

The  double-barrier  diode’s  equivalent  circuit, 
shown  in  Fig.  5,  provides  a  model  for  further 
understanding  of  this  device.  Included  in  the 
equivalent  circuit  are  the  voltage-dependent 
dynamic  conductance  G(V),  the  series  resis¬ 
tance  R^.  and  the  parallel  capacitor  C  that  is 
inherent  in  the  device  structure.  To  a  good 
approximation,  the  capacitance  is  formed 
across  a  combination  of  the  two  barriers  plus 
the  depletion  region  on  the  anode  side  of  the 
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biased  device.  The  slope  of  the  low-frequency 
1-V  curve  gives  the  approximate  conductance, 
ie,  (dI/dV)-i  =  Rg  +  1/G  «  1/G. 

RESONANT-TUNNELING 

OSCILLATORS 

The  NDR  displayed  by  the  double-barrier 
diode  is  the  basis  of  a  fast  and  simple  two- 
terminal  oscillator.  Interest  in  the  oscillator 
application  stems  from  the  need  for  a  solid 
state  oscillation  source  at  frequencies  above 
300  GHz.  In  this  region,  few  fundamental-mode 
solid  state  sources  are  available,  so  the  double¬ 
barrier  diode  would  be  useful  for  applications 
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that  require  only  modest  amounts  of  power. 
Furthermore,  these  oscillations  give  a  direct 
and  unmistakable  proof  of  the  speed  of  these 
devices.  Of  course,  inherent  circuit  elements, 
particvilarly  the  series  resistance  and  the  device 
capacitance  shown  in  Fig.  5,  must  be  consid¬ 
ered  for  circuit  applications  of  these  devices. 

In  the  fvuther  analysis  required  for  circuit 
applications,  the  conductance,  G(Vq),  is  inde¬ 
pendent  of  frequency  but  strongly  dependent 
on  the  voltage  cunplitude,  Vq,  across  the  device. 
The  real  part  of  the  impedance,  Z^,  measured 
across  the  equivalent  circuit  of  Fig.  5  is  nega¬ 
tive  up  to  a  frequency  given  by  the  expression 


'  2,rC  \  Rs 


(2) 


Where  Gmav  is  the  maximum  negative  value  of 
dynamic  conductance  in  the  NDR  region  of  the 
I-V  curve.  For  all  frequencies  above  fmax* 
real  part  of  the  terminal  impedance  will  be  posi¬ 
tive,  making  it  impossible  for  oscillations  to 
occur. 

Brown  et  al  [7]  found  double-barrier  diode 
oscillators  that  cover  the  frequency  range  of  20 
to  200  GHz  (frequencies  below  f^av).  Figure  6 
shows  the  experimental  results  obtained  in 
this  range  with  diodes  from  the  wafers  listed  in 
Table  1 .  The  initial  experiments  were  performed 
with  a  device  from  Wafer  1  in  a  coaxial  resona¬ 
tor,  giving  an  oscillation  frequency  of  20.7  GHz 
and  an  output  power  well  below  1 W.  Attempts 
to  achieve  oscillation  with  this  device  at  fre¬ 
quencies  near  40  GHz  in  a  WR-22  waveguide 
resonator  were  unsuccessful,  consistent  with 
the  theoretical  f^j^  of  35  GHz.  The  first 
millimeter-band  results  were  obtained  in  the 
vicinity  of  30  GHz  and  40  GHz  in  WR-22  and 
WR-15  waveguide  resonators,  respectively, 
using  a  device  from  Wafer  2,  This  device  could 
have  achieved  higher  oscillation  frequencies, 
but  because  of  the  relatively  low  peak-to-valley 
ratios  of  devices  from  Wafer  2,  no  attempt  at 
higher  oscillation  frequencies  was  made.  The 
wafer’s  low  peak-to-valley  ratio  indicates  that 
devices  from  this  wafer  will  provide  limited 
power  output  and  logic-swing  capability. 

The  voltage  range  over  which  NDR  exists  for 
these  devices  limits  the  output  power  of  oscil¬ 
lators  that  use  them.  Kim  and  Brandli  [8]  and 


Fig.  4  —  The  current-voltage  (I-V)  curves  of  a  double¬ 
barrier  diode  show  the  variation  with  frequency  of  the 
device's  negative  differential  resistance  —  it  essentially 
disappears  at  2.5  THz. 


1  I  ®max  2 

W=  ^|_  Rg  ^max 

D  -  ^  D 

"S  “  Pepi  ^  “spread  “c 

A 

C  =  e  - - 

*“depl 


Fig.  5  —  Using  this  equivalent  circuit,  the  double-barrier 
diode's  maximum  frequency  of  oscillation  can  be  easily 
calculated. 
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Fig.  6  —  The  output  power  versus  frequency  of  reso¬ 
nant-tunneling  diodes  for  three  different  wafers  was 
measured  in  five  different  resonators. 

Trambarulo  [9]  have  shown  that  the  maximum 
output  power,  assuming  a  sinusoidal  voltage 
waveform,  is  Pmax  ~  (3/16)AI/AV,  where  A1  and 
AV  are  the  current  and  voltage  ranges  of  the 
negative  resistance  region.  For  the  4-pm- 
diameter  diodes  of  Wafer  3  in  Table  1,  this  rela¬ 
tion  gives  Pmax  “  ^25  /uW.  Diodes  of  larger  area 
can  produce  more  output  power,  but  Gn,av 
scales  linearly  with  diode  area,  making  it  very 
difficult  to  obtain  dc  stability  in  larger  devices. 
All  of  the  4-/im-diameter  devices  tested  to  date 
have  been  stabilized  with  standard  50  fl  coaxi¬ 
al  loads. 

Diodes  from  Wafer  3  have  produced  the  most 
powerful  oscillations  and  the  highest  oscilla¬ 
tion  frequency  to  date.  This  wafer  extended  the 
maximum  observed  oscillation  frequency  of 
resonant-tunneling  diodes  from  56  to  201  GHz. 
The  output  power  produced  at  201  GHz  was 
about  0.2  /uW,  although  it  reached  60  /uW  at 
lower  frequencies.  Several  of  the  diodes  tested 
did  not  oscillate  at  all.  The  difficulty  in  initiat¬ 
ing  oscillations  near  200  GHz  is  consistent 
with  EJq.  2,  which  predicts  that  fmav  270  GHz 
for  this  diode. 

The  Fig.  5  model  is  identical  to  that  used  for 
p-n  junction  timnel  diodes.  Indeed,  the  double- 
barrier  resonant-tuimeling  diode  and  the  p-n 
jimction  tunnel  diode  are  similar,  displaying 
similarities  in  I-V  characteristics  and  in  circuit 
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behavior.  The  major  difference  between  the  two 
t3q>es  of  diodes  is  in  the  magnitude  of  their 
parasitics.  In  the  p-n  junction  tunnel  diode, 
very  high  doping  densities  are  required  on  botli 
sides  of  the  junction  (N  ^  1  X  10^®  cm"^)  to 
achieve  a  high  tunneling  current  density.  This 
high  doping  density  creates  a  short  depletion 
layer  and  therefore  a  relatively  large  specific 
capacitance.  The  double-barrier  diode,  however, 
can  achieve  a  high  current  density  with  much 
lower  doping  densities  on  both  sides  of  the 
structure  (N  =  10^^  to  10^®  cm'®). 

A  good  figure  of  merit  for  comparing  the  speed 
of  different  diodes  is  the  ratio  of  specific  capac¬ 
itance  to  peak  current  density,  y  =  Cg/Jp.  This 
ratio  is  known  as  the  speed  index  (10);  it  is  a 
measure  of  the  current  available  for  charging 
the  device  capacitance.  For  Wafer  3,  this  quan¬ 
tity  is  7  =  1  X  10®pFcm'2/4.0  X  10“*  A  cm  ®  =  3.0 
psA^.  The  fastest  p-n  junction  tunnel  diodes 
ever  reported  (11)  were  made  of  GaAs  and  had 
7  =  14  to  16  psA^,  nearly  a  factor  of  five  poorer 
than  the  double-barrier  devices.  Moreover,  the 
p-n  junction  tunnel  diodes  ( 12)  achieved  a  max¬ 
imum  experimental  oscillation  frequency  of 
only  103  GHz  with  an  output  power  well  below 
1  mW. 

SELF>OSCILLATlNG  MIXERS 

Near  the  NDR  region  in  an  1-V  curve,  a 
resonant-turmeling  diode's  dynamic  conduc¬ 
tance  varies  rapidly  as  a  function  of  the  applied 
signal  voltage.  This  feature  makes  the  resonant 
timneling  diode  an  efficient  mixer.  The  Fourier 
series  of  the  diode’s  dynamic  conductance  has 
large  components  at  the  oscillation  frequency 
ffl  and  its  first  few  harmonics,  particularly  the 
even  harmonics  when  the  I-V  curve  is  antis3nm- 
metric  about  the  bias  point.  The  relative 
strength  of  these  components  will  determine 
the  efficiency  of  power  conversion  from  the 
signal  frequency  fg  to  the  intermediate  fre¬ 
quency  f],  assuming  that  fg  >  fj.  In  the  funda¬ 
mental  mode  of  conversion,  the  signal  has  fre¬ 
quency  fg  =  fo  ±  fp  and  in  the  second-harmonic 
mode,  the  signal  frequency  is  fg  =  2f^  ±  fj.  The 
most  efficient  conversion  should  be  in  the 
second-harmonic  mode  at  the  dc  bias  point  of 
maximum  negative  resistance,  because  this 
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point  has  approximate  antis3mimetry  in  the  I-V 
curve  (ie,  I(V)  =  -I(-V)].  However,  the  fundamen¬ 
tal  mode  should  achieve  its  maximum  efficien¬ 
cy  at  bias  voltages  nearest  the  regions  of  great¬ 
est  ciuvatiu^,  where  the  Fourier  series  of  g(t) 
has  a  predominant  coefficient  at  the  oscilla¬ 
tion  frequency. 

The  resonant-tunneling  self-oscillating  mixer 
has  the  potential  to  displace  the  Schottky  diode 
in  many  millimeter-wave  applications.  To  do 
this,  it  must  demonstrate  a  competitive  noise 
figure,  roughly  3  to  6  dB  in  the  microwave 
region  and  6  to  10  dB  in  the  millimeter  band. 
Although  noise  figru'es  have  not  yet  been  stud¬ 
ied,  measurements  of  stable  resonant-tunneling 
diodes  indicate  that  they  have  very  low  intrin¬ 
sic  noise  [  13|.  In  fact,  the  measured  noise  power 
is  less  than  expected  from  the  full  shot  noise 


expression,  for  reasons  that  are  not  fully 
understood.  Nevertheless,  this  low  noise  figure 
will  undoubtedly  work  to  the  benefit  of  the 
resonant-tunneling-diode-based  self-oscillat¬ 
ing  mixer.  And  the  resonant-tunneling  self- 
oscillating  mixer  provides  the  intrinsic  capa¬ 
bility  to  achieve  conversion  gain.  Mixers  such 
as  the  standard  Schottky  diode  aren't  capable 
of  gain  and  usually  show  losses  of  several 
decibels. 

Using  the  fimdamental  mode  of  coaxial  res¬ 
onant  tunneling  oscillators,  we  have  operated 
them  as  mixers  in  the  microwave  region.  Shown 
in  Fig.  7  is  the  mixing  spectrum  of  an  oscillator 
with  an  oscillation  frequency  of  14.2  GHz,  a 
signal  frequency  of  8.2  GHz,  and  an  interme¬ 
diate  frequency  of  6.0  GHz.  The  single-sideband 
conversion  loss  in  this  case  is  just  under  6  dB, 
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Fig.  7  —  The  frequency  down-conversion  of  a  signal  at  frequency  y  by  a  double-barrier  diode  oscillating  at  the 
frequency  marked  QW  produces  the  difference  frequency  labeled  QW-y. 
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which  is  only  about  3  dB  higher  than  the  best 
results  for  room-temperature  Schottky-diode 
mixers  in  the  same  frequency  region.  The  exper¬ 
iments  on  the  second-harmonic  mode  have  just 
begun,  using  a  waveguide  oscillator  in  the  mil¬ 
limeter  band.  The  best  conversion  loss  achieved 
to  date  is  about  12  dB  for  an  oscillation  fre¬ 
quency  of  50  GHz  and  a  signal  frequency  of  100 
GHz.  This  performance  is  comparable  to  the 
best  results  reported  for  a  pair  of  Schottky 
diodes  in  the  antiparallel  configuration  [14]. 

resistive  multipliers 

Resonant-tunneling  devices  excel  as  resis¬ 
tive  multipliers.  A  resistive  midtiplier,  which  is 
a  form  of  harmonic  mviltiplier,  generates  power 
at  frequencies  that  aren’t  conveniently  avail¬ 
able  from  fundamental  oscillators.  Harmonic 
multipliers  are  often  used  for  radio  astronomy. 
In  radio  astronomy,  heterodyne  receivers,  oper¬ 
ating  at  frequencies  above  100  GHz,  use  power 
from  a  harmonic  multiplier  as  alocal  oscillator 
source.  Harmonic  multipliers  are  also  the  pri¬ 
mary  source  of  power  for  molecular  spectros¬ 
copy  in  the  submillimeter  wavelength  spec¬ 
trum. 

Ciurently,  resistive  multipliers  usually  use 
Schottlty-barrier  diodes,  but  the  advantages  of 
resonant-tunneling  diodes  will  promote  their 
use  as  resistive  multipliers.  The  absence  of 
even  harmonics  in  aresonant-tunneling diode’s 
dynamic  conductance  simplifies  the  design  of 
multiplier  circuits,  particularly  in  the  milli¬ 
meter-wave  region.  Also,  because  of  its  nega¬ 
tive  resistance,  the  maximum  theoretical  gener¬ 
ation  efficiency  of  a  double-barrier  diode  can  be 
significantly  higher  than  an  ideal  diode’s  n*^ 
(where  n  is  the  harmonic  number)  [15]. 

If  a  resonant-tunneling  diode  is  pumped  so 
that  the  peak  amplitude  of  the  voltage  across 
the  diode  occurs  above  the  device’s  resonant- 
current  peak,  at  least  three  local  maxima  in  the 
diode’s  current  waveform  will  be  produced  over 
one  cycle.  (See  Fig.  8  for  a  theoretical  curve.) 
These  maxima  correspond  to  third-  or  higher- 
order  odd-harmonic  generation.  The  ease  with 
which  these  higher-order  odd  harmonics  are 
obtained  is  due  to  the  presence  of  a  peak  and  a 
valley  in  the  diode’s  I-V  curve,  as  well  as  the 
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Fig.  8  —  Voltage  and  current  waveforms  for  a  double¬ 
barrier  diode  reveal  that  the  local  maxima  in  the  current 
waveform,  if  equally  spaced  in  phase,  should  lead  to  a 
strong  fifth  harmonic  component  in  the  current-power 
spectrum. 

overall  antis5nnmetry  of  this  curve  about  the 
origin. 

To  demonstrate  this  odd-harmonic  genera¬ 
tion,  a  diode’s  I-V  curve  can  be  modeled  with  a 
seventh-order  pol3momial.  Using  this  model, 
the  diode  is  assumed  to  be  driven  by  a  source 
that  has  an  internal  impedance  less  than  the 
minimum  negative  resistance  of  the  double- 
barrier  diode.  Thus  the  current  across  the  diode 
is  a  single-valued  function  of  the  drive  voltage. 
The  numerically  determined  current-waveform 
of  cm  0.5  V  source  is  shown  in  Fig.  8.  The  local 
maxima  in  this  waveform,  if  equally  spaced  in 
phase,  should  lead  to  a  strong  fifth  harmonic 
component  in  the  current-power  spectrum. 
This  spectrum,  shown  in  Fig.  9,  predicts  that 
the  available  fifth-harmonic  power  will  be  0.065 
times  that  of  the  fundamental. 

Figure  10  illustrates  the  power  spectrum 
obtained  for  a  resonant-tunneling  diode  that 
was  mounted  in  a  50  O  coaxial  circuit  and 
pumped  at  4  GHz.  As  expected,  even  harmonics 
are  absent  and  the  fifth  harmonic  predominates 
among  the  odd  harmonics.  However,  the  mea- 
siired  efficiency  was  only  about  0.5%,  signifi¬ 
cantly  less  than  the  theoretical  prediction  of 
6.5%.  This  discrepancy  can  possibly  be  attrib¬ 
uted  to  the  test  circuit,  which  does  not  allow 
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Fig.  9  —  A  power  spectrum  corresponding  to  the  volt¬ 
age-current  plot  of  the  Fig.  6. 
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Fig.  10  —  A  strong  fifth-harmonic  component  was  con¬ 
firmed  experimentally  for  a  4.25-GHz  input  signal. 
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Fig.  11  — Persistent  photoconductivity  is  a  characteris¬ 
tic  of  resonant-tunneling  diodes.  The  experimental  re¬ 
sults  were  obtained  from  a  diode  that  was  cooled  to  20  K 
in  the  dark  and  then  exposed  to  progressively  greater 
levels  of  light. 


independent  tuning  of  the  harmonics.  Ideally, 
the  fifth  harmonic  should  be  terminated  with  a 
resistance  greater  than  the  sotirce  resistance. 
But  the  third  harmonic  should  be  terminated 
with  a  reactance,  which  eliminates  power  dis¬ 
sipation  and  reinforces  the  fifth-harmonic  vol¬ 
tage.  Combining  these  features  in  a  one  circuit 
will  allow  independent  tuning  of  the  harmonics 
and  greater  conversion  efficiency  for  the  fifth 
harmonic. 


•I 
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PERSISTENT  PHOTOCONDUCTIVITir 
AND  A  RESONANT-TUNNELING 
TRANSISTOR 

Atransistorbasedon  adouble-barrier  resonant¬ 
tunneling  structure  could  significantly  improve 
upon  the  present  state  of  the  art.  Persistent 
photoconductivity  provides  clues  to  the  opera¬ 
tion  of  such  a  transistor,  wherein  the  well  region 
acts  as  the  device’s  control  electrode. 

A  resonant-tunneling  structure  sometimes 
exhibits  persistent  photoconductivity  [  161,  that 
is,  the  structure’s  I-V  curve  shifts  with  succes¬ 
sive  exposures  to  light.  Figure  11  illustrates 
this  shift  Photoionization  of  DX  centers  in  the 
AIq  3630  7  barriers  in  the  structure  causes  the 
shift  in  the  I-V  curve.  (DX  centers  are  lattice 
vacancies  associated  with  a  nearby  donor  atom. 
For  a  detailed  discussion  of  DX  centers,  see 
Nelson  1 17]  or  Lang  et  al  [  18].)  This  photoioni¬ 
zation  creates  a  dipole  region  with  fixed  posi¬ 
tive  charge  in  the  barriers  and  free  electrons 
outside  the  barriers.  Any  electrons  that  reside 
initially  in  the  well  quickly  escape  because  of 
the  higher  energy  of  the  confined  band  edge 
within  the  well  and  the  short  resonant-state 
lifetime.  The  dipoles  produce  the  band  bending 
shown  in  Fig.  12.  The  well  is  lowered  in  poten¬ 
tial  by  an  amount  A(/>,  and  the  barriers  are  also 
lowered.  The  resvilting  shifted  I-V  curve  shows 


a  peak  transmission  that  has  been  moved  to  a 
lower  voltage  by  the  lowered  weU  and  a  higher 
ciurent  density  caused  by  the  lower  effective 
barriers. 

Figure  12  shows  that  DX  centers  in  the  barri¬ 
ers  provide  a  way  of  adjusting  the  well  potential 
relative  to  the  outer  contact  regions.  This  effect 
is  analogous  to  the  modulation  of  the  base 
potential  of  a  conventional  transistor.  In  this 
case,  the  structure’s  potential  well  serves  as 
the  control  electrode. 

The  transconductance  of  such  a  device,  Gn,, 
is  related  to  the  change  in  collector  current 
density,  AJ^.,  and  the  change  in  well  potential, 
A<f>,  by 

die  A.iJc 


where  A  is  the  emitter  area.  The  shift  in  the  well 
potential  is  approximately  related  to  the  shift 
in  the  current  peak  AVp  by 

AVp  -  4A<t,  (4) 

This  expression  is  valid  only  for  the  specific 
diode  of  SoUner  et  al  1 19].  One  factor  of  two  on 
the  right  side  approximately  accounts  for  the 
depletion  region  on  the  collector  side  when 
biased  to  the  ciurent  peak.  The  other  factor  of 
two  arises  from  the  base-to-emitter  voltage. 


Fig.  12  —  Positive  charge  in  the  barriers  of  a  double-barrier  diode  produces  band  bending  (red  line);  a  device  without 
ionized  centers  has  no  band  bending  (black  line). 
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which  is,  at  most,  half  of  that  applied  between 
the  base  and  collector.  Thus  the  transconduc¬ 
tance  becomes 

4A^c 

'5) 

But  tliere  is  one  complication.  The  leakage 
ciurent  also  changes  with  A<f>,  but  this  change 
is  probably  related  to  the  presence  of  DX  cen¬ 
ters  in  the  barriers  rather  than  to  the  modula¬ 
tion  of  the  well  potential.  To  accoimt  for  this 
contribution,  we  assume  that,  at  voltages  below 
the  current  peak,  the  transistor  output  imped¬ 
ance  is  large  and  positive,  so  we  can  extrapolate 
the  excess  ciurent  at  high  voltage  back  to  zero 
and  keep  the  output  impedance  positive.  These 
extrapolations  are  shown  in  Fig.  13.  Using 
this  figure,  an  emitter  width  of  1  /txm  gives  an 
intrinsic  transconductance  of  4,000  mS/mm. 


Specific  contact  resistances  of  the  order  of 
lO''^  n  cm2  are  necessary  to  realize  such  a  high 
extrinsic  transconductance. 

The  frequency  above  which  the  small-signal 
short-circuit  ciurent  gain  drops  below  unity  is 


where  C  is  the  total  input  capacitance.  Using  a 
l-/im  gate  contact  and  40-nm  n'  regions,  we  find 
fT=  110  GHz. 

The  maximum  frequency  of  oscillation  f^^^ 
will  also  depend  on  the  base  resistance.  For 
negligible  emitter  resistance, 

where  Rb  is  the  base  resistance  and  is  the 
base-to-coUector  capacitance.  As  the  base  is 
very  thin  and  lightly  doped,  it  might  be  expected 


Symbols 


Units  are  shown  in  brackets. 

A  Cross-sectional  area  of  a  device  Icm2) 
C  Ci^pacitance  (F] 

Cc  Collector  capacitance  [FJ 

Cs  Cs^jacitancc  per  unit  area  [F/cm^] 

Ef  Fermi  energy  [J] 

e  Charge  on  the  electron  [C) 

fi  Intermediate  frequency  [Hz] 

fmax  Frequency  above  which  the  terminal 

conductance  is  positive  for  all  vol¬ 
tages,  hence  maximum  oscillation 
"frequency  [Hz] 
fQ  Output  frequency  (HzJ 

fs  Sig^  frequency  [Hz] 

fx  Frequency  above  which  the  small- 

sigiLEd  short-circuit  ciurent  gain  is 
below  unity  [Hz] 

f  AE  Frequency  corresponding  to  resonant- 

state  lifetime  [Hz] 

Gjoax  Maximiun  negative  value  of  voltage- 
dependent  conductance  [S] 

G(V)  Voltage-dependent  conductance  of  a 
resonant-tunneling  diode  [S] 
g(t)  Conductance  of  a  diode  modulated  by 
a  changing  voltage  [S] 

I  Current  [A] 

IF  Intermediate  frequency  [Hz] 

Jc  CoUector  current  density  [A/cafl] 

Jp  Peak  current  density  [A/cm^] 


no  Carrier  concentration  [cm-3] 

Pinax  Maximum  output  power  obainable 

from  a  negative-resistance  amplifier 
[W] 

Rt)  Base  resistance  [fl] 

Rj  Current  responsivity  [A/W] 

Rs  Series  resistance  [fl] 

V  Voltage  [V] 

Vq  Voltage  amplitude  [V] 

Vm  Voltage  at  frequency  m  [V] 

Za  Anteima  impedance  [fl] 

Zo  Device  impcxlance  [fl] 

Z^  Load  impedance  [fl] 

AE  Energy  half-width  of  transmission 
function  [J] 

AI  Current  excursion  of  negative  differ¬ 

ential  resistance  region  [A] 

AE  Change  of  current  [A] 

APac  Change  of  ac  power  [W] 

AVp  Shift  of  the  voltage  of  the  ciurent 

peak  from  ionized  charges  in  the  bar¬ 
riers  [V] 

AV  Voltage  excursion  of  negative  dif¬ 
ferential  resistance  region  [A] 

A<t>  Shift  in  potential  of  the  well  from 
ionized  charges  in  the  barriers  [eV] 
c  Dielectric  constant  [F/cm] 

7  Speed  index  [s/V] 

w  Angular  frequency  [s’l] 

<^(x)  Electrostatic  potential  energy  [eV] 
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that  the  base  resistance  would  be  high.  How¬ 
ever,  the  quantum-well  geometry  reduces  the 
resistance  considerably.  Electron  mobilities  in 
modulation-doped  quantum  wells  are  very  high 
at  temperatures  less  than  100  K.  Although  this 
structime  is  not  a  modulation-doped  well  in  the 
usual  sense,  the  base  layer  near  the  resommt 
current  peak  acquires  charge  from  the  transit¬ 
ing  carriers  without  any  direct  doping  [20].  For 
present  samples,  the  average  well-charge  den¬ 
sity  is  approximately  10^^  cm'®. 

Using  a  mobility  of  5  X  10®  cm^/V-s  for  a 
base  carrier  density  2  X  10*^'^  cm*®  (the  best 
mobility  [21]  at  50  K),  a  1  X  20-Mm  emitter,  and 
a  4.0-nm  well  thickness,  the  base  resistance  is 
less  than  10  H.  This  gives  an  product  of 
about  0.5  ps,  resulting  in  an  f,nnv  of  330  GHz. 
Any  increase  in  current  density  increases  fj 
and  fn,av  proportionately.  Theoretically  pre¬ 
dicted  ciurents  are  several  times  higher  than 
presently  observed  values. 


CONCLUSIONS 

Much  fertile,  yet  imcultivated,  territory  re¬ 
mains  in  the  realm  of  multiple-barrier  tunnel¬ 
ing.  Three-terminal  devices  are  in  their  infancy, 
as  are  devices  for  digital  applications.  The  pos¬ 
sibility  of  electro-optical  effects  has  received 
some  attention,  and  several  novel  devices  have 
resulted  [22] .  Most  structures  to  date  have  con¬ 
tained  only  two  barriers,  but  investigation  of 
the  interaction  of  three  or  more  barriers  may 
yield  even  more  interesting  and  useful  phenom¬ 
ena. 

An  extension  of  multiple-barrier  structures 
to  the  inclusion  of  many  barriers  produces  a 
superlattice.  Because  of  domain  formation, 
superlattices  have  been  impractical  electrical 
devices  in  the  past.  Many  variations  remain 
untested,  even  vmanalyzed.  For  example,  the 
“CHIRP”  superlattice  was  examined  in  some 
detail  [23] ,  but  only  one  type  of  period  variation 
was  investigated.  The  domain  formation  that 
has  limited  this  field  to  date  may  be  overcome 
with  structures  that  are  yet  undeveloped. 

All  of  the  devices  discussed  in  this  paper 
have  been  small  compared  with  their  operating 


Fig.  13  —  These  calculated  l-V  curves  for  a  resonant¬ 
tunneling  transistor  are  based  upon  the  assumption  that 
the  output  impedance  of  the  transistor  is  large  and  posi¬ 
tive  and  that  the  excess  current  at  high  voltage  can  be 
extrapolated  smoothly  back  to  zero. 


-2-10  1  2 
Voltage  (V) 


Fig.  14  —  The  predicted  peak-to-valley  ratio  of  a  reso¬ 
nant-tunneling  device's  l-V  curve  is  much  greater  than 
that  measured  for  any  experimental  devices.  Improve¬ 
ments  in  material  purity  should  bring  the  observed  val¬ 
ues  into  closer  agreement  with  the  theory. 
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wavelength.  Therefore,  they  are  lumped  ele¬ 
ments,  not  transmission  media.  A  bulk  mate¬ 
rial  that  could  provide  large,  fast  impedance 
nonlinearities  would  be  useful  at  high  frequen¬ 
cies.  With  two-barrier  structures,  these  nonlin¬ 
earities  are  difficult  to  create,  but  some  variety 
of  superlattice  may  provide  the  material  neces¬ 
sary  for  large-scale  structures  that  will  find 
application  as  traveling-wave  devices. 

The  materials  growth  that  makes  possible 
the  high  spatial  resolution  in  heterostructiu'e 
formation  is  well  advanced  and  progressing 
rapidly,  but  much  improvement  is  still  possi¬ 
ble.  Material  of  sufficient  quality  for  observing 
reeisonable  peak-to-valley  ratios  at  room  tem¬ 


perature  has  been  available  only  for  two  years. 
Theoretically,  as  shown  in  Fig.  14,  this  p)eak-to- 
valley  ratio  should  be  much  greater.  Materials 
of  greater  purity  could  validate  most  of  the 
assumptions  leading  to  the  calciilated  curve  of 
Fig.  14. 

In  short,  the  future  of  multiple-bcirrier  tun¬ 
neling  structures  appears  bright,  especially  in 
view  of  the  high  level  of  interest  in  the  field. 
Some  of  the  devices  that  are  just  beginning  to 
demonstrate  feasibility  will  be  commercially 
exploited  in  the  next  several  years  and  many 
yet  undevised  structures  will  reach  the  same 
point  in  the  years  to  come. 
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Appendix:  Resonant  Tunneling  Theory 


The  theory  of  resonant  tunneling  is  still  evolving. 
The  full  problem  is  complicated,  involving  self- 
consistent  solution  of  Poisson's  equation  and  the 
time-dependent  Schrodinger  equation.  The  spatial 
quantization  that  creates  the  effects  of  interest  also 
renders  inadequate  many  of  the  conunon  approxi¬ 
mations  used  in  the  solution  of  these  equations.  To 
further  complicate  matters,  scattering  from  impuri¬ 
ties.  defects,  other  carriers,  and  collective  excita¬ 
tions  are  almost  certainly  important,  but  the  inclu¬ 
sion  of  these  effects  in  detail  is  a  difficult  theoretical 
problem.  On  the  positive  side,  however,  the  simplest 
calculations  seem  to  explain  the  general  outline  of 
the  experimental  observations. 


Stationary-State  Calculation 

The  first  simplification  to  make  in  the  analysis  of 
the  resonant-tunneling  structure  is  the  elimination 
of  time-dependence  in  the  solution  of  the  control¬ 
ling  equations.  This  stationary-state  treatment  is 
based  on  the  approach  first  applied  to  finite  super¬ 
lattices  by  Tsu  and  Esaki  [  1 ).  The  calculation  shown 
here  incorporates  some  of  our  own  refinements. 

For  a  first-order  solution  to  the  problem,  we 
assume  that  electrons  are  the  only  charge  carriers 
and  that  they  interact  only  with  potential  discontin¬ 
uities  in  the  conduction  band.  In  addition,  we  use 
the  simplification  of  the  effective-mass  model  for 
the  electrons.  The  matching  of  wave-functions  at 
hetero-interfaces  is  based  on  the  conservation  of 
probability  current,  which  also  accounts  for  the  dif¬ 
ferent  effective  masses  in  different  semiconductors. 
We  use  superposition  to  determine  the  potential 
profile  as  an  electron  makes  a  transit  of  the  struc¬ 
ture;  the  potential  is  taken  as  the  sum  of  the  poten¬ 
tial  due  to  the  conduction-band  offsets,  the  applied 
bias,  and  any  localized  charges  within  the  structure. 
A  further,  and  final,  simplifying  assumption  is  that 
the  regions  outside  the  barriers  house  an  electron 
gas,  the  field-screening  behavior  of  which  is  des¬ 
cribed  by  the  Thomas-Fermi  approximation. 

The  first  step  in  the  calculation  of  a  resonant¬ 
tunneling  I-V  curve  is  the  determination  of  the  spa¬ 
tial  variation  of  the  electrical  potential  experienced 
by  the  electrons  as  they  transit  the  double-barrier 
device  under  consideration.  In  addition  to  the  poten¬ 
tial  change  produced  by  the  barrier  regions,  the  elec¬ 
trostatic  potential  applied  by  the  external  (biasing) 
field  produces  accumulation  and  depletion  regions 
that  will  perturb  the  potential  profile.  The  spatial 
dependence  of  the  electrostatic  potential,  </>(x),  is 


calculated  in  the  barrier  region  using  the  self- 
consistent  Thomas-Fermi  equation  (egs  units). 


<^"(x)  = 


4neno 

< 


(A-1) 


where  is  the  Fermi  energy,  c  is  the  dielectric  con¬ 
stant,  and  nQ  is  the  average  volume-doping  density. 
This  expression  is  solved  analytically  in  the  regions 
with  nonzero  net  charge  and  then  matched  numeri¬ 
cally  at  the  boundaries.  In  practice,  a  parabola  is 
fitted  to  the  solution  in  the  cathode  to  limit  the 
band-bending  region.  The  spatial  extent  of  the  fitted 
region  is  taken  to  be  about  twice  the  Thomas-Fermi 
screening  length.  A  potential  profile  calculated  for  a 
double-barrier  structure  with  doping  carrier  density 
of  10^^  cm'^  in  the  outer  regions  is  shown  in  the 
Figure. 

Once  the  potential  has  been  found,  the  transmis¬ 
sion  coefficient  of  an  electron  with  a  given  energy  is 
determined  by  using  the  transfer-matrix  method  and 
matching  boundary  conditions.  The  ciurent  is  then 
found  by  integrating  over  the  incident  electron  dis¬ 
tribution  and  the  density  of  available  final  states  1 1 ). 
The  transfer  matrix  is  numerically  calculated  for  an 
arbitrary  parabolic  potential  profile.  To  guarantee 
rapid  convergence,  the  potential  is  divided  into  a  few 
(<  10)  spatial  segments  with  transfer  matrices  cal¬ 
culated  for  each  partition. 

Even  though  the  cathode  contains  an  accumula¬ 
tion  region,  we  obtain  better  agreement  between  the 
calculated  and  measured  transmission  coefficients 
if  we  use  the  equilibrium  Fermi  distribution  of  the 
cathode  material  and  ignore  possible  tuimeling  of 
the  accumulated  electrons.  This  agreement  between 
theory  and  experiment  may  be  due  to  the  improba¬ 
bility  of  electrons  scattering  (inelastically)  into  the 
accumulation  region,  which  makes  the  lifetime  for 
timneling  out  of  the  region  shorter  than  the  time  for 
scattering  in.  However,  this  means  that  there  is 
much  less  band  bending  than  has  been  assumed. 


Temporal  Behavior 

Time  evolution  of  the  S3^tem  has  been  eliminated 
from  the  stationary-state  calculation  above.  But  it  is 
possible  to  find  the  lifetime  of  a  transiting  electron 
in  the  resonant  state  that  it  occupies  while  it  is 
between  the  two  barriers.  The  approach,  first  derived 
for  atomic  and  nuclear  resonances,  is  applicable  to 
any  resonant  state  and  can  be  applied  to  resonant 
tunneling.  The  lifetime  is  given  by  the  uncertainty 
relation,  t  =  >6/ AE;  AE  is  the  half-width  of  the  peak  In 
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the  transmission  coefficient  T(E).  Frequencies  cor¬ 
responding  to  =  (27rr)‘^  are  calculated  for  the 
three  structures  and  listed  in  Table  1. 

A  time-dependent  calculation  of  the  electron  dis¬ 
tribution  has  been  performed  by  using  the  Wigner 
function  approach  [20].  The  time  required  to  reach 
equilibrium  after  an  impulse  to  the  system  is  com¬ 
parable  to  that  estimated  by  the  resonant-state  life¬ 
time.  This  calculation,  however,  offers  the  advantage 
of  predicting  the  time  evolution  of  the  wave  func¬ 
tions  for  all  times. 


Scattering 

A  complete  time-dependent  solution  should  take 
account  of  charge  redistribution  from  applied  elec¬ 
tric  fields  and  motion  of  carriers  through  the  deple¬ 
tion  and  accumulation  regions.  Scattering  plays  an 
important  part  in  these  processes,  because  scatter¬ 
ing  in  the  barriers  removes  the  assumed  condition 
of  conservation  of  transverse  momentum  during  tun¬ 
neling,  and  scattering  in  the  well  region  or  regions 
destroys  the  coherence  of  a  carrier  wave  function 
over  the  structure. 

According  to  Lxuyi  [24],  scattering  may  dominate 
the  charge-transport  process  through  a  double¬ 
barrier  structure.  His  proposal  has  led  to  the  devel¬ 
opment  of  a  second,  sequential,  model  of  the  double¬ 
barrier  structure.  In  contrast  to  the  ballistic  model. 


0  20  40  60  80  100  120  140 

Distance  (nm) 


Calculated  potential  energy  of  an  electron  as  it  transits  a 
double-barrier  diode. 


which  maintains  the  coherence  of  the  wave  function 
across  the  incoming  barrier,  the  sequential  model 
attempts  to  incorporate  the  effects  of  scattering  by 
eliminating  the  assumption  of  phase  coherence.  This 
model  is  called  the  sequential  model  because  each 
electron  is  presumed  to  enter  the  structure  and 
tunnel  through  the  first  barrier,  achieve  a  resonant 
state  within  the  structure’s  potential  well,  and  then 
exit  the  structure  sequentially  by  tunneling  through 
the  remaining  barrier.  Luryi  proposed  that  all  phase; 
coherence  is  lost  before  each  electron  leaves  the; 
structure’s  potential  well  and  that  the  negative; 
resistance  exhibited  by  the  double-barrier  device; 
results  from  the  conservation  of  momentum.  When 
the  peak  in  the  density  of  states  in  the  well  (for 
momentum  perpendicular  to  the  well)  is  below  the 
conduction  band  edge  of  the  cathode,  electrons  can 
not  tunnel  into  the  well. 

Weil  and  Vinter  [25]  have  recently  calculated  thal. 
the  peak  current  densities  should  be  the  same  for 
both  the  ballistic  and  the  sequential  models.  It  may 
be,  however,  that  the  two  models  are  equivalent  as 
they  are  presently  defined.  In  their  sequential- 
tunneling  calculations,  Weil  and  Vinter  use  the  same 
boundary  conditions  as  those  used  for  the  resonant 
model  (see  Ref.  8  in  their  paper);  they  maintain  the 
coherence  of  the  wave  function  across  the  incoming 
barrier.  The  effects  of  the  second  barrier  enter  via 
the  width  of  the  quasibound  state  in  the  well,  again 
with  a  calculation  that  assumes  coherence  across 
the  barrier.  Perhaps  it  is  not  surprising  that,  with  no 
scattering  in  the  well,  the  calculations  of  ballistic 
and  sequential  tunneling  give  the  same  result. 

Weil  and  Vinter  do  attempt  to  include  the  effects  of 
scattering  in  the  well,  via  a  broadening  of  the  quasi¬ 
bound  state.  They  conclude  that,  as  long  as  the 
energy  width  of  this  state  is  small  compared  to  the 
energy  width  of  the  incoming  carriers  (the  usual 
experimental  situation),  the  lower  peak  transmis¬ 
sion  is  offset  by  the  additional  width.  Thus  the  inte¬ 
gral  determining  the  peak  ciurent  is  unaffected. 
However,  the  technique  of  including  scattering  still 
implicitly  assumes  the  phase  coherence  discussed 
above,  even  though  scattering  will  generally  alter  the 
phase  after  the  scattering  event. 

This  calculation  does  not  appear  to  be  the  same  as 
the  one  originally  suggested  by  Liuyi  [24]  because  it 
does  not  perturb  the  phase  coherence  across  the 
structure.  There  are  also  important  scattering  pro¬ 
cesses,  such  as  scattering  in  the  barriers,  that  are 
not  in  Luryi’s  model.  Thus  the  role  of  scattering 
remains  one  of  several  important  imanswered  ques¬ 
tions  in  resonant  tunneling. 
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APPENDIX  B 

Large  room-temperature  effects  from  resonant  tunneling  through  AlAs 
barriers 

W.  D.  Goodhue,  T.  C.  L  G.  Sollner,  H.  Q.  Le,  E.  R.  Brown,  and  B.  A.  Vojak** 

Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  Massachusetts  02173 

(Received  3  June  1986;  accepted  for  publication  2  September  1986) 

At  room  temperature,  we  have  observed  negative  differential  resistance  in  AlAs  double-barrier 
structures  and  a  large  hysteresis  in  the  current-voltage  characteristic  of  a  stack  of  five  AlAs 
double-barrier  structures.  The  peak-to-valley  ratio  of  the  current  was  as  high  as  3.5:1  in  a 
double-barrier  structure.  To  the  best  of  our  knowledge,  this  is  the  largest  room-temperature 
peak-to-valley  ratio  observed  to  date  in  a  double-barrier  structure  and  the  first  report  of  a 
room-temperature  hysteresis  in  a  stacked  structure.  These  structures  were  grown  by  molecular 
beam  epitaxy  using  thin  AlAs  barriers  in  GaAs.  Both  the  first  and  second  resonances  were 
observed,  and  are  well  explained  by  simple  tunneling  theory  assuming  a  value  of  1.0  ±  0. 1  eV 
for  the  GaAs-AlAs  conduction-band  discontinuity  seen  by  the  tunneling  electrons.  This  value 
is  very  close  to  the  difference  in  conduction-band  energy  at  the  T  points  found  by  using  the 
accepted  values  of  GaAs  and  AlAs  band  gaps  with  65%  of  the  band-gap  difference  appearing 
in  the  conduction  band.  This  suggests  that  negligibly  few  electrons  relax  to  the  lower  AlAs  A' 
valley  as  they  tunnel  through  the  1.5-2.5-nm-thick  AlAs  barriers.  These  results  indicate  that 
AlAs  should  be  a  high  quality  barrier  material  for  a  variety  of  heterojunction  devices. 


Resonant  tunneling  through  double-barrier  semicon¬ 
ductor  structures  has  received  much  interest  in  the  past  few 
years.  These  structures  consist  of  two  semiconductor  bar¬ 
riers  roughly  10  nm  or  less  thick,  embedded  in  a  semicon¬ 
ductor  having  a  smaller  band-gap  energy  than  that  of  the 
barrier  material,  and  separated  by  less  than  about  10  nm.  In 
this  way  the  carriers  are  confined  to  the  lower  band-gap  ma¬ 
terial  but  can  readily  tunnel  through  the  thin  barriers  under 
applied  bias  voltage.  In  1970Tsu  and  Esaki*  first  proposed 
that  a  large  peak  in  the  tunneling  current  through  the  struc¬ 
tures  should  occur  when  the  injected  carriers  have  certain 
resonant  energies.  Both  Chang  et  al?  and  Vojak  et  al?  ob¬ 
served  slight  negative  resistance  from  single  quantum  well 
and  multilayer  heterostructures  grown  by  molecular  beam 
epitaxy  (MBE)  and  metalorganic  chemical  vapor  deposi¬ 
tion  (MOeVD),  respectively.  In  1983  Sollner  cra/.^  report¬ 
ed  the  first  observations  of  resonant  tunneling  at  room  tem¬ 
perature  and  of  large  regions  of  negative  differential 
resistance  at  lower  temperatures.  Their  experiments  indicat¬ 
ed  that  the  intrinsic  response  time  is  shorter  than  1 X  lO”  s. 
Sollner  et  al?  also  reported  microwave  oscillations  at  low 
temperatures.  Recently,  Bonnefoi  et  al?  have  observed  reso¬ 
nant  tunneling  in  a  MOCVD-grown  GaAs/AlAs  structure, 
Shewchuk  et  al?  have  observed  room-temperature  reso¬ 
nant-tunneling  oscillations  and  negative  differential  resis¬ 
tance  in  MBE-grown  GaAs/GaAlAs  structures,  and  Tsu- 
chiya  et  al?  have  observed  room-temperature  negative 
differential  resistance  from  MBE-grown  GaAs/ AlAs  struc¬ 
tures. 

In  this  letter  we  report  a  variety  of  resonant  tunneling 
effects  observed  in  GaAs  quantum  well  structures  with  AlAs 
barriers.  All  of  the  structures  were  grown  by  MBE  on  ( 100) 
oriented  n  ^  -GaAs  substrates  and  utilized  at  least  a  l-/im- 
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thick  temperature  stabilizing  GaAs  buffer  layer.  The  growth 
temperature  was  560  *C,  and  the  growth  rates  were  about  1 .0 
/im/h  and  0.2  /im/h  for  the  GaAs  and  AlAs,  respectively. 
Computer-controlled  shutters  were  used  during  the  growth 
of  the  barriers  and  wells  of  the  samples.  The  wells  and  bar¬ 
riers  were  unintentionally  doped p  type  with  a  net  concentra¬ 
tion  5x10*^  cm“^,  while  the  remaining  GaAs 

was  doped  n  type  with  silicon  to  a  concentration 
Njy  —  A^'^lXlO*^  cm“^.  Low  resistance  contacts  were 
obtained  by  evaporating  and  alloying  a  standard  Ni/Ge/Au 
metallization,  and  individual  mesas  were  defined  on  the  wa¬ 
fers  by  ion  beam  etching.  Further  details  of  the  fabrication 
procedure  are  described  elsewhere.'* 

Figure  1  shows  the  curves  of  room-temperature  current 
density  versus  voltage  (J-V)  for  two-terminal  devices 
formed  from  both  a  single  double-barrier  structure  with  1.5- 
nm-thick  barriers  and  a  4.5-nm-thick  well,  and  a  stack  struc¬ 
ture  consisting  of  five  double-barrier  structures  with  1.5-nm 
barriers  and  a  4.5-nm  well,  separated  by  80  nm  of  /i-rype 
GaAs  having  a  donor  concentration  —  1  X  10*^ 

cm“^.  As  shown  in  Fig.  1(a),  the  double-barrier  structure 
exhibited  large  current  densities  and  high  peak-to-valley  ra¬ 
tios  at  room  temperature.  The  current  density  at  the  peak  of 
the  first  resonance  is  about  4 X  10^  A  cm"^  while  the  peak- 
to-valley  ratio  is  about  3.5:1.  At  77  K  the  current  density  of 
this  sample  remained  practically  the  same  as  above  but  the 
peak-to-valley  ratio  increased  to  about  10:1.  In  addition, 
both  the  first  and  second  resonances  are  observed.  As  shown 
in  Fig.  1  (b),  the  stack  also  exhibited  large  current  densities 
at  room  temperature  although  a  hysteresis  effect  was  ob¬ 
served.  This  hysteresis  effect  is  presently  under  further  in¬ 
vestigation.  The  transmission  peak  in  the  stack  was  approxi¬ 
mately  five  times  higher  in  voltage  than  for  a  single 
double-barrier  structure.  This  is  the  exj)ccted  result  if  the 
double-barrier  structures  in  the  stack  are  far  enough  apart 
that  the  electrons  can  relax  in  the  buffer  layers  between 
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(a)  VOLTAGE  (V) 


VOLTAGE  (V) 

FIG.  1.  ( a)  y-FcharacteristicofanAlAs  double-barrier  Structure  with  two 

1.5-nm  barriers  and  a  4.5-nm  GaAs  well.  Both  the  first  and  the  second  re¬ 
sonances  are  clearly  delineated  at  room  temperature.  The  chairlike  form  of 
the  y*  Feurve  in  the  negative  resistance  region  of  the  first  resonance  is  indi¬ 
cative  of  electrical  oscillation  with  the  measurement  circuit,  (b)  7-  K  charac¬ 
teristic  of  a  stack  of  five  AlAs  double  barriers,  each  similar  to  that  of  (a), 
with  80.0  nm  of /i-type  GaAs  between  the  pairs.  The  dashed  portion  of  the 
curve  indicates  switching  or  bistable  behavior.  The  arrows  are  located  on 
that  part  of  the  7-  Kcurve  which  is  traced  if  the  hysteretic  region  is  entered 
with  a  direction  defined  by  the  arrow. 

them.  Each  double  barrier  then  presents  the  same  tunneling 
transmission  coefficient  and  the  structures  behave  like  inde¬ 
pendent  scries  elements.  Further  consideration  of  the  stack 
structure  will  be  given  in  a  future  letter. 

We  grew  several  double-barrier  structures  having  bar¬ 
rier  thicknesses  between  1.5  and  2.5  nm  and  well  widths 
between  4.5  and  5.6  nm  for  analyzing  the  conduction-band 
discontinuity  between  AlAs  and  GaAs.  The  dimensions 
were  confirmed  to  within  10%  using  transmission  electron 
microscopy  (TEM).  These  devices  exhibited  resonances  at 
room  temperature  and  a  higher  voltage  resonance  became 
visible  at  low  temperatures.  Because  of  the  high  current  den- 


VOLTAGE  (V) 


FIG.  2.  /-  ^characteristics  at  77  K  of  an  AlAs  double-barrier  structure  with 
2.5-nm  barriers  and  a  4.5-nm  well.  The  second,  higher  voltage  resonance, 
which  appears  as  a  bump  at  room  temperaure,  is  clearly  delineated  at  77  K. 
The  energy  levels  of  the  well  were  calculated  from  theory  using  both  the  /-  V 
characteristic  and  TEM  measured  thicknesses. 

sides  generated  when  operating  the  devices  at  the  second, 
higher  voltage  resonance,  the  current-voltage  (/-lOcurves 
for  these  devices  were  generated  in  a  pulsed  mode.  Figure  2 
shows  the  /-  V  characteristic  at  77  K  of  a  device  with  2.5-nm 
AlAs  barriers  and  a  4.5-nm  GaAs  well.  Shown  in  Fig.  3  is  a 
TEM  cross-sectional  micrograph  of  this  device.  The  second 
resonance,  which  manifests  itself  only  as  a  small  change  of 
curvature  at  room  temperature,  becomes  a  clearly  delineated 
kink  in  the  7- F  curve  at  77  K. 

To  determine  the  conduction-band  discontinuity 
between  GaAs  and  AlAs,  we  analyzed  the  above  7- F  curves 
using  the  resonant  tunneling  theory  described  in  Ref  9.  Fig¬ 
ure  4  shows  one  important  result  of  this  theory:  the  position 


FIG.  3.  TEM  cross-scctional  micrograph  of  the  device  of  Fig.  2. 
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BARRIER  HEIGHT  (eV) 

FIG.  4.  Diagram  showing  how  sensitive  the  second,  higher  voltage  reson¬ 
ances  of  double-barrier  structures  are  to  barrier  height. 


of  the  second  resonance  depends  much  more  strongly  on  the 
barrier  height  than  does  the  position  of  the  first  resonance, 
which  depends  strongly  on  the  well  width.  This  allows  us  to 
use  the  position  of  the  first  resonance  to  fix  the  well  width 
and  the  position  of  the  second  to  determine  the  barrier 
height.  For  the  four  samples  measured  with  barrier  thick¬ 
nesses  between  1.5  and  2.5  nm  and  wells  between  4.5  and  5.6 
nm,  the  best  agreement  with  the  theoretical  /-K  curves  in 
each  case  required  a  1.0  ±  0.1  eV  conduction-band  discon¬ 
tinuity. 

Using  the  known  F-point  energy  gaps  of  1.424  eV  for 
GaAs  and  3.018  eV  for  AlAs^®  with  65%  of  the  band-gap 
discontinuity  in  the  conduction  band,  we  calculate  the  F- 
point  discontinuity  to  be  1.036  eV.  This  is  close  enough  to 
our  1.0  eV  estimate  to  strongly  suggest  that  direct  F-valley 
tunneling  processes  are  favored  over  indirect  processes  in  the 
thin  AlAs.  This  is  not  surprising  since  the  indirect  processes 
require  phonon  interactions  to  conserve  transverse  crystal 
momentum  or  other  scattering  processes  to  allow  tunnehng 
without  momentum  conservation.  Few  if  any  phonons  can 
be  created  or  absorbed  in  the  thin  AlAs  barriers  because  the 
time  that  the  electron  wave  packet  spends  in  these  regions 
during  tunneling  is  prohibitively  small.  On  the  other  hand, 
impurity  or  crystal-defect  assisted  tunneling  is  certainly 
present  and  is  thought  to  be  responsible  for  the  rapidly  in¬ 
creasing  current  in  the  region  of  the  second  resonance.  The 
uncertainty  in  our  barrier  height  determination  includes  the 
small  shift  of  the  second  resonance  peak  due  to  the  presence 
of  this  current  component.  Our  results  provide  independent 
support  of  the  65-35  rule”  over  the  older  85-15  rule‘s  per¬ 
taining  to  the  percentage  of  F-point  band-gap  discontinuity 
that  appears  across  the  conduction  and  valence  bands,  re¬ 
spectively. 

Another  interesting  effect,  shown  in  Fig.  5,  is  the  room- 
temperature  /-  V characteristic  of  an  asymmetric  double-bar¬ 
rier  AlAs  structure  with  a  2.5-nm  barrier,  4.0-nm  well  and 
5.0-nm  barrier.  The  asymmetry  observed  in  the  I-V curve  is 
consistent  with  the  theory.  It  is  apparent  that  many  different 
/-  V  curves  can  be  generated  for  a  variety  of  devices  by  vary¬ 
ing  double-barrier  parameters. 


FIG.  5.  Room-temperature  /- ^characteristic  of  an  asymmetric  AlAs  dou¬ 
ble-barrier  structure.  The  barrierthicknessesare2.5nm  and  5.0  nm,  respec¬ 
tively,  and  the  well  width  is  4.0  nm. 

In  summary,  we  have  observed  large  resonant  tunneling 
effects  at  room  temperature  in  AlAs  double-barrier  struc¬ 
tures  and  in  a  stack  of  five  AlAs  double-barrier  structures. 
The  peak  of  the  tunneling  current  is  raised  to  the  2-5  V  range 
in  the  stacked  structure,  thereby  demonstrating  a  means  of 
increasing  the  power  handling  capacity  of  devices  based  on 
this  effect.  From  the  observed  first  and  second  resonances  in 
the  tunneling  current,  we  have  inferred  a  1.0  ±  0.1  eV  F- 
valley  conduction-band  discontinuity  at  the  interface 
between  GaAs  and  thin  AlAs  barriers.  The  Jf-valley  and  L- 
valley  band  gaps  are  not  important  in  determining  the  bar¬ 
rier  height  for  resonant  tunneling  of  electrons  through  thin 
AlAs  barriers.  We  have  fabricated  an  asymmetric  AlAs  dou¬ 
ble-barrier  structure  and  observed  the  expected  asymmetric 
/-  V characteristic.  Our  resonant  tunneling  model  agrees  well 
with  both  I-V characteristics  and  TEM  thickness  measure¬ 
ments  of  the  structures.  Our  experiments  indicate  that 
GaAs-AlAs  double-barrier  structures  can  be  grown  with 
high  quality  interfaces  and  with  properties  that  appear  po¬ 
tentially  useful  for  a  variey  of  devices. 
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AtPENDiX  C 

Fundamental  oscillations  up  to  200  GHz  in  resonant  tunneling  diodes 
and  new  estimates  of  their  maximum  oscillation  frequency 
from  stationary-state  tunneling  theory 
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Fundamental  oscillations  have  been  measured  up  to  200  GHz  in  resonant-tunneling  diodes  at 
room  temperature.  Oscillations  in  the  range  102-112  GHz  were  achieved  with  diodes  mounted 
in  a  WR-6  waveguide  resonator,  and  the  peak  output  power  in  this  range  was  approximately  5 
The  same  diodes  oscillated  between  192  and  201  GHz  and  generated  about  0.2  //W  when 
mounted  in  a  WR-3  resonator.  The  estimated  maximum  oscillation  frequency  ( )  for  these 
devices  is  244  GHz,  assuming  the  average  drift  velocity  across  the  depletion  layer  to  be  4  X  10^ 
cm  s“ '.  This  estimate  has  been  obtained  from  a  new  phenomenological  theory  of  the  negative 
differential  conductance  which  accounts  for  the  frequency-dependent  spreading  resistance  and 
transit-time  delay.  The  theory  is  also  used  to  show  that  diodes  having exceeding  600  GHz 
are  feasible  simply  by  modifying  the  doping  profile  in  the  regions  on  either  side  of  the  double¬ 
barrier  structure. 


I.  INTRODUCTION 

Several  experimental  and  theoretical  studies  of  reso¬ 
nant-tunneling  oscillators  have  been  conducted  since  the 
first  microwave'  and  room-temperature^  oscillation  results 
were  obtained.  In  a  recent  letter  we  reported  the  first  milli¬ 
meter-band  oscillations  based  on  the  negative  differential 
conductance  (NDC)  of  double-barrier  diodes  at  room  tem¬ 
perature.^  We  also  emphasized  that  two  of  the  important 
requirements  for  fast  diodes,  high  NDC  and  low  specific 
capacitance,  would  be  achieved  by  using  thin  A1  As  barriers 
with  moderate  doping  outside  of  these  barriers.  AlAs  was 
the  material  of  choice  for  the  barriers  because  it  yielded  a 
large  room-temperature  peak-to-valley  ratio,  an  important 
requirement  for  large  NDC.  The  barriers  were  made  thin 
(1.7  nm)  to  compensate  for  the  large  F-point  AlAs-GaAs 
conduction-band  offset  (sl.O  eV)  and  to  allow  for  high 
peak  current  density/  Thin  barriers  also  decrease  the  life¬ 
time  of  the  tunneling  electron  in  the  quasibound  state  of  the 
region  (quantum  well)  between  the  barriers,  which  further 
increases  the  frequency  limit  imposed  by  this  lifetime.^  Our 
previous  results  led  to  the  prediction  that  the  fastest  device 
could  oscillate  up  to  186  GHz.  The  measurement  of  oscilla¬ 
tions  up  to  201  GHz,  as  reported  here,  has  compelled  us  to 
derive  a  new  estimate  of/n,,  using  the  more  sophisticated 
analysis  of  the  double-barrer  diode  described  here. 

II.  EXPERIMENT 

The  diodes  studied  here  were  taken  from  the  same  wafer 
that  yielded  the  highest  frequency  oscillations  in  our  pre¬ 
vious  work.^  In  brief,  these  diodes  consist  of  two  1.7-nm- 
thick  (three  lattice  constants),  undoped- AlAs  barriers  sepa¬ 
rated  by  a  4.5-nm-thick  quantum  well  made  of  undoped 
GaAs.  Outside  of  each  barrier  is  a  0.5-/zm-thick  layer  of  n- 
GaAs,  doped  to  2x10'^  cm“^.  The  entire  structure  was 
grown  by  molecular-beam  epitaxy  (MBE)  on  an  sub¬ 
strate  at  560 ’C.^  This  growth  temperature  is  lower  than 
those  usually  applied  in  MBE  and  was  used  here  to  suppress 


the  rapid  diffusion  of  A1  and  Si.  This  yields  sharper  inter¬ 
faces  and  cleaner  barrier  material,^  and  is  likely  the  reason 
for  the  large  room-temperature  peak-to-valley  ratio  (3.5:1 ) 
of  the  present  device.  Devices  were  fabricated  with  a  stan¬ 
dard  sequence  of  steps,  including  patterning  of  Ni/Ge/Au 
layers  into  dots  on  the  epitaxial  side  of  the  wafer,  followed  by 
alloying  of  these  layers  to  make  ohmic  contacts.  Individual 
mesas  were  defined  by  chlorine  ion-beam-assisted  etching 
using  the  patterned  metal  as  a  mask.^  The  wafer  was  then 
diced  into  chips  ( '-0.018  cm  squares)  and  mounted  in  the 
millimeter-wave  resonator. 

The  resonator  structures  used  to  obtain  the  oscillations 
above  100  GHz  are  represented  schematically  in  Fig.  1(a). 
In  this  structure  the  device  chip  is  mounted  in  the  gap 
between  the  bottom  of  a  whisker  post  and  the  floor  of  a  stan¬ 
dard-height  rectangular  waveguide.  The  gap  is  a  radial 
transmission  line  that  is  terminated  at  its  periphery  by  a  high 
impedance,  and  the  diameter  of  the  post  is  such  that  the  first 
radial-line  resonance  occurs  slightly  above  the  cutoff  fre¬ 
quency  (/c)  of  the  waveguide.  In  the  vicinity  of  this  reso¬ 
nance  the  behavior  of  the  radial  line  is  well  approximated  by 
the  series  combination  of  L  and  shown  in  Fig.  1  (b).  A 
backshort  placed  behind  the  diode  provides  limited  tuning  of 
the  oscillation  power  and  frequency.  The  dc  bias  is  provided 
through  a  coaxial  line  that  also  serves  to  stabilize  the  diode  at 
all  frequencies  below  the  band  of  interest. 

Some  oscillation  results  obtained  in  these  resonators  are 
shown  in  Fig.  2,  along  with  the  results  obtained  at  lower 
frequencies  with  diodes  from  the  same  wafer.  Oscillations 
between  102  and  1 12  GHz  were  obtained  with  a  single  di(xle 
in  a  WR-6  waveguide  ( /^  =  9 1  GHz ) .  In  this  structure,  the 
post  diameter  was  0.063  cm  and  the  height  of  the  gap  was 
about  0.035  cm.  A  50%  scaled-down  version  of  this  resona¬ 
tor  in  WR-3  waveguide  (/^  =  174  GHz)  then  yielded  oscil¬ 
lations  from  191  to  201  GHz.  For  both  of  these  resonators, 
the  frequency  was  varied  primarily  by  changing  the  bias  vol¬ 
tage  within  the  NDC  region.  The  frequency  was  measured 
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FIG.  1.  (a)  Schematic  diagram  of  rectangular 
waveguide  resonator  used  to  generate  100-  and 
200-GHz  oscillations  with  resonant-tunneling  di¬ 
odes.  (b)  Resonator  equivalent  circuit  with  L  and 
Cq  determined  by  the  radial  waveguide  below  the 
whisker  post. 


with  wavemeters,  and  the  power  was  measured  with  calibra¬ 
ted  Schottky  diode  detectors. 

The  frequency  of  all  of  these  oscillations  could  also  be 
varied  somewhat  by  moving  the  backshort,  even  from  posi¬ 
tions  located  many  guide  wavelengths  behind  the  diode.  This 
provides  strong  evidence  that  fundamental  oscillations  have 
been  observed,  as  opposed  to  harmonics  of  an  oscillation  at  a 
lower  frequency.  If  the  oscillations  were  second  or  higher 
harmonics,  then  the  fundamental  would  be  well  below 
and  thus  would  be  unaffected  by  the  position  of  the  back- 
short.  This  argument  rests  on  the  fact  that  the  output  power 
of  the  second  and  all  higher  harmonics  is  considerably 
weaker  (at  least  10  dB  or  so)  than  that  of  the  fundamental. 


FIG.  2.  Room-temperature  oscillation  results  obtained  to  date  with  present 
samples. 


The  oscillation  frequency  is  then  determined  almost  entirely 
by  the  circuit  impedance  at  the  fundamenal  frequency. 

Although  the  power  falls  very  rapidly  with  frequency,  it 
is  clear  from  the  highest  frequency  data  point  in  Fig.  2  that 
they^,,  of  this  device  must  be  above  200  GHz.  This  is  higher 
than  our  previous,  conservative  estimate  for  the  same  device 
(187  GHz).  Therefore,  we  will  now  calculate  more 
carefully  and  consider  ways  to  further  increase  it. 

III.  THEORY 

As  previously  described,/^.,  is  the  highest  frequency  at 
which  a  particular  diode  can  oscillate,  independent  of  the 
circuit.  An  upper  bound  on  can  be  found  using  the 
lumped  equivalent  circuit  shown  in  Fig.  1(b),  consisting  of  a 
(negative)  differential  conductance  G,  a  parasitic  series  re¬ 
sistance  Rs ,  and  a  capacitance  C.  If  we  assume  that  all  of  the 
elements  are  frequency  independent,  the  real  part  of  the  ter¬ 
minal  impedance  of  this  circuit  is  negative  up  to  a  frequency 
given  by 

In  all  double-barrier  diodes  of  interest,  both  Rs  and  C  are 
also  nearly  independent  of  bias  voltage  throughout  the  NDC 
region,  but  G  varies  rapidly  from  a  maximum  negative  value 
near  the  center  of  the  NDC  region  to  zero  at  the  peak 
and  valley  points.  Therefore,  we  consider //i  as  a  function  of 
G  only,  and  find  that  it  gives  a  maximum  value 
/max  =  (47rCR5)”‘  for  G=s  (  —  2/^5)“*.  Ifthe  diode  of  in¬ 
terest  is  such  that  G^,  <  (  —  2/?^ )  ”  *  (as  is  often  the  case  in 
practice),  then/„.,  =/^  (G„., ). 

The  circuit  element  that  varies  the  most  across  a  wafer  is 
Rs-  For  the  most  part  this  reflects  nonuniformities  in  the 
contact  metallization.  A  more  consistent  component  of  is 


1520 


J.  Appl.  Phys.,  Vol.  64.  No.  3. 1  August  1988 


Brown.  Goodhue,  and  Sollner 


1520 


the  resistance  of  the  undepleted  eipilayers  (total  thickness 
)  on  both  sides  of  the  heterostructure.  A  third  important 
component  of/J^  is  the  resistance  caused  by  the  spreading  of 
the  current  from  the  mesa  into  the  much  wider  substrate.  At 
low  frequencies  and  in  a  mesa  that  is  etched  well  into  the 
low-resistivity  substrate  material,  this  spreading  resistance 
can  be  well  approximated  by  p^/(2d)y  where is  the  sub¬ 
strate  resistivity  and  d  is  the  mesa  diameter.®  This  leads  to 
the  following  expression  for  the  total  series  resistance, 

Rs  =  [iPc  ^PeL,)/A  ]  4-  (A/2rf),  (2) 

where is  the  mesa  area,/^^  is  the  specific  contact  resistance, 
and/7,  is  the  epilayer  resistivity.  At  frequencies  high  enough 
that  the  skin  depth  5  in  the  substate  is  much  less  than  the 
effective  diameter  b  of  the  chip,  the  spreading  resistance  is 
increased.  Assuming  also  that  the  skin  depth  is  much  less 
than  the  chip  thickness  h  but  much  greater  than  the  mesa 
diameter,  one  can  write’ 


A  1TO 


•‘Mi)*! 


(3) 


The  geometrical  parameters  in  this  expression  are  shown  in 
Fig.  3.  Equation  (3)  is  frequency  dependent  through  the 
behavior  of  the  skin  depth.  We  will  apply  the  standard 
planar  formula,  S=  [2p/(pco)y^^,  where//  is  the  perme¬ 
ability  and  0)  is  the  angular  frequency.  In  any  case,  the  zero- 
resistance  frequency //j  is  now  the  solution  to  the  equation 


The  solution  will  be  determined  later  in  our  analysis  of  diode 
yinax  • 

The  differential  conductance  of  the  present  device  has 
not  been  measured  accurately  because  the  observed  current- 
voltage  (I-V)  curve  of  this  device  in  the  NDC  region  is  al¬ 
ways  discontinuous,  even  when  it  is  not  oscillating  at  a  fre¬ 
quency  above  The  most  likely  cause  of  this  is  a  spurious 
oscillation  in  the  dc  bias  line,  which  can  occur  if  a  stronger 
oscillation  does  not  exist  in  the  desired  frequency  range.  Be¬ 


cause  a  true  /-  V curve  could  not  be  measured,  we  previously 
estimated  G  by  scaling  the  value  obtained  from  the  contin¬ 
uous  I-V  curve  of  a  much  slower  device.  Shown  in  Fig.  4  is 
the  /-  V  curve  of  such  a  diode  that  was  stabilized  in  a  coaxial 
mount  (the  material  parameters  of  this  diode  are  those  of 
wafer  1,  Ref.  3).  This  is  also  a  diode  for  which  we  know  the 
series  resistance  to  be  negligible  (i.e.,  ).  Mea¬ 

surement  of  the  slope  yields  =  —  2.1  A//AF,  where 
=  Ip  lyy  ^V=  Vp  —  Vy,  and  the  subscripts  P  and  V 
refer  to  the  peak  and  valley  points,  respectively.  Applying 
this  same  prefactor  to  the  present  device  gives  =—15 
mS  for  the  best  A/  /A  V ratio  observed,  and  G,„.,  =  —  1 1  mS 
for  the  worst  ratio  observed. 

The  differential  conductance  can  also  be  obtained  theo¬ 
retically  from  the  single-particle,  stationary-state  model  of 
tunneling.  In  this  model  the  electrical  current  through  the 
structure  is  found  by  summing  the  probability-current  trans¬ 
mission  coefficient  over  the  occupied  states  on  the  cathode 
side  of  the  structure.  Each  state  is  weighted  by  the  longitudi¬ 
nal  velocity  in  that  state,  v  =  li^'dE/dkl,  where  E  is  the 
kinetic  energy  and  k*i'\s  the  longitudinal  component  of  crys¬ 
tal  momentum  on  the  cathode  side.  This  procedure  leads  to 
the  following  expression  for  the  magnitude  of  the  current 
density 

J=0  r  dE^!^T*nE^,y^) 

Jo  kl 

\  1  -h  exp(£J.  —  eV^y  —  eV^  —  Ei^)/kT)  ’ 

where  E^  is  the  longitudinal  component  of  kinetic  energy  of 
electrons  measured  just  to  the  cathode  side  of  the  double¬ 
barrier  structure,  /:  is  the  longitudinal  crystal  momentum 
on  the  anode  side  of  the  structure,  y?  =  em*kT  and  as 

shown  in  Fig.  5,  is  the  potential  difference  across  the 
structure,  Fp  is  the  potential  difference  across  the  depletion 
layer  on  the  anode  side,  and  £5r(£’  J-)  is  the  Fermi  energy  on 
the  cathode  (anode)  side  measured  relative  to  the  conduc¬ 
tion-band  edge  in  the  neutral  region  on  that  side.  The  quanti- 


FIG.  3.  Cross-sectional  view  of  present  diode 
mesa.  The  dimensions  d  i>,  and  h  are  not  drawn 
to  scale  for  clarity. 
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ty  T*T is  the  particle  transmission  coefficient  through  the 
structure,  and  is  assumed  to  be  independent  of  the  transverse 
energy  and  crystal  momentum  in  the  derivation  of  Eq.  (5). 
In  the  case  of  resonant  tunneling,  T*T  has  one  or  more  sharp 
peaks  as  a  function  of  ,  and  the  position  of  these  peaks 
generally  falls  to  lower  Ej^  as  is  increased.  Each  of  these 
peaks  is  conveniently  represented  by  a  (quasibound)  level  in 
the  double-barrier  structure,  as  is  done  for  the  lowest  T^T 
peak  in  Fig.  5.  At  the  cathode  end  of  the  structure,  this  level 
occurs  at  the  value  of  equal  to  the  energy  Ep  of  the  center 
of  the  peak.  The  extension  of  this  level  through  the  entire 
double-barrier  structure  delineates  the  fact  that  the  total  lon¬ 
gitudinal  energy  is  conserved  in  the  resonant  tunneling  pro¬ 
cess. 

When  increasing  bias  voltage  causes  the  center  of  the 
first  r*  r peak  to  fall  near  the  conduction-band  edge  or  near 
the  upper  edge  of  any  other  gap  in  the  density  of  states  on  the 
cathode  side,  the  current  density  given  by  Eq.  (5)  should 
display  a  corresponding  maximum  with  an  NDC  region  on 
its  high-voltage  side.  To  accurately  predict  the  position  of 
this  maximum  and  the  magnitude  of  the  negative  conduc¬ 
tance,  we  must  account  for  the  band  bending  in  the  active 
regions  of  the  device.  The  diagram  in  Fig.  5  shows  qualita¬ 
tively  the  spatial  variation  of  the  conduction-band  edge  at  a 
bias  voltage  near  the  first  current-density  peak.  In  the  deple¬ 
tion  layer  on  the  anode  side  of  the  structure,  the  electron 
potential  energy  is  quadratic,  reflecting  the  fact  that  the 
present  devices  have  a  uniform  doping  concentration  on  the 
anode  side  over  a  distance  significantly  longer  than  the  de¬ 
pletion  length  under  bias  />.  The  potential  in  the  double¬ 
barrier  structure  is  linear  since  this  region  is  nominally  un¬ 
doped  and  is  assumed  to  have  negligible  mobile  space  charge 
density  from  tunneling  electrons.  This  assumption  will  be 
verified  later.  The  potential  profile  in  the  accumulation  layer 
is  not  so  straightforward.  A  numerical  solution  to  Poisson 
equation  indicates  that  the  length  of  this  layer  is  approxi¬ 
mately  30  nm  and  that  is  close  to  100  meV  for  the  device 

shown  in  Fig.  3  at  room  temperature  and  at  bias  voltages  in 
the  NDC  region.  This  solution  was  obtained  by  assuming 


FIG.  4.  /-F  curve  of  scpmic  device  subilized  throughout  the  negative- 
conductance  region. 


that  the  electrons  are  distributed  in  a  three-dimensional  con¬ 
tinuum  throughout  the  layer.  However,  the  solution  to  the 
single  particle  Schrodinger’s  equation  indicates  that  an  accu¬ 
mulation  layer  having  these  dimensions  would  yield  a  maxi¬ 
mum  single-electron  quasibinding  energy  of  only  about  10 
meV.  This  leads  us  to  suppose  that  there  is  a  minimum  ener¬ 
gy  level  below  which  electron  states  do  not  exist,  and  that 
this  level  is  only  slightly  less  in  energy  than  the  neutral  band 
edge  on  the  cathode  side.  For  simplicity  in  the  present  calcu¬ 
lations,  we  will  assume  that  the  minimum  energy  occurs  ex- 
actly  at  the  neutral  band  edge  shown  in  Fig.  5.  If  we  make  the 
additional,  reasonable  assumptions  that  e(  Fj*,  +  F£) )  >£?• 
and  e(  F^/  -f  F^, )  >/:r,  the  current  density  integral  can  be 
written 

J=)5  r  dE^{kl/kl)TmE^,V^) 

JeV, 

Xln  [1  +exp(£^^-heF^  ^E^)/kTY  (6) 

For  the  particle  transmission  probability,  we  adopt  the 
Breit-Wigner  resonant  form, 

T*T=  (kl/k  1  )n  [ - £,)"  +  rg ] 
where  F  is  the  half-width  at  the  half-maximum  value. 
Strictly  speaking,  this  expression  is  only  an  upper  limit  on 
the  transmission  probability  as  required  by  unitarity.'*  Nev¬ 
ertheless,  we  apply  it  in  lieu  of  a  more  accurate  numerical 
method  because  it  leads  to  a  closed-form  expression  for  the 
current  density,  and  because  it  can  analytically  account  for 
the  effects  of  band  bending  and  scattering.  It  we  use  the  mod¬ 
el  shown  in  Fig.  5,  the  band  bending  causes  Ep  to  decrease 
with  Vfp  almost  exactly  as  eF^//2,  and  therefore  causes  it  to 
decrease  with  Vp  asEp  =  £„  —  aV  {^,  where  Eq  is  the  peak 
of  the  r*rcurve  (i.e.,  the  quasibound  state  energy)  at  zero 
bias,  a  =  eW(eNp/2€)*^^,  fV  is  the  width  of  the  double¬ 
barrier  structure,  and  e  is  the  dielectric  constant  in  the  anode 
region.  In  deriving  this  expression  the  dielectric  constant  in 
the  barriers  was  set  equal  to  e.  The  error  caused  by  this  as- 


FIG.  5.  Electron  potential  energy  at  the  conduction-band  edge  of  the  pres¬ 
ent  device  near  the  bias  voltage  required  for  the  first  current  peak.  The  po¬ 
tential-energy  drop  eV^  across  the  accumulation  layer  is  exaggerated  for 
clarity. 
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sumption  will  be  negligible  in  the  present  devices  because  of 
the  thinness  of  the  barriers  in  these  devices. 

Scattering  can  be  accounted  for  by  generalizing  the  re¬ 
cent  result  of  Stone  and  Lee/"* 


111 

n 


Io__EL_. 

r 7*  —  Ep)^  T\ 


(7) 


where  =  To  -f  r5,  and  F^  is  related  to  the  scattering 
time  through  the  relation  F^  =  fi/r^ .  For  the  present 
double-barrier  structures  at  room  temperature,  so 

that  we  can  treat  the  logarithmic  term  as  a  constant,  evaluate 
it  at  the  peak  of  F  ♦  F,  and  approximate  the  current  integral  in 
Eq.  (6)  as 


J  =  pTo{n/l  +  tan- '  [  (E^  -  eV^  )/Tr\) 


X\n{\  +  eV^  —  Ep)/kT\).  (8) 

The  differential  conductance  is  related  to  the  current  density 
\)y  G  =  AdJ  /dVj,  where 

yj  =  Vo  +  V^+V^  +  (£>  -  E'’p)/e. 

Upon  differentiating  Eq.  (8)  we  find 


G=pr^ 


dEp 


dVj 


X 


•rrln{l  +  exp[(£'^f +  -£:;.)Ar]} 

r\.+  (Ep-eV^)^ 


7r/2  +  t^n-'[(Ep-eV^)/rr]  1 

kT{exp[(Ep  -  E'ir  -  eVJ/kT]  +  1}J  ‘ 

As  has  been  noted  previously,  the  current  density  is  largely 
independent  of  the  scattering  width  F5  provided  that 
FsKEf-  However,  in  the  NDC  region,  the  expression  for  G 
depends  strongly  on  this  parameter.  This  is  to  be  expected, 
since  in  this  model  the  NDC  occurs  as  the  bias  voltage  drives 
Ep  below  eV^ .  Thus,  the  /-F  curve  in  the  NDC  region  re¬ 
flects  the  shape  of  F*F and  how  it  decreases  with  Vj. 

The  calculation  of  the  capacitance  in  the  lumped-circuit 
model  is  straightforward.  If  we  assume  that  no  charge  exists 
in  the  double-barrier  structure,  that  the  depletion  approxi¬ 
mation  is  valid  on  the  anode  side,  and  that  the  dielectric 
constant  €  is  uniform  throughout  the  active  region,  the  solu¬ 
tion  to  Poisson’s  equation  yields 

D+W=[W^+ae/eNi,)(yj-y^)]''^-  (10) 

The  capacitance  of  these  regions  is  given  by  =  €A  / 
(Z)  -f  HO,  and  the  total  device  capacitance  is  C  = 

(Cpiy  -f  ),  where  is  attributable  to  the  accumulation 
layer.  Because  should  exceed  by  about  a  factor  of 
10,  its  contribution  will  henceforth  be  ignored  and  we  will 
assume  the  total  capacitance  to  be  equal  toCpiy.  Confidence 
in  the  depletion  approximation  is  provided  by  the  fact  that 
the  peak  current  density  Jp  through  the  present  devices  is 
typically  less  than  4  X  10^  A  cm“^,  so  that  the  mobile  space- 
charge  density  under  bias  {Jp/eu^  =  1.2X10'^  cm"^  for 
Uj  =  2  X  10’  cm  s“ ' )  is  far  less  than  the  background  doping 
concentration,  Np  =  2X  10*’  cm“^.  The  assumption  of  an 
uncharged  quantum  well  is  definitely  good  near  the  current 
valley,  but  as  the  bias  point  approaches  the  current  peak,  the 
well  should  be  charged  maximally  by  tunneling  electrons. 


The  magnitude  of  this  charge  will  be  evaluated  later  for  one 
particular  diode. 

At  oscillation  frequencies  for  which  the  period  is  com¬ 
parable  to  the  time  of  electron  transit  across  the  total  deple¬ 
tion  length,  the  lumped  circuit  analysis  becomes  inaccurate 
and  should  be  replaced  by  a  distributed  theory.  Since  we  are 
only  interested  here  in  maximum  oscillation  frequencies 
(near  which  the  oscillation  power  is  necessarily  weak),  we 
adopt  the  small-signal  transit  time  theory  first  applied  to 
BARITT  diodes.'^  The  real  part  of  the  terminal  impedance 
is  given  by 

^  oCp  -f 

+  (11) 

where  cr  is  the  injection  conductance  {dJ  /dF)  of  the  double¬ 
barrier  structure,  Fis  the  electric  field  across  the  structure,  6 
is  the  transit  angle  and  equals  coD  /v^ ,  is  the  average  drift 
velocity  across  the  depletion  layer  on  the  anode  side  of  the 
barriers,  Rsif)  is  given  by  Eq.  (3),  and  Cp  =  eA  /D.  The 
frequency  at  which  R  p  vanishes  defines  the  transit-time-lim¬ 
ited,  zero-resistance  frequency.  Because  the  double-barrier 
structure  is  assumed  to  be  uncharged,  we  can  write 

<j={D+]V)-^=(D+W)^, 

dV j  A 

where  G  is  given  by  Eq.  (9). 


IV.  DEVICE  PARAMETER  EVALUATION 

The  double-barrier  diode  that  yielded  the  best  perfor¬ 
mance  has  the  material  composition  shown  in  Fig.  3.  From 
this  information  and  the  measured  /-  V  curve  of  the  device, 
all  of  the  desired  device  parameters  can  be  evaluated.  We 
first  calculate  the  total  series  resistance  of  Eq.  ( 3 )  for  a  A-fim 
mesa.  Separate  measurements  indicate  that  for 
Np  =  2x  10”  cm“^  in  the  region  below  the  ohmic  contact, 
can  be  no  less  than  10“^  ft  cm’.  Using  this  lower  limit  we 
find  a  contact  resistance  of  8  ft.  The  epilayer  resistance  Repi 
is  found  to  be  about  4  ft  for  =  500  nm  and  a  0.009 

ft  cm  on  each  side  of  the  barriers.  The  high-frequency 
spreading  resistance  component  in  Eq.  ( 3 )  can  be  simplified 
to  6x  10”^  ft  if  we  assume  p^  =  0.0015  ft  cm  (for  the 

=  2X  10'**  cm”’  substrate),  b  =  0.02  cm,  and  h  ==  0.01 
cm.  The  sum  of  these  components  is  Rs  =  14  ft  at  200  GHz 
and  16  ft  at  600  GHz,  respectively.  However,  nonuniformi¬ 
ties  in  the  contact  quality  often  resulted  in  much  larger  val¬ 
ues  of p^  than  the  minimal  value  assumed  here.  Higher  series 
resistance  causes  the  current  peak  to  shift  to  higher  voltage. 
Devices  having  peak  voltages  shifted  by  about  0.25  V  or 
more  (R^  ^  50  ft)  did  not  oscillate  in  the  WR-3  resonator. 

Shown  in  Fig.  6  is  the  7-  V  curve  for  the  present  device 
obtained  by  applying  Eq.  ( 8 )  and  the  expression  for  the  total 
voltage  across  the  device  F  7-  =  -^JARs.  We  assume 
here  that  R^  equals  the  minimum  value  14  ft.  In  applying 
this  expression,  Fq  was  obtained  by  demanding  that  the  peak 
current  be  the  same  as  that  observed  experimentally.  The 
value  of  F()  =  1.65  meV  obtained  by  that  procedure  is  much 
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greater  than  the  value  of  =  0.5  meV  predicted  by  a  sepa¬ 
rate  (numerical)  calculation  oiT*T  for  the  present  struc¬ 
ture  in  which  the  effective  mass  in  the  AlAs  barriers  was 
assumed  to  be  equal  to  the  bulk- AlAs  value,  =  0.15  mo- 
It  is  much  closer  to  the  numerical  result  of  Tq  =  2.1  meV 
obtained  with  an  effective  mass  in  the  barriers  equal  to  the 
bulk-GaAs  value,  m*  =  0.067  rriQ.  The  latter  calculation 
also  determined  Eq  =  140  meV.  The  value  of  Ts  =  2.5  meV 
was  calculated  from  data  on  the  mobility  of  quantized  inver¬ 
sion  layers  in  MODFETs,’^-'**  assuming  Ts—H/t, 
r  =  and  =  7000  cm^  s“'.  The  value  for 

==  —  18  meV  corresponds  to  fully  ionized  but  neutral 
2  X  10'^  cm”^  material  at  295  K.  Finally,  =  0.065  V  was 
obtained  by  demanding  that  Eq.  (8)  yield  the  same  peak 
voltage  as  the  experimental  curve.  Since  we  are  only  interest¬ 
ed  in  the  rather  narrow  range  of  the  NDC  region,  was 
kept  fixed  in  this  calculation. 

The  determination  of  the  two  quantities,  Fq  and  , 
allows  us  to  calculate  the  capacitance,  the  differential  con¬ 
ductance,  and  the  charge  distribution  in  the  double-barrier 
structure.  Substituting  =  65  mV,  W=^19  nm,  and 
Vj  =  0.85  V  in  Eq.  ( 10),  we  find  Z>  -f-  75.5  nm,  so  that 

Ca  19  fF  for  a  4-//m-diam  mesa.  This  parameter  was  found 
to  be  21  fF  in  our  previous  estimation  when  the  quantum 
well  was  assumed  to  be  charged.^  The  value  Fq  =  1.65  meV 
and  the  assumed  value  of  F^  allow  us  to  obtain  the  curve  G 
vs  Vj  shown  in  Fig.  7.  Note  that  the  =  —  26  mS  is 
significantly  larger  than  the  value  —  13  mS  obtained  by  ap¬ 
plying  the  rule  =  —  2.1  A/ /A  F  to  the  experimental /- 
V curve  in  Fig.  4.  However,  it  is  much  smaller  in  magnitude 

than  the -  1000  mS  predicted  by  a  previous  theory.'^ 

Band  bending  on  the  anode  side  of  the  barriers  reduces  the 
latter  prediction  by  an  order  of  magnitude,  and  scattering 
reduces  it  further  by  about  a  factor  of  3. 

To  obtain  the  charge  distribution  in  the  present  double¬ 
barrier  structure,  we  employ  a  variation  of  the  model  that 
was  used  to  derive  Eq.  ( 5 ) .  Here  we  sum  the  electron  station¬ 


ary-state  probability  density  Pe^(^)  = 

where  is  the  wave  function  at  a  longitudinal  energy  Ei^ , 

over  the  occupied  states  on  the  cathode  side  of  the  structure 

without  weighting  these  states  by  their  longitudinal  velocity. 

If  we  adopt  the  model  of  band  bending  shown  in  Fig.  5,  this 

procedure  leads  to  the  following  expression  for  the  charge 

density: 


n(x)  =  r  dEi^ 

€  JeV^ 


Pe,M 

El^ 


Xln[l  -fexp(£J<-h^F^  —  £^)A7'].  (12) 

The  probability  density  is  obtained  by  first  modeling  the  real 
potential  energy  profile  as  a  series  of  0. 1-nm-wide,  uniform- 
potential  layers.  Starting  at  the  accumulation  layer  bound¬ 
ary,  we  evaluate (x)  in  each  successive  layer  by  applying 
the  transfer-matrix  method."  Finally,  Eq.  (12)  is  numeri¬ 
cally  integrated  using  the  value  for  found  above  and  the 
value  of  F,  at  the  current  peak.  The  resulting  charge  distri¬ 
bution  is  shown  in  Fig.  8  along  with  the  assumed  (linear) 
potential  energy  profile  in  the  double-barrier  structure.  No¬ 
tice  that  the  charge  peaks  near  the  center  of  the  quantum 
well  as  expected,  but  the  maximum  magnitude  is  only  about 
1  X  10"  cm”^.  When  this  distribution  is  integrated  over  the 
double-barrier  structure,  a  sheet  charge  density  of  3.9  X  10'*^ 
cm"^  is  obtained.  This  is  so  much  less  than  the  total  sheet 
density,  s  1.5x10"  cm-%  in  the  depletion  layer  that  its 
effect  on  the  resonant-tunneling  process  is  negligible.  A  de¬ 
tailed  calculation  verifies  this  assertion  by  showing  the 
charge  to  depress  the  assumed  potential  profile  by  less  than  1 
meV  across  the  entire  structure. 


V.  ANALYSIS  OF  MAXIMUM  OSCILLATION 
FREQUENCY 

We  can  now  determine  for  the  present  diode  struc¬ 
ture  and  for  a  feasible  variation  of  it.  In  relation  to  the  theory 


VOLTAGE  (V) 


FIG.  6.  Experimental  (solid)  and  theoretical  (dashed)  /-^curves  for  pres¬ 
ent  4-^m-diam  device  at  room  temperature.  The  arrows  denote  the  direc¬ 
tion  of  switching  that  occurs  in  the  experimental  hysteresis  loop. 


FIG.  7.  Theoretical  /-Kand  G-Kcurves  of  the  present  4-/im-diam  device  at 
room  temperature. 
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FIG.  8.  Electron  potential  energy  and  charge  distri¬ 
bution  across  the  present  double-barrier  structure 
with  0.9-V  bias  applied  at  room  temperature. 


developed  in  Sec.  Ill,  is  the  peak  value  of  the  zero- 
resistance  frequency  over  the  possible  range  of  differen¬ 
tial  conductance  values,  G  <  .  It  is  instructive  to  calcu- 

late/nax  for  each  of  the  three  theoretical  models  considered: 
the  lumped-circuit  model  [Eq.  (1)],  the  lumped-circuit 
model  modified  by  the  frequency-dependent  skin  effect  [Eq. 
(4)  ],  and  the  combined  transit-time,  skin-effect  model  [Eq. 
(11)].  We  also  suspect  that/^.^^  will  depend  strongly  on  the 
diode  geometry,  so  that  three  different  mesa  diameters  will 
be  anayzed. 

Shown  in  Figs.  9(a),  9(b),  and  9(c)  are  curves  of for 
mesa  diameters  of  4, 2,  and  1  //m,  respectively.  The  range  of 
possible  operating  conditions  for  each  diode  size  lies  left  of 
the  dashed  line  G  =  G^.^ .  The  resulting  values  of are 
denoted  (/„„)s,  and  (/„„)r  for  the  lumped- 

circuit,  skin-effect,  and  transit-time  models,  respectively, 
and  are  summarized  in  Table  I.  Note  that  curves  of  ( )  j- 
are  given  for  two  different  values  of  drift  velocity, 
Vj  =2X  10^  cm  s”*  and  4x  10^  cm  s"*.  Although  we  do 
not  precisely  know  the  average  drift  velocity  across  the  75- 
nm  depletion  layer,  the  two  assumed  values  are  thought  to  be 
representative  of  the  possible  range.  They  are  both  greater 
than  the  saturation  velocity  in  GaAs  (-- 1 X 10^  cm  s"  * )  be¬ 
cause  the  double-barrier  structure  injects  the  elecrons  with  a 
relatively  large  kinetic  energy  {^200  meV),  which  results 
in  a  very  large  initial  group  velocity. 

Several  aspects  of  these  theoretical  results  deserve  atten¬ 
tion.  First,  it  is  gratifying  that  the  range  of  ( )  7.,  214  to 
244  GHz,  for  the  4-/im-diam  mesa  lies  above  our  highest 
experimental  result.  However,  notice  that  the  transit-time 
delay  lowers  /^.x  significantly  compared  with  the  lumped- 
circuit  and  skin-effect  models.  In  addition,  varies  no¬ 
ticeably  with  the  mesa  diameter.  For  example,  when  the  di¬ 
ameter  is  reduced  from  4  to  2  //m,  all  values  increase 
somewhat  because  the  spreading  resistance  is  a  much 
smaller  fraction  ofT?^  in  the  2-/im-diam  device.  This  reflects 
the  fact  that  the  spreading  resistance  varies  with  area  much 


slower  than  inversely.  However,  there  is  very  little  difference 
i^yinax  between  the  2-  and  l-//m-diam  devices  because  the 
relative  magnitude  of  the  spreading  resistance  is  becoming 
negligible  at  these  diameters.  Finally,  it  is  clear  that  a  change 
in  the  material  parameters  of  this  device  that  increased  G^.,^ 
but  kept  Rs  and  C  fixed  would  do  little  or  nothing  to  en- 
hance/„,.,i .  This  is  because  the  Gn,,,^  calculated  for  our  pres¬ 
ent  device  is  very  close  to  the  value  of  G  which  produces  the 
highest  A  more  promising  material  change  would  be 
one  that  decreases  Rs- 

A  simple  variation  of  the  present  device  structure  that 
should  reduce  Rs  significantly  is  shown  diagrammatically  in 
Fig.  10.  The  doping  concentration  on  the  cathode  side  is 
rapidly  increased  to  the  largest  practical  value  (Six  10'** 
cm~^)  after  a  few  accumulation-layer  lengths  on  the  cath¬ 
ode  side,  and  is  increased  to  the  same  level  after  about  100 
nm  on  the  anode  side.  This  will  reduce  the  contact  specific 
resistance  to  about  1  X  10"^^  ft  cm^  and  should  cut  R^^^  in 
half,  thereby  lowering  /?5  to  about  5  ft  at  200  GHz  and  7  ft 
at  600  GHz.  We  assume  that  the  other  parameters  are  identi¬ 
cal  to  those  of  the  present  device. 

Theoretical  curves  of  the  zero-resistance  frequency  for 
the  proposed  device  are  given  in  Fig.  11.  First  notice  the 
large  increase  in  all/^.^^  values  compared  with  present  de¬ 
vices  of  the  same  diameter.  For  example,  the  4-//m-diam 
device  now  gives  (/^.^  )l  =  593  GHz  and  a  transit-time- 
limited  range  (/^.x  )7-  =  369-426  GHz.  Because  the  rela¬ 
tive  contribution  of  the  spreading  resistance  is  now  much 
larger,  there  is  a  great  increase  in  all values  as  the  mesa 
diameter  is  decreased.  For  example,  the  upper  limit  of  the 
( /max )  T  range  increases  to  568  and  641  GHz  in  the  2-  and  1- 
fim  mesas,  respectively.  Also  notice  that  the  transit -time- 
limited  curves  all  peak  at  a  value  of  G  greater  than  in 
contrast  to  the  present  device.  This  means  that  for  these  low 
values  of  Rsy  significantly  faster  devices  can  be  obtained  by 
increasing 

In  principle  could  be  increased  in  one  of  three 
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FIG.  9.  Zero>resistance  frequency  curves  for  ( a)  ( b)  2-^ni-,  and  ( c) 

l-/ini-diam  mesas  of  the  present  sample.  The  three  quantities 

and  (/Mii)r  maximum  oscillation  frequencies  for  the 

lumped-circuit  model,  the  lumped-circuit  model  with  skin  effect,  and  the 
transit-time  model  with  skin  effect,  respectively.  The  quantity  is  the 
maximum  differential  conductance. 


TABLE  I.  Values  of  for  the  three  models  under  consideration. 


Mesa  diameter 
(//m) 

(GHz) 

t/m*»  Is 
(GHz) 

C/in*»  )r 

(GHz) 

iv^  =  2-4 X  10^ cm  s“‘) 

4 

293 

Present  device: 
274 

214-244 

2 

310 

310 

237-275 

1 

317 

317 

244-283 

4 

593 

Proposed  device: 
467 

369-426 

2 

676 

634 

475-568 

1 

758 

730 

535-641 

ways:  ( 1 )  by  material  changes  that  increase  the  resonant- 
tunneling  current  density  by  some  fraction  while  enhancing 
the  leakage  current  density  by  a  lesser  fraction;  (2)  material 
changes  that  reduce  only  the  leakage  current  density;  or  ( 3 ) 
operational  changes  (such  as  a  reduction  in  temperature) 
that  reduce  the  leakage  current  density.  We  will  consider 
only  the  first  method  because  it  is  the  simplest  to  compute  in 
the  present  theoretical  formalism,  and  because  it  has  the 
practical  appeal  of  increasing  the  oscillation  power  density. 
The  most  plausible  material  change  is  to  increase  the  doping 
concentration,  and  hence  the  Fermi  level,  on  the  cathode 
side  of  the  heterostructure.  For  example,  suppose  that 
Np  =  lx  10***  cm“^  throughout  the  cathode  side  of  the  bar¬ 
riers.  In  this  case  £  J.  =  -f  39  meV,  and  if  we  assume  that  the 

remains  the  same  as  for  the  lower  doping,  Eqs.  (8)  and 
(9)  predict  the /-FandG-K  curves  of  Fig.  12.  Notice  that  7;, 
has  increased  to  four  times  its  value  in  Fig.  7,  and 
Gmmx  =  —  125  mS.  All  of  the  peaks  in  the //^  curves  of  Fig. 
11  could  then  be  achieved.  Thus,  the  upp>er  limit  of  the 
( /max )  T  range  would  be  453, 642,  and  750  GHz  for  the  4-,  2- 
,  and  l-//m  devices,  respectively. 

VI.  DISCUSSION 

This  analysis  is  not  intended  to  provide  a  complete  opti¬ 
mization  for  the  double-barrier  diode.  Instead,  we  have 
pointed  out  some  simple  changes  in  the  diode  material  struc¬ 
ture  that  should  increase  the  maximum  oscillation  frequency 
by  at  least  a  factor  of  2.  A  complete  optimization  would 
require  consideration  of  other  material  parameters,  such  as 
the  barrier  and  quantum-well  thickness,  and  the  doping  con¬ 
centration  on  the  anode  side  of  the  structure.  Very  little  is 
known  about  the  effect  of  these  parameters  on  the  diode  be¬ 
havior,  especially  their  influence  on  the  excess  current  [i.e., 
the  current  beyond  that  predicted  by  Eq.  (8)).  As  seen  in 
Fig.  6,  the  excess  current  is  the  major  component  of  the  total 
current  at  and  above  the  valley.  It  is  possibly  associated  with 
scattering  from  impurities  in  the  barriers,  a  mechanism  for 
which  transverse  crystal  momentum  is  not  conserved  in  the 
tunneling  process.  The  component  should  also  be  very  fjist, 
so  that  it  would  affect  the  device  behavior  by  reducing  the 
peak-to-valley  ratio  and  by  reducing  . 
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FIG.  10.  Doping  concentration  profile  for 
present  (dotted)  and  proposed  (solid) 
doublC'barrier  device. 


FIG.  11.  Zeixyresisunce  frequency  curves  for  (a)  (b)  2-/iin%  and 

(c)  l-/im-diani  mesas  of  the  proposed  sample. 
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FIG.  12.  Theoretical  /-Fand  G-Fcurves  of  pro¬ 
posed  sample  at  room  temperature.  This  sample 
would  have  a  uniform  doping  concentration  of 
No  ==  2x  10'’cm“^and  1 X  10’"  cm“^  on  the  an¬ 
ode  and  cathode  sides  of  the  double-barrier  struc¬ 
ture. 


Another  factor  that  we  have  ignored  is  the  finite  lifetime 
of  electrons  in  the  quasibound  states  of  the  quantum  well. 
According  to  the  uncertainty  principle,  each  electron  wave- 
packet  builds  up  and  decays  in  the  well  in  a  time  r>^ro, 
where  Fq  is  the  half  width  of  the  T*T  curve.  It  has  often  been 
assumed  that  this  corresponds  to  a  zero-resistance  frequen¬ 
cy, =  (2^rT)"',  which  would  be  near  400  GHz  for  the 
present  device.  While  the  finite  lifetime  will  definitely  slow 
down  the  tunneling  process,  its  effect  of  is  not  clear. 
Intuitively  we  expect  that  the  lifetime  in  the  well  should 
cause  the  ac  conduction  current  to  lag  behind  the  ac  voltage. 
This  is  clearly  an  inductive  effect.  A  double-barrier  diode 
equivalent  circuit  containing  an  intrinsic  inductance  has 
been  proposed  by  Gering  et 

All  of  the  oscillations  measured  to  date  above  30  GHz 
have  had  output  power  less  than  1  mW,  and  it  is  unlikely  that 
this  level  will  ever  be  greatly  exceeded  in  a  single  diode.  To 
do  so  would  require  a  much  larger  device  than  those  consid¬ 
ered  here.  Such  a  device  would  be  practically  impossible  to 
stabilize  at  low  frequencies  because  its  would  be  at  least 

—  100  mS.  A  more  promising  way  to  increase  the  power 
would  be  to  use  the  double-barrier  structure  as  an  injector 
for  a  transit-time  oscillator.^®  In  this  case  both  the  material 
parameters  of  the  structure  and  the  doping  level  on  the  an¬ 
ode  side  are  tailored  to  achieve  a  negative  terminal  resistance 
according  to  Eq.  (11).  This  can  be  achieved  with  the  device 
biased  in  the  p)Ositive  resistance  region  if  the  depletion  length 
yields  a  transit  angle  6  between  tt  and  27r  in  the  frequency 
range  of  interest.  A  significant  increase  in  the  device  area  will 
then  enhance  the  output  power  without  causing  low-fre¬ 
quency  instabilities.  Another  approach  to  increasing  the 
power  is  to  combine  the  output  of  several  diodes  that  operate 
in  the  NDC  region  but  are  independently  biased  and  stabi¬ 
lized.  For  example,  one  could  design  a  monolithic  circuit  in 
which  each  of  the  diodes  had  its  own  resonator,  bias 
network,  and  a  planar  antenna  to  couple  out  the  radiation. 


The  analysis  of  such  a  circuit  has  been  given  recently  by 
Mink.^‘ 

VII.  SUMMARY 

We  have  demonstrated  the  first  fundamental  oscilla¬ 
tions  in  a  resonant  tunneling  diode  above  1(X)  GHz,  with  our 
highest  frequency  occurring  at  201  GHz.  This  has  led  to  a 
new  estimate  of  293  GHz  for  the  upper  limit  of in  the 
present  4-/im-diam  device.  It  should  be  compared  with  the 
value  =  187  GHz  obtained  previously  from  empirical 
estimates  of  The  negative  conductance  is  now  ob¬ 
tained  with  a  phenomenological  theory  containing  two  fit¬ 
ting  parameters.  This  upper  limit  is  reduced  by  skin  effects 
and  transit-time  delay  to  244  GHz,  if  we  assume  the  drift 
velocity  equals  4  X 10^  cm  s”  *  across  the  depletion  layer.  We 
propose  that  a  simple  change  in  the  doping  profile  of  the 
present  device  which  reduces  the  series  resistance  should  in¬ 
crease  the  transit-time-limited/^,;^  to  about  430,  570,  and 
640  GHz  for  a  4-,  2-,  and  l-/im-diam  device,  respectively.  If 
the  current  density  could  also  be  increased  without  signifi¬ 
cantly  diminishing  the  pcak-to- valley  ratio,  these values 
would  increase  further  by  about  10%. 
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APPENDIX  D 


Harmonic  multiplication  using  resonant  tunneling 

T.  C.  L  G.  Sollner,  E.  R.  Brown,  W.  D.  Goodhue,  and  C.  A.  Correa 

Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  Massachusetts  02 173 
(Received  4  March  1988;  accepted  for  publication  1  July  1988) 

We  propose  and  demonstrate  an  odd-harmonic  resistive  multiplier  that  exploits  the  unique 
current-voltage  (/-F)  relation  of  resonant-tunneling  diodes.  It  is  shown  that  efficient  odd- 
harmonic  conversion  is  possible  and  that  even  harmonics  do  not  appear  because  of  the 
antisymmetry  of  the  /-K  curve.  In  addition,  with  the  proper  choice  of  resonant-tunneling 
structure  and  pump  amplitude,  most  of  the  harmonic  output  power  can  be  confined  to  a  single, 
odd-harmonic  frequency.  Fifth-harmonic  multiplication  has  been  demonstrated  with  output  at 
21.75  GHz  and  a  power  conversion  efficiency  of  0.5%  in  a  coaxial  circuit.  Circuit  simulation 
using  an  improved  /-Fcurve  yields  a  fifth-harmonic  efficiency  of  2.7%,  suggesting  that  useful 
harmonic  output  in  the  millimeter  and  submillimeter  spectrum  should  be  achievable. 


Harmonic  multipliers  are  often  used  for  the  generation 
of  electromagnetic  waves  at  frequencies  above  those  con¬ 
veniently  available  from  fundamental  oscillators.  For  exam¬ 
ple,  radio  astronomy  receivers  at  frequencies  above  50  GHz 
often  use  power  from  a  multiplier  as  the  local  oscillator,  *  and 
multipliers  are  a  primary  source  of  power  for  molecular 
spectroscopy  in  the  submillimeter  spectrum.^ 

Harmonic  multipliers  have  traditionally  used  either  the 
voltage-dependent  resistance  (a  varistor)  or  capacitance  (a 
varactor)  associated  with  either  metal-semiconductor 
(Schottky  diodes)  or  p-n  junctions.  Capacitive  nonlineari¬ 
ties  have  been  favored  because  of  the  higher  available  con¬ 
version  efficiency  between  the  power  delivered  by  the  pump 
and  that  available  at  the  desired  harmonic.  It  was  shown  by 
Page^  that  a  resistive  nonlinearity  containing  no  negative 
differential  resistances  can  produce  a  conversion  efficiency 
of  at  most  \/n^j  where  n  is  the  harmonic  number.  On  the 
other  hand,  varactors  are  limited  by  the  Manley-Rowe  rela¬ 
tions,^  which  in  this  case  allow  a  maximum  conversion  effi¬ 
ciency  of  unity. 

In  this  communication,  we  propose  and  demonstrate 
the  use  of  resonant-tunneling  diodes  as  varistors  for  harmon¬ 
ic  generation.  (Some  of  the  basic  ideas  and  preliminary  re¬ 
sults  were  first  published  by  Sollner,  Brown,  and  Goodhue^ 
in  1987.)  One  advantage  of  these  devices  over  Schottky  di¬ 
odes  is  the  presence  of  negative  differential  resistance  regions 
in  the  current-voltage  (/-K)  curve,  so  it  is  possible,  as  dis¬ 
cussed  by  Andreyev,^  to  surpass  the  limits  found  by  Page  for 
resistive  mixers.  Another  advantage  is  that,  under  the  proper 
conditions,  it  is  possible  to  confine  most  of  the  harmonic 
output  power  to  a  single  harmonic  frequency.  For  Schottky- 
diode  multipliers,  larger  voltages  are  always  generated  at 
lower  harmonics,  which  must  be  properly  reactively  termin¬ 
ated.  Avoiding  these  additional  “idler”  circuits  will  greatly 
simplify  circuit  design  and  should  therefore  increase  effi¬ 
ciency,  especially  in  the  millimeter  and  submillimeter  wave¬ 
length  range  where  it  is  difficult  to  avoid  significant  resistive 
loss.  Also,  no  dc  bias  is  required,  further  simplifying  the 
circuit.  For  resonant-tunneling  diodes  available  today,  out¬ 
put  frequencies  of  several  hundred  gigahertz  are  expected  to 
be  feasible,  and  perhaps  even  higher  frequencies  will  be  avail¬ 
able  with  optimized  diodes.  This  part  of  the  frequency  spec¬ 


trum  now  suffers  from  a  lack  of  available  solid-state  sources, 
so  resonant-tunneling  multipliers  could  fill  this  need,  at  least 
for  low-power  sources. 

Resonant  tunneling  occurs  when  two  or  more  semi¬ 
transparent  barriers  are  placed  closely  enough  together  that 
charge  carriers  interact  coherently  with  them.  The  resonant 
states  formed  between  the  barriers  greatly  increase  the  tun¬ 
neling  probability  of  carriers  with  energy  equal  to  that  of  the 
states,  resulting  in  peaks  and  valleys  in  the  current  as  the 
voltage  across  the  structure  is  increased.  Typically,  the  bar¬ 
riers  are  made  from  an  epitaxially  grown  material  of  larger 
band  gap,  such  as  AlAs  grown  on  GaAs.  A  typical  /-Feurve 
for  a  double-barrier  resonant-tunneling  diode  is  shown  in 
Fig.  1.  More  details  on  the  fabrication  of  resonant-tunneling 
diodes  have  been  published  by  Goodhue  et  al? 


VOLTAGE  (V) 


FIG.  1.  Current-voltage  (/-F)  curve  of  the  resonant-tunneling  diode  that 
was  used  as  a  resistive  harmonic  multiplier.  Oscillations  were  occurring  in 
the  negative  differential  resistance  regions,  causing  the  departures  from  a 
smooth  curve.  These  oscillations,  at  about  10  GHz,  apparently  did  not  begin 
when  the  device  was  pumped  at  frequencies  above  1  GHz. 
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The  presence  of  a  peak  and  valley  in  the  /-  Kcurve,  com¬ 
bined  with  the  overall  antisymmetry  of  the  /-  V  curve  about 
the  origin  [i.e.,  I(V)  =  —  ^  offer  the  potential  for 

efficient  odd-harmonic  generation  with  an  unbiased  reso¬ 
nant-tunneling  diode.  If  the  diode  is  pumped  with  an  ac  sig¬ 
nal,  so  that  the  peak  amplitude  of  the  voltage  across  the 
diode  occurs  above  the  resonant  current  peak,  at  least  three 
local  maxima  occur  in  the  diode  current  waveform  over  one 
cycle,  corresponding  to  third  or  higher  odd-harmonic  gener¬ 
ation. 

It  is  interesting  to  note  that  the  velocity-field  relation  of 
n-type  GaAs  is  similar  to  the  current-voltage  relation  shown 
in  Fig.  1.  This  velocity-field  relation  arises  from  the  scatter¬ 
ing  of  hot  electrons  between  the  T  and  L  valleys,  and  a  num¬ 
ber  of  workers  have  used  this  effect  for  frequency  multiplica¬ 
tion.*’  The  main  disadvantage  is  that  electrons  require  time 
of  the  order  of  a  picosecond  to  return  to  the  T  valley,  limiting 
the  useful  frequencies  to  below  100  GHz.  The  frequency  lim¬ 
its  for  resonant  tunneling  will  be  discussed  below. 

A  resonant-tunneling  diode  with  the  /-  Fcurve  shown  in 
Fig.  1  was  mounted  in  a  coaxial,  whisker-contacted  package 
and  placed  in  the  circuit  shown  in  Fig.  2.  All  components 
were  connected  together  with  50-ft  coaxial  line,  and  the  cir¬ 
cuit  was  relatively  free  of  resonances  to  22  GHz.  The  oscilla¬ 
tor  power  at  4.25  GHz  was  adjusted  to  give  optimum  conver¬ 
sion  efficiency  at  the  fifth  harmonic,  resulting  in  the 
spectrum  shown  in  Fig.  3.  The  power  available  at  the  funda¬ 
mental  (4.25  GHz)  was  measured  by  replacing  the  reso¬ 
nant-tunneling  diode  with  a  short.  The  best  fifth-harmonic 
conversion  efficiency  was  0.5%.  As  expected,  the  power 
converted  to  even  harmonics  is  more  than  40  dB  below  the 
odd-harmonic  output. 

To  estimate  the  conversion  efficiency  expected  from  a 
resonant-tunneling  diode,  the  /-  V  curve  of  Fig.  1  was  fitted 
to  a  piecewise  continuous  curve,  and  the  circuit  of  Fig.  2  was 
simulated  using  lumped  circuit  elements.  We  have  assumed 
that  the  only  power  source  occurs  at  the  fundamental  fre¬ 
quency  of  the  oscillator;  however,  the  solution  to  the  full 
nonlinear  problem  requires  sources  that  account  self-consis- 
tently  for  the  harmonics  generated  by  the  nonlinearities  of 
the  diode.  The  neglect  of  these  terms  greatly  simplifies  the 


BROADBAND  RESONANT 

MICROWAVE  DIRECTIONAL  TUNNEUNG 

OSCILLATOR  COUPLER  DIODE 


FIG.  2.  Schematic  circuit  used  to  measure  the  harmonic  spectrum  shown  in 
Fig.  3. 


FIG.  3.  Power  spectrum  generated  by  a  resonant-tunneling  diode  in  the 
circuit  of  Fig.  2.  The  amplitude  of  the  pump  at  4.25  GHz  has  been  adjusted 
to  optimize  the  fifth-harmonic  output  power  at  21.25  GHz.  The  absence  of 
even  harmonics  is  an  indication  of  the  high  degree  of  antisymmetry  of  the 
/-  V  curve. 

calculation,  but  causes  some  quantitative  discrepancy.  The 
calculated  fifth-harmonic  conversion  efficiency  is  0.3% 
rather  than  the  observed  value  of  0.5%,  and  the  calculated 
third-harmonic  output  power  is  below  that  observed  by  a 
similar  amount.  It  is  possible  that  the  differences  may  arise 
from  our  neglect  of  the  harmonic  sources  mentioned  above, 
or  perhaps  from  ignoring  the  voltage  dependence  of  the  ca¬ 
pacitance  that  occurs  in  the  resonant-tunneling  diode.  Still, 
the  agreement  with  the  measured  spectrum  of  Fig.  3  is  suffi¬ 
cient  to  give  us  confidence  that  the  numerical  model  is  a 
good  approximation  to  the  experimental  measurement. 

The  model  described  above  can  be  used  to  predict  the 
performance  of  a  multiplier  when  the  /-F  curve  of  the  reso¬ 
nant-tunneling  device  has  been  arranged  to  increase  har¬ 
monic  generation.  The  characteristics  that  are  desirable  are 
a  large  peak-to-valley  ratio  and  a  shape  that  provides  more 
equal  spacing  in  time  between  the  maxima  in  the  current 
waveform.  The  simulated  /-Feurve  shown  in  Fig.  4  has  not 
been  optimized  for  these  characteristics,  but  it  is  a  consider¬ 
able  improvement  over  the  curve  of  Fig.  1.  Figure  5  shows 
the  power  spectrum  at  the  terminals  of  the  diode  when  the 
voltage  amplitude  and  source  resistance  have  been  opti¬ 
mized  for  fifth-harmonic  conversion.  This  simulation  yields 
d  fifth-harmonic  conversion  efficiency  of  2.7%,  and  addi¬ 
tional  optimization  is  still  possible.  The  calculated  efficiency 
is  near  the  limit  of  \/n^  derived  by  Page^  for  positive-con¬ 
ductance  nonlinearities,  and  perhaps  further  improvements 
could  surpass  it.  For  example,  the  filters  could  be  arranged 
so  that  the  fundamental  and  harmonics  could  be  terminated 
with  different  impedances. 

The  efficiency  of  multipliers  using  /-  V  curves  similar  to 
that  of  Fig.  1  drops  quickly  for  harmonics  above  the  fifth 
because  there  is  only  one  peak  and  valley  in  each  quadrant  of 
the  curve;  however,  it  is  possible  to  have  several  peaks  in  a 
resonant-tunneling  /-  Feurve.  Two  peaks  would  make  possi¬ 
ble  the  efficient  generation  of  up  to  the  ninth  harmonic,  three 
peaks  allow  the  thirteenth  harmonic,  and  so  on.  (Recently, 
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FIG.  4.  Hypothetical  /-Fcurvc  that  gives  better  harmonic  generation  effi¬ 
ciency  than  the  experimental  curve  of  Fig.  1. 
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FIG.  5.  Power  spectrum  calculated  using  the  /-  V curve  of  Fig.  4.  The  con¬ 
version  efficiency  to  the  fifth  harmonic  is  nearly  an  order  of  magnitude  larg¬ 
er  than  that  calculated  using  the  experimental  /-Feurve  of  Fig.  1. 


Sen  et  al."  have  suggested  that  two  resonant-tunneling  di¬ 
odes  in  series  can  produce  a  large  fifth-harmonic  output  in  a 
somewhat  different  application. )  It  should  also  be  feasible  to 
adjust  the  pump  amplitude  to  produce  most  of  the  harmonic 
power  in  a  single  harmonic  if  the  current  in  the  higher-vol¬ 
tage  valleys  can  be  kept  sufficiently  small.  Thus,  a  single 
harmonic  multiplier  could  convert  power  below  100  GHz  to 
the  terahertz  range. 

The  upper  frequency  limit  for  a  given  resonant-tunnel¬ 
ing  diode  oscillator  has  been  studied  in  some  detail,*^  but  the 
multiplier  has  not  been  carefully  analyzed  to  date.  An  esti¬ 
mate  can  be  reached  through  reasoning  leading  to  the  maxi¬ 
mum  oscillation  frequency  .  This  is  the  frequency  above 
which  there  is  no  negative  resistance  observable  from  the 
device  terminals,  so  the  peaks  in  the  current  waveform,  so 
essential  to  the  concepts  presented  here,  would  not  be  pres¬ 
ent.  This  frequency  is  given  by 

A.,  =  ( l/27rC)  L,  ] 

Here  C  is  the  device  capacitance  parallel  with  the  maximum 
(negative)  conductance  and  is  the  resistance  in 

series  with  this  parallel  combination.  Currently  available 
resonant-tunneling  diodes  have  values  approaching  300 
GHz.*^  One  of  these  diodes  has  recently  been  used  as  a 
triplet  by  Batelaan  and  Frerking*^  to  produce  250^  W  at  191 
GHz  with  0.61%  efficiency. 

In  summary,  a  promising  new  approach  for  resistive 
harmonic  multipliers  has  been  demonstrated  with  resonant¬ 
tunneling  diodes.  Several  advantages  for  high-harmonic 
multipliers  are  available,  and  it  is  expected  that  these  ideas 


can  be  extended  to  obtain  usable  power  in  the  terahertz  re¬ 
gion. 
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APPENDIX  E 


Observation  of  millimeter-wave  oscillations  from  resonant  tunneling  diodes 
and  some  theoretical  considerations  of  ultimate  frequency  limits 

T.  C.  L.  G.  Sollner,  E.  R.  Brown,  W.  D.  Goodhue,  and  H.  Q.  Le 

Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  Massachusetts  02173 

(Received  30  October  1986;  accepted  for  publication  8  December  1986) 

Recent  observations  of  oscillation  frequencies  up  to  56  GHz  in  resonant  tunneling  structures 
are  discussed  in  relation  to  calculations  by  several  authors  of  the  ultimate  frequency  limits  of 
these  devices.  We  find  that  calculations  relying  on  the  Wentzel-Kramers-Brillouin  (WKB) 
approximation  give  limits  well  below  the  observed  oscillation  frequencies.  Two  other 
techniques  for  calculating  the  upper  frequency  limit  were  found  to  give  more  reasonable 
results.  In  one  method  we  use  the  solution  of  the  time-dependent  Schrodinger  equation 
obtained  by  Kundrotas  and  Dargys  [Phys.  Status  Solidi  B  134,  267  ( 1986)  ],  while  in  the  other 
we  use  the  energy  width  of  the  transmission  function  for  electrons  through  the  double-barrier 
structure.  This  last  technique  is  believed  to  be  the  most  accurate  since  it  is  based  on  general 
results  for  the  lifetime  of  any  resonant  state.  It  gives  frequency  limits  on  the  order  of  1  THz  for 
two  recently  fabricated  structures.  It  appears  that  the  primary  limitation  of  the  oscillation 
frequency  for  double-barrier  resonant  tunneling  diodes  will  be  imposed  by  intrinsic  device 
circuit  parameters  and  by  the  transit  time  of  the  depletion  layer  rather  than  by  time  delays 
encountered  in  the  double-barrier  region. 


There  has  been  considerable  debate  on  the  subject  of  the 
speed  of  devices  based  on  resonant  tunneling.  By  comparing 
high-frequency  current  response  measurements  with  the  ob- 
.served  dc  characteristics,  we  have  previously  established' 
that  the  response  time  is  as  short  as  100  fs.  Luryi"  subse¬ 
quently  pointed  out  that  the  oscillation  frequency  of  a  dou¬ 
ble-barrier  diode  (DBD)  employing  resonant  tunneling 
should  be  limited  by  the  well  charging  and  discharging  time. 
He  calculated  this  time  to  be  40  ps  for  the  structure  mea¬ 
sured  in  Ref.  1,  corresponding  to  a  frequency  limit  of  ap¬ 
proximately  /  =  l/(27rr)  =  4  GHz.  We  have  recently  ob¬ 
served  oscillations  with  r  as  short  as  2.8  ps  (corresponding 
to  56  GHz).  This  frequency  is  more  than  an  order  of  magni¬ 
tude  greater  than  predicted  by  Luryi’s  calculation.  We  will 
present  evidence  that  this  discrepancy  lies  in  the  Wentzel- 
Kramers-Brillouin  (WKB)  approximation  used  by  Luryi. 

The  complete  problem  of  electron  kinetics  through  a 
DBD  is  complex,  involving  self-consistent  solution  to  the 
Poisson  equation  and  the  time-dependent  Schrodinger  equa¬ 
tion,  including  elastic  and  inelastic  scattering.  Although  sev¬ 
eral  groups  are  currently  attempting  solutions,  no  results  are 
yet  available.  Since  guidelines  are  needed  to  help  direct  ex¬ 
perimental  programs,  we  will  here  attempt  to  set  a  bound  on 
the  maximum  frequency  of  resp)onse  by  calculating  parts  of 
the  problem  separately  and  finding  a  characteristic  time  for 
each  part. 

The  structure  of  a  DBD  and  the  resulting  current-vol¬ 
tage  {I-V)  curve  are  shown  in  Fig.  1.  Oscillations  are  ob¬ 
tained  by  biasing  the  device  into  the  negative  differential 
resistance  region  of  the  /-  V  curve  while  it  is  embedded  in  a 
suitable  high-frequency  circuit.  Figure  2  shows  a  sketch  of 
the  waveguide  circuit  used  and  its  equivalent  circuit  with  the 
DBD.  The  details  of  this  circuit  and  the  material  parameters 
of  the  DBD  have  been  reported  elsewhere.^ 

Table  I  summarizes  the  material  and  electrical  param¬ 
eters  for  three  waveguide  oscillators.  The  capacitance  is  cal¬ 
culated  including  the  size  of  the  depletion  region  when  the 


device  is  biased  in  the  negative  differential  resistance  region. 
The  series  resistance  is  also  calculated,  but  it  has  been  mea¬ 
sured  in  one  sample  and  found  to  agree  with  the  calculation. 
The  actual  oscillation  frequency/^^  is  determined  primarily 
by  the  waveguide  circuit.  The  maximum  oscillation  frequen¬ 
cy  consistent  with  the  equivalent  circuit  shown  in  Fig.  2  oc¬ 
curs  when  the  conductance  in  parallel  with  the  DBD  from 
the  external  circuit  is  made  arbitrarily  small.  This  maximum 
frequency  determined  by  the  intrinsic  circuit  elements  of  the 
DBD  is 

=  (l/27rC)(  -  -  GL.  . 

where  is  the  maximum  negative  conductance. 

The  transit  time  across  the  depletion  region  is  closely 
related  to  the  capacitance  in  the  equivalent  circuit.  A  wider 
depletion  region  reduces  the  capacitance  of  the  device,  but 
also  introduces  a  phase  shift  in  the  current  relative  to  the 
applied  voltage  just  as  would  storage  in  the  well.  While  this 
phase  shift  could  be  used  to  achieve  a  negative  conductance 
as  in  a  barrier  injection  and  transit  time  diodc,^  we  will  as¬ 
sume  here  that  the  phase  shift  must  be  kept  small  to  utilize 
the  intrinsic  negative  conductance  of  the  DBD.  This  time 
delay  yields  a  frequency  limit  that  we  arbitrarily  define 
as  /depi  =  (I/2^depi  )•  This  limit  is  calculated  in  Table  I 
assuming  the  saturated  carrier  velocity  of  10^  cm/s.  Since 
the  scattering  mean  free  path  is  comparable  to  the  depletion 
region  thickness,  the  actual  velocity  of  carriers  in  the  deple¬ 
tion  region  will  probably  be  higher  than  this  value,  so  this 
estimate  of /depi  is  conservative. 

Another  effect  that  could  limit  the  maximum  oscillation 
frequency  is  charge  storage  delay  in  the  quantum  wcU.  In 
other  words,  if  each  carrier  must  spend  a  certain  period  of 
time  in  a  quantum  state  in  the  well,  then  this  state  lifetime 
introduces  a  phase  difference  between  the  current  and  vol¬ 
tage  and  would  limit  the  intrinsic  response  of  the  DBD.  At 
first  glance  this  state  lifetime  would  appear  to  be  rather  long. 
Bohm^  has  solved  this  problem  in  the  WKB  approximation 
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FIG.  1.  Energy-band  diagrams  showing  the  origins  of  negative  resistance  in 
adoDble-barrierdiode.  Thebarriersarecausedby  the  higher  band  gap  of  the 
AlGaAs  alloy.  These  barriers  create  resonant  states  in  the  well  at  energies 
Et ,  and  the  current  has  local  maxima  when  these  resonances  are  at  the  same 
energy  as  most  of  the  elections  incident  from  the  emitter  on  the  left.  The 
lifetime  of  an  electron  in  one  of  these  resonant  states  and  the  transit  time 
through  the  depletion  region  are  possible  limits  on  the  speed  of  the  negative 
resistance. 


and  found  the  time  to  be  given  by  the  classical  period  of  a 
carrier  in  the  well  divided  by  the  transmission  probability  of 
the  barrier,  i.e., 

A/=(^^+l/16^^)r 

a«^  V, 

where  r  is  the  classical  period, 

t=i4={  <2m(E-V)dx, 
oE 

and  the  barrier  transmission  probability  is  6  where 

«/b«fTier  H 

Here  £  is  the  energy  of  the  carrier,  F  b  the  potential  of  the 
well  or  barrier,  and  m  is  the  carrier  effective  mass.  In  Table  I, 
the  frequencies  corresponding  to  A/  are  given  as/wKB  =  1/ 
(IvEt)  for  the  several  DBD’s  fabricated.  These  frequencies 
are  well  below  the  observed  oscillation  frequencies,  so  it  is 
clear  that  there  must  be  something  wrong  with  this  simple 
approach.  Of  course,  the  WKB  approximation  does  not 
strictly  apply  since  the  potential  in  the  DBD  changes  rapidly 
on  the  sc^e  of  the  carrier  wavelength.  But  probably  the  more 
severe  limitation  stems  from  the  fact  that  the  carrier  is  con- 


FIG.  2.  Waveguide  structure  used  to  make  oscillators  from  double-liarrier 
diodes  (DBD's)  and  the  associated  equivalent  circuit.  The  circuit  intrinsic 
to  the  DBD  also  limits  the  maximum  oscillation  frequency. 

tained  between  two  barriers  separated  by  less  than  the  car¬ 
rier  wavelength.  Thus  it  is  no  longer  reasonable  to  think  of  a 
carrier  as  being  reflected  between  two  barriers  since  it  is  lo¬ 
cated  partially  outside  both  barriers  at  any  time.  This  greatly 
increases  the  escape  probability,  as  we  shall  see. 

So  far  we  have  considered  only  the  motion  of  electrons 
in  the  well.  However,  Luryi^  calculates  the  current  through  a 
single  barrier  in  the  WKB  approximation  to  And  a  resistance 
associated  with  the  barrier,  and  from  the  capacitance  across 
the  barrier  a  time  constant  results, 

/2C=  (6a“U  (cgs  units) , 

where  A  =  h  /V2m^,  (f>  is  the  average  barrier  height,  m  is  the 
electron  efiective  mass,  €  is  the  barrier  dielectric  constant, 
a”*  =  ykr/e^a  137,  d  is  the  barrier  thickness,  and  c  is  the 
speed  of  light.  Luryi  asserts  that  this  time  constant  repre¬ 

sents  the  lifetime  of  carriers  in  the  well  and  hence  the  limit  to 
oscillation  frequencies.  This  limit  is  labeled  f  ^  in  Table  I. 
This  limit  is  also  considerably  less  than  the  observed  oscilla¬ 
tion  frequencies. 

The  lifetime  of  carriers  in  the  well  can  also  be  estimated 
from  a  recent  exact  (numerical)  solution  to  the  time-depen- 
dent  Schrodinger  equation  for  a  thin  quantum  well  by  Kun- 
drotas  and  Dargys.^  The  authors  solved  for  the  probability 
of  escape  from  a  delta-function  quantum  well  after  the  in¬ 
stantaneous  q)plication  of  an  electric  fleld.  Their  barriers 
were  infinitely  thick,  so  the  escape  time  will  be  much  larger 
than  that  from  our  thin  barriers;  hence  this  calculation  pro¬ 
vides  an  upper  bound  on  the  DBD  escape  time.  They  found 
the  time  required  for  the  escape  probability  to  reach  near- 
unity  values  was  significantly  shorter  than  that  predicted  by 
standard  approximations.  Their  results  yield  a  cut-oflT  fre¬ 
quency  labeled /^d  in  Table  I.  Note  that /kd  is  much  higher 
than  /wKB  though  the  WKB  calculation  was  for  bar¬ 

riers  of  finite  thickness. 

Perhaps  the  most  accurate  method  for  estimating  the 
lifetime  of  an  electron  in  the  well  is  to  calculate  the  energy 
width  of  the  transmission  resonance  through  the  double-bar¬ 
rier  structure.  Then  the  lifetime  of  this  resonance  is  given  by 
an  uncertainty  relation,  r  =  as  recently  pointed  out 
by  Ricco  and  Azbef  and  by  Coon  and  Liu.®  This  very  gen¬ 
eral  result,  originally  derived  for  atomic  and  nuclear  transi¬ 
tions,  applies  to  any  resonant  state.^  We  have  calculated  the 
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TABLE  I.  Measured  and  calculated  parameters  for  three  different  wafers  of  double-barrier  diodes. 


Sample 

1 

2 

3 

Material  parameters 

Barrier  material 

AlAs 

Ga^  7  Alfl  3  As 

AlAs 

Barrier  thickness  (nm) 

2.5 

3.0 

1.5 

Well  thickness  (nm) 

4.5 

4.5 

4.5 

Doping  outside  barriers  (cm"') 

lx  10'" 

2X10'^ 

2X10'" 

Electrical  parameters 

Peak-to- valley  ratio,  300  K 

1.7/1 

1.3/1 

3.5/1 

Peak  current  density  (  X  lO*  A  cm"^) 

0.8 

1.2 

4.0 

Depletion  layer  at  bias  (nm) 

15 

30 

70 

Capacitance  (fF)* 

100 

50 

20 

Max.  negative  conductance  (mS)" 

5.0 

8.0 

13.0 

Series  resistance  (H)* 

10 

15 

15 

Oscillation  characteristics 

dcbias/f,F^  (mA,V) 

0.7,0.40 

2.7,0.32 

3.0,0.95 

/^.(GHz)‘ 

20.7 

43.7 

56.4 

THEORETICAL 

Max.  oscillation  frequency 

/c.k(GH2)' 

35 

70 

200 

/*^(GHz)‘' 

1000 

500 

230 

y*WKB  (GHz)* 

0.01 

220 

6.0 

A  (GHz)' 

0.01 

60 

6.0 

Ad  (GHz)* 

130 

300 

100 

At  (GHz)* 

80 

2100 

800 

"Typical  values  for  a  circular  mesa  of  4/im  diameter. 

**  Maximum  observed  fundamental  oscillation  frequency. 

Xrc  =  (2;rC)-'(  -GLJ’''. 

**  From  depletion  layer  drift  time. 

*  WKB  estimate  of  charge  storage  time  in  well. 

'^From  Ref.  2  for  charge  storage  time  in  well. 

•  From  Ref.  6,  for  electron  lifetime  in  well  with  infinitely  thick  barriers. 

**  From  calculation  of  energy  width  of  transmission  through  double-barrier  structure. 


transmission  as  a  function  of  incident  electron  energy  for  the 
structures  of  Table  I  using  essentially  the  method  outlined 
by  Tsu  and  Esaki.*®  A  frequency  corresponding  to  the  life¬ 
time  is  listed  in  Table  I  as  This  frequency  is  generally 
higher  than / kd  because  of  the  finite  barriers.  ( It  seems  to  be 
slightly  lower  for  sample  1.) 

An  important  conclusion  from  this  work  is  that  the 
WKB  method  appears  to  consistently  overestimate  the  es¬ 
cape  time  of  an  electron  in  a  narrow  quantum  well  with  an 
electric  field  applied.  The  most  satisfactory  estimate  of  the 
well  storage  time  comes  from  considering  the  width  of  the 
transmission  resonance  (leading  to  )  since  this  deriva¬ 
tion  is  firmly  grounded  in  the  theory  of  scattering  from  reso¬ 
nant  states.  Transit  time  effects  such  as  /^epi  must  also  be 
considered,  as  must  those  limitations  arising  purely  from  the 
device  equivalent  circuit,  At  this  point  in  the  develop¬ 
ment  of  DBD  oscillators,  these  intrinsic  circuit  elements 
present  the  first  hurdle  to  higher  frequencies  and  must  be 
optimized  with  the  depletion  layer  drift  time  in  mind.  In 
view  of  the  results  discussed  here,  there  appear  to  be  no  fun¬ 


damental  or  practical  impediments  to  DBD  oscillators  up  to 
several  hundred  gigahertz. 
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Persistent  photoconductivity  in  quantum  well  resonators 

T.  C.  L.  G.  Sollner,  H.  Q.  Le,  C.  A.  Correa,  and  W.  D.  Goodhue 

Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  Massachusetts  02173 

(Received  8  February  1985;  accepted  for  publication  8  April  1985) 

We  have  made  the  first  observation  of  persistent  photoconductivity  in  resonant  tunneling 
structures.  The  spectral  dependence  suggests  that  it  arises  from  Si  DX  centers  in  the  AlGaAs 
barriers,  and  calculations  based  on  a  simple  model  using  this  assumption  agree  well  with  the 
observations.  This  effect  has  been  useful  in  investigating  the  charge  distribution  and  electric  fields 
near  the  heterojunction  interface  and  in  determining  the  barrier  parameters.  The  model  should 
help  in  the  design  and  fabrication  of  optimized  resonant  tunneling  devices. 


Resonant  tunneling  structures  (quantum  well  resona¬ 
tors)  have  recently  been  shown  to  possess  some  intriguing 
characteristics,  including  high-speed  charge  transport 
(r  <  10~  s)  and  broad  regions  of  the  negative  resistance  that 
can  be  used  to  produce  oscillations  at  microwave,  and  per¬ 
haps  higher,  frequencies.’*^  These  heterostructure  devices 
contain  GaAs-AlGaAs  interfaces  that  are  capable  of  macro¬ 
scopic  segregation  of  carriers  because  of  the  difference  in 
band  gaps.  Carriers  in  the  AlGaAs  diffuse  to  the  GaAs  with 
lower  band  gap,  leaving  behind  an  ionized  impurity,  and 
establish  a  dipole  layer  across  the  interface.  In  addition  to 
hydrogenic  impurities,  it  is  well  known^  that  in  AlGaAs, 
electron  traps  called  DX  centers  exist  which,  once  optically 
ionized,  remain  so  for  long  times  at  low  temperatures,  giving 
rise  to  a  persistent  photoconductivity.  This  effect  has  been 
used  in  quantum  wells  to  study  charge  transport  along  the 
interface,"’  and  we  will  show  here  that  it  can  be  helpful  in 
clarifying  resonant  tunneling  as  well. 

In  this  letter  we  report  the  first  observation  of  persistent 
photoconductivity  in  resonant  tunneling  structures.  We 
have  observed  a  large  effect  upon  the  current-voltage  (/-F) 
curve,  and  the  spectral  dependence  of  the  effect  confirms 
that  it  is  due  to  Si  DX  centers  in  the  AlGaAs  barrier  regions. 
A  simple  model  agrees  well  with  the  observations  and  pro¬ 
vides  further  insight  into  the  physics  of  heterojunctions.  It  is 
possible  for  the  first  time  from  /-  Fcurves  alone  to  determine 
the  conduction-band  discontinuity  and  the  barrier  thickness 
uniquely. 

A  resonant  tunneling  structure  consists  of  thin  barrier 
regions  of  Al^^  Ga,  _  As  separated  by  thin  regions  of  GaAs. 
As  Fig.  1  shows,  the  barriers  to  electron  flow  occur  because 
of  the  higher  energy  of  the  conduction  band  of  the 
Al^Ga,_^As.  These  barriers  act  as  partially  transparent 
mirrors  to  electrons,  the  charge  transport  taking  place  by 
tunneling  through  the  thin  barriers.  Two  of  these  barriers, 
placed  sufficiently  closely  together  that  electron  wave  func¬ 
tions  are  coherent  across  them,  form  a  sort  of  electron 
Fabry-Perot  resonator  with  peaks  in  the  electron  transmis¬ 
sion  (current)  as  a  function  of  incident  electron  energy  (vol¬ 
tage).  Current  peaks  occur  when  the  voltage  across  the  dou¬ 
ble-barrier  region  is  roughly  equal  to  lEf/e,  where  E,  is  one 
of  the  allowed  energies,  measured  from  the  Fermi  level,  of 
electron  motion  normal  to  the  barriers  for  carriers  confined 
to  the  well.  The  lifetime  of  a  carrier  in  the  well  is  shortened 
by  the  possibility  of  tunneling  out,  so  this  level  is  broadened, 
as  shown  in  Fig.  I .  In  our  structures,  only  one  level  occurs, 
which  is  marked  E^  in  the  figure. 


DX  centers  in  AlGaAs  are  believed  to  be  electron  traps 
that  are  neutral  until  the  electron  is  ionized.  With  the  elec¬ 
tron  removed,  a  configurational  change  occurs  so  that  re¬ 
combination  is  unlikely  without  the  assistance  of  a  phonon 
of  large  momentum;  thus,  at  low  temperatures,  relaxation 
times  are  extremely  long.  (Some  repopulation  may  occur 
from  long  wavelength  photons,  as  discussed  by  Nathan  et 
alf  and  by  Kastalski  and  Hwang;^  however,  we  have  not 
observed  repopulation  with  either  1.06-  or  1.9-//m  pulses.) 
The  effect  of  the  presence  of  positively  charged  centers  in  the 
barriers  on  the  electron  energy  is  shown  by  the  solid  lines  in 
Fig.  1,  whereas  the  dashed  line  shows  the  energy  in  the  ab¬ 
sence  of  ionized  centers.  If  the  conduction  band  in  the  well  is 
below  the  Fermi  level,  a  significant  number  of  electrons  ac¬ 
cumulate  in  the  well  at  zero  bias.  But  at  bias  levels  above  a 
few  tens  of  millivolts,  they  quickly  tunnel  out,  so  we  have 
ignored  this  accumulation  in  our  calculations  and  in  Fig.  1 . 
The  net  effect  is  that  the  well  is  lowered  relative  to  the  Fermi 
level  and  the  effective  barriers  are  lowered.  The  conse¬ 
quences  of  this  potential  shift  will  be  discussed  in  more  detail 
below. 

The  resonant  tunneling  structures  were  grown  by  mo¬ 
lecular  beam  epitaxy  on  an  wafer  of  GaAs.  The  barriers 
and  the  well  widths  are  about  40  A.  The  regions  outside  the 
barriers  were  doped  with  Si  to  1  X  10”’  cm”*',  whereas  the 
barriers  and  well  were  nominally  undoped  (total  impurity 
content  <  10’'  cm”').  The  A1  mole  fraction  in  the  barriers 
was  20  ±  5%.  Although  no  Si  was  grown  in  the  barriers, 
appreciable  diffusion  of  Si  (but  not  Al)  does  take  place  during 
overgrowth  of  the  top  5000-A  layer  of  GaAs.  Using  the  dif¬ 
fusion  constants  of  Tatsuta  ei  al.,^  we  estimate  an  average  Si 
concentration  in  the  barriers  of  2-5  X  10’^  cm”'.  We  note, 
however,  that  the  Si  diffusion  in  GaAs  is  about  half  as  fast  as 
in  the  Alo.3Gao.7As  measured  by  Tatsuta  et  al.,  so  there 


FIG.  1.  Electron  energy  as  a  function  of  position  in  a  quantum  well  struc¬ 
ture  with  positive  charge  in  the  barriers  and  the  electrons  outside  the  barrier 
well  region.  The  dashed  lines  show  the  conduction  band  in  the  absence  of 
ionized  centers. 
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FIG.  2.  (a)  Observed  and  (b)  calculated  /♦  V curves  for  the  structure  shown  in 
Fig.  1  for  several  exposures  to0.85-/im  light.  The  long-dashed  curve  is  unex¬ 
posed,  the  solid  curve  has  had  an  arbitrary  exposure,  Q,  the  short-dashed 
curve  has  exposure  10  Q,  and  the  broken  top  curve  has  exposure  100  Q. 

should  be  significantly  less  Si  in  the  well.  The  heterojunction 
interface  may  also  inhibit  diffusion.  Ohmic  contacts  are  pro¬ 
vided  on  the  top  and  bottom  of  the  wafer,  and  isolation  is 
accomplished  by  the  wet  etching  of  3-/im-diam  mesas. 
Whisker  contact  is  made  to  the  small  mesas. 

Illumination  was  provided  by  a  quartz-halogen  lamp 
and  a  grating  monochromator  with  the  resolution  set  to  25 
nm.  A  low-pass  filter  was  used  to  eliminate  higher  order 
dispersion,  and  an  electronic  shutter  controlled  the  expo¬ 
sure.  The  sample  was  mounted  in  a  mechanical  refrigerator 
with  the  mesa  at  45*  to  the  optical  axis.  The  sample  was 
cooled  to  20  K  in  the  dark.  I-V curves  were  recorded  after 
every  exposure  to  produce  a  sequence  of  curves,  some  of 
which  are  shown  in  Fig.  2(a).  Aside  from  the  curve  taken 
before  exposure,  each  curve  shown  has  an  exposure  ten  times 
greater  than  the  curve  below  it.  The  top  curve  has  had  nearly 
all  XhtDX  centers  ionized.  After  the  persistent  photoconduc- 
tive  effect  had  saturated,  the  sample  was  warmed  to  room 
temperature  and  recooled  for  another  experiment  at  a  differ¬ 
ent  wavelength. 

The  effects  observed  in  Fig.  2(a)  can  be  explained  by  the 
model  shown  in  Fig.  1  and  a  calculation  of  the  /-  V  curves  at 
finite  temperature  in  a  way  similar  to  that  suggested  by  Tsu 
and  Esaki,^  except  that  accumulation  and  depletion  regions 
are  included.  Figure  2(b)  shows  /-T  curves  calculated  using 
this  model  for  the  same  relative  exposures  as  the  experimen¬ 
tal  curves  of  Fig.  2(a)  and  the  DX  center  density  as  deter¬ 
mined  below.  The  parameters  used  in  the  model  were  consis¬ 
tent  with  the  growth  constants  of  the  structure,  as  described 
above,  except  that  the  temperature  used  in  the  model  was  60 
K  rather  than  the  physical  temperature  of  the  device  of  20  K. 
The  fact  that  little  change  in  the  observed  /-  V  curves  is  seen 
below  about  100  K  supports  this  evidence  that  some  electron 


heating  is  occurring,  although  the  heating  mechanism  is  un¬ 
known.  The  agreement  with  the  observed  curves  is  good, 
considering  the  simplicity  of  the  model.  The  higher  predict¬ 
ed  current  density  is  believed  to  arise  from  the  omission  of 
scattering  in  the  theory.  The  current  increase  at  higher  vol¬ 
tages  is  not  completely  understood,  but  if  it  is  due  to  impuri¬ 
ty-assisted  tunneling,  then  the  observed  increase  with  in¬ 
creasing  DX  center  density  would  seem  reasonable.  More 
details  of  this  model  will  be  published  in  the  future. 

It  remains  to  relate  the  observed  shift  of  peak  current  at 
resonant  bias  voltage  to  the  incident  photon  number.  I-V 
curve  calculations  similar  to  those  of  Fig.  2(b)  show  that  the 
shift  in  current  and  voltage  of  the  resonant  peak  is  a  linear 
function  of  the  shift  in  the  well  potential,  for  A^'  <  £, 
(see  Fig.  1).  The  voltage  shift  can  be  understood  from  the 
resonant  energy  of  the  structure,  £„  which  follows  linearly 
the  level  of  the  well  bottom.  The  peak  current  is  also  calcu¬ 
lated  to  vary  linearly  with  A(t>,  although  it  is  not  so  obvious 
why  this  should  occur. 

It  is  easy  to  show  from  Poisson’s  equation  that  the  shift 
in  the  well  potential  because  of  positive  charges  in  the  bar¬ 
riers  is  linearly  related  to  the  volume  density,  n,  of  ionized 
centers  in  the  barriers  by 

=  (1) 

l€ 

where  d  is  the  barrier  thickness,  e  is  the  electron  charge,  €  is 
the  dielectric  constant  of  the  barrier  region,  and  /tf  is  the 
Thomas-Fermi  screening  length  in  the  GaAs  surrounding 
regions.  Thus,  the  change  in  resonant  current,  7^^,  and  vol¬ 
tage,  A  Vpp ,  are  expected  to  vary  linearly  with  the  ionized  DX 
center  density  n.  (We  measure  and  as  the  total  cur¬ 
rent  and  voltage,  respectively,  between  the  positive  and  neg¬ 
ative  resonance  peaks.)  We  will  assume  for  simplicity  that 
the  charge  density  in  the  two  barriers  is  the  same,  although 
the  slight  asymmetry  about  zero  in  the  /-  V  curve  suggests 
some  difference.  This  effect  is  minimized  by  measuring  the 
sum  of  the  changes  in  the  /-K  curve  for  both  positive  and 
negative  directions. 


PHOTON  NUMBER  (arbitrary  units) 

FIG.  3.  DifTerence  in  current  between  positive  and  negative  resonant  peaks 
in  the  l-V curve  as  a  function  of  exposure  (in  arbitrary  units)  to  0.85-/im 
photons.  The  solid  curve  fit  to  the  data  points  is  given  by  Eq.  (3). 
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FIG.  4.  Points  show  the  observed  photon  energy  dependence  of  the  relative 
ionization  cross  section  for  the  persistent  photoconductive  effect  described 
in  this  letter.  A  typical  error  bar  is  shown  on  one  data  point.  The  solid  line  is 
from  Ref.  9  for  Si  centers  in  AI04  Gao  6  As.  The  agreement  confirms  that 
the  effect  observed  here  is  caused  by  Si  DX  centers  in  the  barrier  regions  of 

Fig.  I. 

The  density  of  ionized  centers  can  be  related  to  the 
number  of  incident  photons,  Q,  by  the  usual  exponential  ap¬ 
proach  to  saturation, 

^  ^  A^oO/(Av)|  1  -  exp  [  -  (C7y  +  a, ]Q /A]] 

oy(Av)  +  <7,  (/tv) 

where  is  the  total  volume  density  of  DX  centers,  A  is  the 
illuminated  area,  and  af(hy)  and  a,(/rv)  are  the  photon-in¬ 
duced  ionization  and  trapping  cross  sections,  respectively. 
The  latter  can  be  expected  to  be  negligible,  except  perhaps  at 
longer  wavelengths  where  some  quenching  may  occur.^’^ 
From  Eq.  (2)  and  the  linearity  discussed  above,  we  arrive  at 
an  expression  for  the  change  in  peak  current. 

AI^  =^/“' [  1  -  exp(  -a^Q/A)].  (3) 

A  fit  to  the  observed  as  a  function  of  Q  yields  ^4/ and 

the  relative  ionization  cross  section  0/(Av).  Figure  3  shows 
the  fit  of  Eq.  (3)  to  the  observations  at  a  wavelength  of  0.85 
//m.  A  similar  expression  holds  for  A  and  the  fit  to  ex¬ 
perimental  data  is  good,  if  not  quite  so  close  as  that  shown  in 
Fig.  3.  It  is  clear  from  these  data  that  the  model  describes  the 
observed  effects  very  well. 

The  sequence  of  /-F  curves  for  different  exposures,  as 
shown  in  Fig.  2,  permits  a  determination  of  the  material 
parameters  that  is  not  possible  with  a  single  I-V  curve. 
Changing  parameters  in  the  /-  V  curve  calculations  shows 
that  Vpp  is  determined  primarily  by  the  well  width,  but  Ipp  is 
sensitive  to  both  barrier  height  and  barrier  width,  so  that  the 
latter  two  quantities  cannot  be  determined  uniquely  from 
the  dark  curve  alone.  With  the  added  variable  of  charge  in 
the  barrier,  however,  the  additional  information  from  AIpp 
and  A  Vpp  resolves  this  ambiguity,  allowing  barrier  height 
and  width  to  be  determined  independently. 


The  total  DX  center  density  in  the  barriers 
=  3.8X  10‘^  cm“'^)  is  given  by  A  F“*,  using  Eq.  (1)  and 
the  computed  /-  V curves.  The  well  width  (38  ±  3  A),  barrier 
height  (160  ±  20  meV),  and  barrier  width  (40  ±  5  A)  are  the 
parameters  that  give  best  agreement  between  the  observed 
and  calculated  curves  of  Fig.  2.  These  structural  parameters 
agree  with  estimates  from  growth  conditions,  whereas  A^o 
compares  well  with  our  estimate  of  total  Si  concentration  in 
the  barriers.  This  Nq  measurement  indicates  that  a  large 
fraction  of  the  Si  in  AlGaAs  occurs  as  DX  centers,  in  agree¬ 
ment  with  most  other  studies.  The  value  of  barrier  height 
gives  an  aluminum  fraction  of  0.18  ±  0.02,  in  good  agree¬ 
ment  with  the  target  value  of  0.2  ±  0.05,  assuming  the  pres¬ 
ently  accepted  value  of  60%  of  the  band-gap  difference  in  the 
conduction  band. 

Repeating  the  experiment  at  different  wavelengths 
yields  the  frequency  dependence  of  the  ionization  cross  sec¬ 
tion.  Figure  4  shows  the  relative  ionization  cross  section  in 
the  measured  range.  The  solid  line  is  that  found  for  Si  DX 
centers  in  AlGaAs  by  Lang  and  Logan.’  There  seems  to  be 
little  doubt  that  the  persistent  photoconductive  effect  we  ob¬ 
serve  arises  from  Si  DX  centers  in  the  AlGaAs  barriers. 

We  have  shown  that  a  simple  model  can  explain  persis¬ 
tent  photoconductivity  in  resonant  tunneling  structures. 
This  gives  us  some  faith  that  the  characteristics  crucial  to  the 
design  of  optimum  device  structures  can  now  be  obtained 
from  straightforward  measurements,  aiding  in  the  design 
and  fabrication  of  useful  devices.  The  persistent  photocon¬ 
ductivity  demonstrated  here  is  just  one  of  several  possibly 
useful  phenomena  (e.g.,  for  image  storage  or  memory  ele¬ 
ments).  We  have  also  observed  fast  photoconductivity  in 
quantum  well  resonators,  which  will  be  reported  in  the  fu¬ 
ture. 

Help  in  fabrication  of  these  devices  was  graciously  pro¬ 
vided  by  J.  J.  Lambert,  I.  H.  Mroczkowski,  and  N.  Usiak. 
We  are  also  indebted  to  S.  J.  Allen  of  Bell  Communications 
Research  and  to  A.  R.  Calawa,  B.  A.  Vojak,  and  D.  D.  Peck 
for  advice  and  assistance.  We  would  particularly  like  to 
thank  P.  E.  Tannenwald  for  useful  discussions.  This  work 
was  supported  by  the  U.S.  Army  Research  Office. 
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Resonant  tunneling  structures  have  been  shown  to  be  capable  of  high 
speed  charge  transport  as  well  as  microwave  oscillations. ^>2  More  recently^ 
persistent  photoconductivity  has  been  observed  and  has  been  used  to  improve 
model  calculations  of  the  electrostatic  potential  and  hence  of  the 
current-voltage  (I-V)  curves.  We  will  show  here  that  the  model  can  be  used 
with  the  persistent  photoconductivity  measurements  to  derive  the 
transconductance  of  a  transistor  based  on  the  double-barrier  resonant 
tunneling  structure. 

Figure  1  shows  the  conduction  band  of  a  typical  double-barrier  resonant 
tunneling  structure.  AlGaAs  is  used  to  provide  thin  (-50  A)  tunneling 
barriers  in  GaAs.  Between  these  barriers  is  a  thin  (-50  A)  layer  of  GaAs. 
Electrons  in  this  layer  are  confined  by  the  barriers  and  hence  are  resonant 
levels.  When  a  voltage  is  applied,  a  peak  in  the  transmitted  current  occurs 
when  the  incident  energy  of  an  electron  in  the  outer  electrode  region 
coincides  with  one  of  these  resonant  levels.  This  process  is  very  similar 
to  the  action  of  a  Fabry-Perot  etalon,  in  which  the  AlGaAs  tunnel  barriers 
serve  as  partially  transmitting  mirrors,  and  the  energy  of  the  electrons 
plays  the  role  of  the  photon  energy. 

It  is  intriguing  to  speculate  on  the  possibility  of  attaching  an 
electrical  connection  to  the  central  GaAs  well,  so  that  a  three-terminal 
device,  a  resonant  tunneling  transistor,  is  formed.  Figure  2  shows  one 
conceivable  geometry  for  such  a  device.  We  will  assume  that  the  emitter  is 
on  the  top  surface,  that  the  well  acts  as  the  base  by  controlling  current 
flow  between  the  emitter  and  collector,  and  that  the  collector  is  the  bottom 
electrode.  Tb  reduce  base  resistance,  the  spacing  between  the  base  pads 
should  be  minimized.  The  n"  regions  are  essentially  intrinsic,  thus 
reducing  base-to-emitter  and  base-to-col lector  capacitances.  If  these 
regions  are  kept  thin,  the  current  density  will  remain  unchanged  and  the 
carrier  transit  time  can  be  kept  short  compared  to  other  time  constants  of 
the  device.  The  conduction  band  at  zero  bias  is  shown  on  the  left  of  the 
figure.  We  note  that  contact  to  the  very  thin  base  region  will  probably  be 
a  major  technological  challenge,  but  may  be  worth  the  effort  in  view  of 
results  presented  later  in  this  paper.  We  will  estimate  the  operating 
characteristics  of  this  type  of  transistor  from  the  two-terminal  results 
using  the  effects  of  persistent  photoconductivity  to  simulate  the  third 
terminal . 


*This  work  was  supported  by  the  U.S.  Army  Research  Office. 


The  persistent  effect  of  light  on  the  quantum  well  I-V  curves  is  shown 
in  Figure  3.  At  temperatures  below  20  K  many  hours  are  required  for  the 
current  of  an  illuminated  device  to  decay  to  that  of  an  unexposed  device. 

The  solid  I-V  curve  is  unexposed,  the  long-dashed  curve  has  had  an  arbitrary 
exposure,  Q,  the  short-dashed  curve  has  exposure  10  Q,  and  the  broken  top 
curve  has  exposure  100  Q  which  nearly  saturates  the  effect.  After  exposure 
it  is  necessary  to  warm  the  sample  to  nearly  room  temperature  and  then 
recool  to  restore  the  unexposed  state.  The  model  which  best  explains  these 
results  is  sketched  in  Fig.  4.  It  is  well  known  that  complexes  called  DX 
centers  exist  in  n-type  AlGaAs  which  exhibit  a  persistent  ionization  upon 
illumination.  These  ionized  DX  centers  produce  positive  space  charge  in  the 
barriers  and  the  electrons  which  have  been  freed  in  the  process  find 
themselves  in  the  electrode  regions  because  their  energy  is  lower  there  than 
either  in  the  barriers  or  well.  The  dipoles  thus  formed  have  the  effect  of 
lowering  the  barriers  and  potential  of  the  resonant  energy  level  in  the  well 
with  respect  to  the  electrons  in  the  electrodes,  which  is  similar  to  the 
result  of  applying  a  potential  directly  to  the  well  (base).  This  shifts  the 
peak  current  to  lower  applied  voltages  and  increases  the  transparency  of  the 
barriers. 

Figure  5  shows  a  comparison  between  the  observed  and  calculated  I-V 
curves.  The  shift  of  the  peaks  to  lower  voltages  is  evident,  and  the 
increase  in  peak  current  and  zero-voltage  conductance  are  a  result  of  the 
lower  effective  barriers.  Notice  that  the  observed  current  has  a  component 
which  continues  to  increase  at  higher  voltages.  A  further  complication  is 
that  this  component  increases  with  increasing  concentration  of  ionized 
centers.  This  leakage  current  probably  arises  from  impurity  assisted 
tunneling  through  the  barriers.  It  will  be  a  source  of  uncertainty  in  the 
calculation  of  transconductance  which  follows.  It  should  also,  be  noted  that 
the  calculated  current  density  is  much  higher  than  that  observed.  Higher 
current  densities  can  be  expected  if  the  actual  fabrication  more  closely 
approaches  the  ideal  structure  on  which  our  model  is  based. 

The  transconductance  of  a  three-terminal  device  can  be  calculated  from 
our  measured*  I-V  curves.  The  central  idea  involves  the  fact  that  the 
potential  of  the  well  (base)  is  lowered  relative  to  the  electrodes  (emitter 
and  collector)  upon  exposure  to  light  as  shown  in  Fig.  5.  The  shift  in  this 
potential,  A<|>,  is  easily  calculated  from  the  voltage  shift,  AVp,  in  the 
current  peak  and  is  given  by 

AVp  =  2A<fr  (1) 

The  factor  of  2  arises  from  the  fact  that  the  base-to-emitter  applied 
voltage  is  half  of  that  applied  between  the  base  and  collector.  The  change 
in  collector  current  density,  AJc»  due  to  this  change  in  base-to-emitter 
voltage,  A<(),  is  given  by  the  difference  in  current  between  the  two  curves. 
The  transconductance,  G^n,  is  therefore  given  by 

Gjn  -  A  AJ^/A^  =  2A  Aj^/AVp 


(2) 


where  A  is  the  base  area.  There  is  one  complication,  however.  Since  the 
leakage  current  discussed  above  also  changes  with  a<(>,  it  is  difficult  to 
determine  and  account  for  this  contribution.  Nevertheless,  it  seems 
reasonable  to  assume  that  at  voltages  below  the  resonance  the  output 
impedance  is  large  and  positive.  We  have  therefore  extrapolated  the  excess 
current  at  high  voltage  smoothly  back  toward  zero  keeping  the  output 
impedance  positive.  These  extrapolations  are  shown  in  Fig.  6.  Using  the 
I-V  curves  of  Fig.  6  we  obtain  an  intrinsic  transconductance  of  2000  mS/mm 
for  an  emitter  width  of  1  un.  Specific  contact  resistances  on  the  order  of 
10"7  ohm  cm2  gpg  necessary  to  realize  an  extrinsic  transconductance  of  this 
order. 


It  is  interesting  to  calculate  the  transit  time  through  the  active 
region.  Suppose  we  use  intrinsic  regions  near  the  emitter  and  collector  of 
400  A  each  to  minimize  capacitance.  Then  at  a  saturation  velocity  of 
10^  cm/s  an  electron  will  cross  the  device  in  about  1  ps.  In  fact,  the 
velocity  would  probably  be  much  higher  since  transport  across  such  short 
distances  will  be  nearly  ballistic. 

Another  parameter  of  interest  is  the  small  signal  delay  time 

t  =  ^/C 

where  C  is  the  total  input  capacitance.  Using  the  geometry  and  parameters 
described  above,  we  find  a  gate  delay  time  of  about  3  ps,  or  a  cutoff 
frequency  of  about  55  GHz. 

Finally,  we  will  examine  the  base  resistance.  Since  the  base  is  very 
thin  and  lightly  doped,  it  might  be  expected  that  the  base  resistance  would 
be  quite  high.  However,  there  are  some  fortunate  aspects  to  the  quantum 
well  geometry  which  reduce  its  resistance  considerably.  It  is  well  known 
that  electron  mobilities  in  modulation-doped  quantum  wells  are  very  high  at 
temperatures  of  100  K  and  below.  Although  this  structure  is  not  a 
modulation-doped  well  in  the  usual  sense,  the  base  layer  acquires  charge 
from  the  tra'nsiting  carriers  without  any  direct  doping.  Just  as  in  a 
Fabry-Perot  resonator,  near  the  resonant  transmission  peak  carriers  build  up 
between  the  barriers  to  a  level  much  higher  than  the  average  transiting 
carrier  density.  For  present  samples  we  estimate  an  average  well  charge 
density  of  approximately  lO^^  cm"3.  Using  the  best  reported^  mobility  at 
50  K  of  5  X  10^  cm2/V-s,  a  base  carrier  density  of  2  x  10l7  cm"^^  an  emitter 
width  of  1  »JT>  and  length  of  20  im,  and  a  40-A  well,  we  find  a  base 
resistance  of  less  than  10  ohms.  This  gives  an  RC  product  for  base 
resistance  and  collector  capacitance  of  about  0.5  ps,  so  base  resistance 
does  not  seem  to  be  a  serious  limitation  on  performance. 

The  charge  in  the  base  layer  will  build  up  or  decay  by  electrons 
tunneling  in  or  out  of  the  well.  The  characteristic  time  for  this  process 
is  not  easily  calculated  because  the  well  size  is  comparable  to  the  electron 


de  Broglie  wavelength.  However,  an  exact  calculation  by  Dargys  and 
Kundrotas®  for  an  electron  in  a  delta-function  potential  with  an  applied 
electric  field  gives  a  value  of  this  characteristic  time  of  order  s 

for  a  GaAs  well  in  AlGaAs  and  electric  fields  comparable  to  those  in  our 
structures  near  the  resonance  peak.  The  actual  time  required  for  our 
structure  may  even  be  less  since  our  barriers  have  finite  thickness. 

In  view  of  the  promising  estimations  which  have  been  made  possible  by 
the  observation  of  persistent  photoconductive  effects,  it  seems  reasonable 
to  pursue  the  idea  of  a  resonant  tunneling  transistor  further.  Contact  to 
the  thin  base  region  appears  to  be  the  most  difficult  task  at  the  moment, 
but  if  this  could  be  overcome,  the  structure  suggested  here,  as  well  as 
other  similar  designs,  could  be  fruitfully  investigated  experimentally. 
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Figure  1.  Electron  energy  as  a  function  of  position  in  the  quantum 
well  resonator.  The  energy  level  Ej  occurs  above  the 
bottom  of  the  conduction  band  because  of  confinement  in 
the  X  direction. 


Figure  2.  Proposed  geometry  for  a  resonant  tunneling  transistor. 
The  spatial  dependence  of  the  conduction  band  energy 
is  also  shown. 


Figure  3.  Current-voltage  curves  of  a  quantum  well  resonator  after 
exposure  to  light,  but  after  the  light  has  been  turned 
off.  Four  different  exposures  are  shown. 
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Figure  4.  Electron  energy  as  a  function  of  position  in  a  quantum 
well  resonator  showing  the  effect  of  positive  charge 
in  the  barriers.  The  dashed  lines  show  the  conduction 
band  in  the  absence  of  ionized  centers. 


Figure  5.  Observed  and  calculated  I-V  curves  for  the  structure 
shown  in  Figure  4.  The  exposures  are  the  same  as 
in  Figure  4. 
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Figure  6 


Method  of  calculating  the  transconductance  of  the 
proposed  transistor  from  two  terminal  I-V  measurements 
(see  text). 
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Picc^ond  bistable  operation  has  been  experimentally  observed  for  the  first  time  in  a  double¬ 
barrier  resonant  tunneling  diode.  A  rise  time  of  2  ps  was  measured  using  the  electro-optic 
sampling  technique;  this  is  the  fastest  switching  event  yet  observed  for  an  electronic  device. 

This  time  domain  measurement  adds  necessary  information  to  the  understanding  of  the 
transport  mechanisms  in  the  resonant  tunneling  diode  and  is  consistent  with  switching  time 
limiUtions  computed  for  the  device.  It  also  demonstrates  that  appropriately  designed  double- 
barrier  quantum  well  diodes  have  a  response  time  comparable  to  that  of  the  fastest  all-optical 
logic  elements,  and  that  they  may  be  very  useful  in  high-speed  logic  applications. 


Since  its  inception  by  Tsu  and  Esaki,*  the  notion  of 
quantum  mechanical  tunneling  through  double-barrier  he¬ 
terostructures  has  raised  considerable  interest.  The  negative 
differential  resistance  (NDR)  exhibited  by  this  device 
makes  it  suitable  for  a  wide  variety  of  applications;  among 
them  are  oscillators,^  mixers,*^  and  transistors/-^  Especially 
noteworthy  is  the  experimental*  and  theoreticaP  evidence 
suggesting  that  these  devices  could  operate  in  the  terahenz 
frequency  regime.  We  report  here,  for  the  first  time,  the  mea¬ 
surement  of  a  double-barrier  resonant  tunneling  diode 
( RTD )  in  its  bistable  or  switching  mode  with  sufficient  tem¬ 
poral  resolution  to  resolve  the  switching  time.  We  have  mea¬ 
sured  this  switching  time  to  be  less  than  2  ps.  This  time  com¬ 
pares  favorably  with  the  theoretical  considerations  of  Liu 
and  Coon*  for  the  switching  time  of  a  heterojunction  double¬ 
barrier  diode  circuit.  Practical  electronic  devices  that  re¬ 
spond  on  such  a  short  time  scale  would  rival  the  speed  of  the 
fastest  optical  switching  elements,  leading  to  their  utility  in 
electronic  logic  circuits  of  the  future.  These  time  domain 
measurements  may  also  provide  some  insight  concerning  the 
mechanisms  of  electron  transport  through  double-barrier 
heterostructures,  processes  that  arc  far  from  being  complete¬ 
ly  understood. 

The  RTD  used  in  these  measurements  was  grown  by 
molecular  beam  epitaxy  to  have  two  barriers  of  1.5-nm-thick 
AlAs  separated  by  a  GaAs  well  4.5  nm  thick.  The  outer 
regions  of  GaAs  contained  Si  doping  to  an  electron  density 
2x  10*^  cm”^  This  resulted  in  a  peak  current  density  of 
about  4x10^  A/cm^  For  the  4-/im-diam  RTD  tested  here, 
the  capacitance  when  biased  near  the  current  peak  was  cal¬ 
culated  to  be  20  fF,  and  the  scries  resistance  was  15  n.  More 
details  of  the  material  growth  and  fabrication  can  be  found  in 
Ref.  9. 

The  device  was  tested  using  an  instrument  with  the  sub¬ 
picosecond  resolution  necessary  to  measure  such  brief 
events,  the  electro-optic  sampling  system.***  A  schematic 
diagram  of  the  circuit  used  to  test  the  time  response  of  the 
RTD  is  shown  in  Fig.  1.  Coplanar  electrodes  of  gold  acting 
as  a  transmission  line  were  deposited  on  both  a  GaAs  switch 
and  a  lithium  tantalate  (LiTaOj)  sampling  crystal,  and  the 


chip  containing  the  array  of  resonant  tunneling  diodes  was 
mounted  so  that  the  device  faced  the  switching  clement  as 
shown.  An  ohmic  contact  on  one  of  the  4-/im-diam  diodes 
was  connected  to  an  electrode  using  a  gold-coated  phosphor- 
bronze  wire  that  had  been  etched  to  a  1  /im  diameter  at  its 
tip.  Another  ohmic  contact  on  the  back  of  the  RTD  was 
fixed  to  one  coplanar  electrode  on  the  LiTaO,  crystal  with 
silver  epoxy,  and  the  contact  was  completed  with  the  pointed 
wire  (whisker)  such  that  the  device  was  in  series  with  the 
transmission  line.  Synchronized  optical  pulses  from  a  dye 
laser  of  about  80  fs  duration  were  used  to  activate  a  photo- 
conductive  switch  at  the  input  to  the  device  and  to  measure 
the  electric  fields  across  the  transmission  line  by  probing  the 
change  in  index  of  refraction  in  the  LiTaO,  induced  by  the 
output  of  the  device.  All  the  measurements  reported  were 
made  at  room  temperature. 

The  equivalent  circuit  of  the  RTD  can  be  represented  as 
a  capacitor  in  parallel  with  a  variable  resistor  (which  can  be 
negative ) ,  and  these  in  series  with  a  small  positive  resistor,  as 
shown  in  Fig.  2.  The  RTD  was  connected  in  scries  with  a 
load  resistance,  a  whisker  inductance,  and  a  time-varying 
photoconductivc  switch  resistance  R{t),  as  shown.  The  dc 


eicitation 


FIG.  I.  Experimental  configuration  for  electro-optic  sampling  of  fesonani* 
tunneling  diode  A  12-/im-diam  wire  with  a  I /im  tip  contacts  a  mesa  on  the 
resonant-tunneling  duxJe  chip. 
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no.  2.  Equivalent  circuit  for  resonant*tunneiing  diode  and  test  circuit. 


bias  to  the  RTD  was  applied  through  a  large  inductor,  and 
the  bias  from  the  photoconductive  switch  was  considered  as 
the  input  to  the  network.  Any  other  parasitic  circuit  ele¬ 
ments  have  been  ignored.  When  an  additional  current  A/ 
was  presented  to  the  RTD  as  the  switch  was  illuminated  and 
closed,  a  potential  difference  developed  across  the  diode  and 
its  operating  point  was  changed.  As  shown  in  the  second  and 
third  frames  of  Fig.  3(a),  the  load  line  was  thus  shifted  up¬ 
wards  on  the  current-voltage  (I-V)  characteristic  from 
point  A.  When  the  load  line  exceeded  the  current  peak  of  the 
I-V  curve,  the  stable  operating  point  became  px)inl  C  (see 
Fig.  3(a)  ],  and  switching  was  initiated.  The  switching  time 
depended  on  how  rapidly  the  device  capacitance  could  be 
charged,  resulting  in  a  change  of  voltage  from  that  of  point  A 
to  point  C. 

The  method  used  to  determine  the  switching  time  of  the 


device  will  now  be  presented.  The  region  between  the  two 
load  lines  labeled  •‘1”  in  each  frame  of  Fig.  3(a)  represents 
the  shift  in  load  line  resulting  from  closing  the  switch  with 
the  optical  pulse,  when  starting  from  different  dc  bias  condi¬ 
tions.  The  region  labeled  “2,”  shown  in  the  top  frame  only, 
indicates  the  shift  in  load  line  along  a  relatively  linear  por¬ 
tion  of  the  curve.  The  resulting  waveform  from  region  2 
at  the  output  of  the  diode  was  used  as  a  reference  so  that  a 
comparison  could  be  made  to  each  of  the  other  waveforms 
from  region  1  and  the  switching  response  more  easily  re¬ 
solved.  This  waveform  2  is  depicted  in  each  of  the  frames  of 
Fig.  3(b).  Displayed  as  a  dashed  line  is  the  waveform  from 
each  region  1  in  Fig.  3(a)  and,  below  the  waveforms,  the 
temporal  variation  of  the  difference  of  the  two  waveforms. 
The  waveforms  in  Fig.  3(b)  are  exaggerated  to  schematical¬ 
ly  show  the  effects  of  switching.  The  experimental  result  of 
subtracting  the  two  waveforms  from  regions  1  and  2  is  pre¬ 
sented  in  Fig.  3(c).  In  the  first  frame,  the  cancellation  is 
nearly  complete.  The  discrepancies  observed  are  due  to  noise 
and  to  small  deviations  in  the  regions  of  the  /-K  curve  tra¬ 
versed.  The  switching  process  can  be  observed  in  the  second 
and  third  frames  of  Figs.  3(b)  and3(c).  In  the  second  frame, 
curve  1,  the  dashed  line,  follows  curve  2  only  from  points  A 
to  B,  where  it  coincides  with  the  rising  portion  of  the  I-V 
curve.  When  the  load  line  exceeds  the  current  peak  of  the  /-  V 
curve,  the  diode  switches  to  point  C  and  a  rapid  drop  in 
current  is  observed.  In  the  second  frame,  as  the  input  de¬ 
creases,  the  optical  switch  current  also  decreases  so  that  the 
load  line  drops  below  point  E,  and  the  operating  point  re¬ 
turns  to  F  and  relaxes  to  A,  where  waveforms  1  and  2  are 
again  nearly  identical.  This  is  indicated  in  the  difference 
waveform  (1—2)  of  Fig.  3(b)  and  can  also  be  seen  in  the 
experimental  difference  waveform  in  Fig.  3(c). 


FIG.  3.  Rcsonant'lunneling  diode 
switching  as  a  function  of  bias:  (a)  /- 
Feurve  with  load  lines,  (b)  analytic 
representation  of  waveforms  result¬ 
ing  from  movement  of  load  lines  in 
(a)  and  the  difference  of  these  wave¬ 
forms.  and  (c)  experimental  signal 
difference. 
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The  third  frame  of  Fig.  3  displays  the  same  switching 
action,  but  with  the  device  attaining  a  latched  state.  The 
operation  from  points  C  ioD  to  £  is  along  a  path  of  the  /-  V 
curve  that  would  approximately  cancel  with  that  from  re¬ 
gion  2,  leading  to  a  long  plateau  on  the  difference  waveform 
and  a  device  that  has  been  switched  into  the  lower  current 
state.  This  was  distinctly  observed  in  frame  3  of  Fig.  3(c). 
Since  the  peak  of  the  /-  V  curve  was  reached  by  the  load  line 
earlier  than  in  the  frame  above,  the  switching  took  place  at 
an  earlier  time,  as  evidenced  by  the  position  of  the  wave¬ 
forms  relative  to  the  thin  vertical  line  of  Fig.  3(c)  that  repre¬ 
sents  a  constant  time  reference.  In  frame  2,  the  onset  of  the 
switching  began  at  the  marker,  while  in  frame  3  the  switch¬ 
ing  event  was  already  about  halfway  completed.  The  rise 
time  of  the  switching  event,  as  measured  from  the  third 
frame  of  Fig.  3(c),  was  discovered  to  be  1.9  ps  in  duration 
between  the  10%  and  90%  points.  During  the  relatively  long 
time  between  optical  pulses,  the  bias  point  returns  to  point  A. 
This  is  possible  because  the  impedance  determining  the  low- 
frequency  load  line  is  much  smaller  than  the  high-frequency 
impedance.  In  frame  4  of  this  figure,  the  dc  bias  point  is 
extended  above  the  area  of  negative  differential  resistance, 
where  region  1  traverses  virtually  the  same  slope  as  region  2, 
and  cancellation  is  again  expected  and  found. 

As  mentioned  earlier,  the  switching  time  should  depend 
on  the  time  to  charge  the  device  capacitance  between  the 
voltages  corresponding  to  points  A  and  C  in  Fig.  3(a)^’*: 

ts^C(Vc-VB)/(Ip-U.  (I) 

where  and  are  the  peak  and  valley  currents  in  the  nega¬ 
tive  differential  resistance  region.  As  measured  experimen- 
^  Charge  stor¬ 

age  occurs  predominantly  outside  the  double-barrier  region, 
and  so  to  compute  the  capacitance,  we  must  know  the  dis¬ 
tance  between  the  conductive  regions  on  each  side  of  the 
quantum  well.  The  GaAs  well  width  was  4.5  nm,  while  the 
AlAs  barriers  were  1.5  nm  thick,  and  the  depletion  width 
into  the  doped  collector  region  was  approximately  70  nm. 
Using  a  relative  permittivity  for  GaAs  of  13.1,  the  capaci¬ 
tance  was  found  to  be  about  20  fF,  and  the  rise  time  for 
switching  was  calculated  to  be  2.9  ps.  Bearing  in  mind  that 
Eq.(l)  is  only  approximately  valid  and  that  the  instrumen¬ 
tal  contribution  was  not  measured,  this  is  in  reasonable 
agreement  with  the  measured  value  of  1.9  ps. 

The  tunneling  current  in  a  RTD  can  be  a  result  of  a 
resonant,  coherent  tunneling  current,  as  in  the  optical  ana¬ 
log  of  a  Fabry-Perot  resonator,  and/or  a  sequential,  inco¬ 
herent  tunneling  process  stemming  from  elastic  and  inelastic 
collisions.  Originally,  the  idea  of  electrons  scattering  and 
sequentially  tunneling  was  used  by  Luryi”  to  explain  the 
discrepancy  between  the  measured  and  expected  current  of 
the  device.  It  was  pointed  out  by  Weil  and  Vinter*^  and  Jon- 
son  and  Grinewajg*^  that,  in  most  cases,  the  tunneling  cur¬ 
rent  is  proportional  only  to  the  resonant  level  linewidth 
and  insensitive  to  the  scattering  energy  width  .  The  device 
response  time,  however,  is  inversely  proportional  to  the  total 
energy  width  of  the  system,  F,  =  F,.  -f  F^.  The  measured 
switching  time  can  be  looked  at  as  an  upper  limit  of  the 
response  time  of  the  structure,  since  the  measurement  circuit 


could  only  increase  the  time  of  response.  In  our  experiments, 
the  device  tested  had  a  geometry  such  that  F^  was  about  2 
meV.  The  resulting  device  response  from  resonant  tunneling 
alone  (r:::^  2^/F, )  could  therefore  be  of  the  order  of  a  single 
picosecond.  Thus  we  can  make  no  statement  about  the  im¬ 
portance  of  sequential  tunneling  from  this  measurement. 
However,  if  a  RTD  with  a  much  smaller  F,  were  to  show  a 
switching  time  of  a  picosecond  or  so,  then  we  could  conclude 
that  sequential  tunneling  was  the  dominant  transport  mech¬ 
anism,  since  this  process  would  give  us  the  very  fast  response 
commensurate  with  the  collisionally  broadened  F^ . 

In  conclusion,  we  have  used  the  electro-optic  sampling 
technique  to  measure  the  fastest  switching  event  yet  for  an 
electronic  device.  The  switching  speed  of  the  RTD  is  also 
comparable  to  the  one  observed  in  the  fastest  optical  bistable 
devices.  The  similarity  in  the  observed  performance  of  the 
electronic  and  optical  bistable  devices  is  not  surprising  con¬ 
sidering  the  fact  that  the  switching  time  is  limited  by  the 
round-trip  time  in  the  Fabry-Perot  resonant  cavity.  That  is, 
the  nanometer  dimensions  and  slow  electron  velocity  of  a 
RTD  provide  comparable  round-trip  times  to  those  obtaina¬ 
ble  with  the  submicron  dimensions  but  higher  photon  veloc¬ 
ities  of  multiple  quantum  well  optical  devices.  This  illus¬ 
trates  just  how  similar  electrical  and  optical  resonators  are, 
and  how  the  response  of  electronic  devices  is  catching  up  to 
the  supposedly  faster  optical  devices. 
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Abstract 

The  excellent  sensitivity  of  single  planar  mixers  at  140  GHz  fTa  •  4900K  DSB)  has  led  us  to  investi¬ 
gate  arrays  of  mixers »  each  fed  by  planar  dipole  antennas.  Model  measurements  of  impedances  and 
mutual  coupling,  as  well  as  fabrication  and  tests  at  millimeter  and  submillimeter  wavelengths,  are 
described. 


Introduction 

We  have  been  investigating  planar  mixers  for  focal  plane  imaging 
at  millimeter  and  submillimeter  wavelengths  We  report  here 
the  progress  in  that  effort,  with  new  results  at  140  GHz  and  at 
1760  GHz,  as  well  as  new  scale  model  results  at  3  GHz. 

Our  design  approach,  described  in  References  1  and  2,  has  been 
the  use  of  two  microstrip  dipole  antennas  on  top  of  a  grounded 
crystal  quartz  substrate,  whi(^  feed  a  nonlinear  mixing  elemenL 
Crystal  quartz  was  selected  because  of  its  availability  in  thin 
wafers  (^150  pm)  and  its  relatively  low  dielectric  constaru.  Re¬ 
cently  we  have  been  using  high  performance  GaAs  beam-lead 
Schotiky  diodes  from  AIL 

The  basic  circuit  is  shown  in  Figure  1.  The  antennas  are  two 
full-wavelength  dipoles  spaced  by  approximately  X/2  with  the 
diode  bonded  midway  between  the  antennas.  The  antenna  pair 
acts  as  a  2  element  array  which  radiates  primarily  in  the  direction 
perpendicular  to  the  substrate.  A  low-pass  filter,  designed  to 
reduce  fundamental  and  second  harmonic  leakage,  is  partiaUy 
shown  at  the  right  edge  of  Figure  1.  Not  shown  is  the  50  Q 
coupled-microstrip  transmission  line  from  the  filter  to  the  edge  of 
the  substrate  which  carries  the  DC  bias  and  IF  signal.  Both  the 
filter  and  transmission  line  have  been  redesigned  based  on  an  im¬ 
proved  theoretical  model. 


Scale  Model  Measurements 

In  order  to  continue  to  improve  the  performance  of  this  circuit, 
scale  models  were  constructed  and  tested  at  S-band  (2-4  GHz). 
Measurements  were  first  made  of  the  impedance  seen  at  the  diode 
position  of  the  double-dipole  circuit  This  measurement  was  made 
through  a  V4  balun  with  no  diode  in  the  drcuit  A  value  of  38 
j23  Q  was  obtained  at  3.2  GHz,  the  scaled  center  frequency  of  the 
140  GHz  circuit  The  measured  impedance  can  be  underwood  as 
the  high,  real  impedarKe  of  a  pnllel  pair  of  full  wavelength  di¬ 
poles  transformed  by  a  length  of  100  Q  coupled  microstrip 
nansmission  line.  According  to  theory,  at  3.2  GHz  the  phase  shift 
benveen  the  diode  port  and  each  dipole  is  1 12*.  The  extra  phase 
shift  beyoixl  90*  is  responsible  for  40%  of  the  iixluctive  reactance 
observed.  The  additional  inductance  may  be  due  to  the  balun  or 
to  operation  at  a  frequency  slightly  off  resonance. 


A  measurement  of  the  impedance  seen  at  the  diode  pon  is  impor¬ 
tant  for  improving  the  march  between  the  feed  and  the  diode.  For 
optimum  coupling  id  this  diode  at  millimeter  wavelengths,  the 
resistance  presented  by  the  feed  should  be  approximately  150  D. 
This  assumes  that  the  diode  resistance  will  be  100  to  200  D  when 
it  is  biased  for  minimum  conversion  loss  and  noise  The  op¬ 
timum  value  of  feed  inductance  depends  on  the  diode  capacitance 
and  the  frequency  of  operation.  Minor  changes  in  the  circuit  di¬ 
mensions  such  as  an  increase  in  the  dipole  feed  gap  or  a  decrease 
in  the  coupled  microstrip  line  length  could  provide  a  significantly 
better  impedance  match. 

Measurements  were  also  made  of  the  coupling  between  pairs  of 
closely  spaced  double-dipole  antennas  to  evaluate  the  suiubility  of 
the  structure  for  diffraction-limited  imaging.  Any  coupling 
between  pain  of  antennas  in  an  array  lessens  the  contrast  available 
when  imaging.  Two  configurations  were  tested  for  coupling 
between  pain:  the  H-ptane  orientation  (i.e.,  parallel  pain  of  di¬ 
poles)  and  the  E-plane  orientation  (x.e.,  coUinear  pain  of  dipoles). 
In  both  cases  the  center-to-ceroer  separation  was  Xo,  the  free-space 
wavelength.  We  measured  ISnl  -  -26  dB  and  -35  dB  in  Cie  E- 
and  H-planes,  respectively.  The  results  are  consistent  with  the 
value  predicted  by  Pozar  ^  for  pain  of  single  resonant  antemias  in 
the  E-ptane  configuration  (ISnl  *  -30  dB). 


Measurements  at  140  GHz 

We  have  fabricated  the  improved  circuit  on  150  pm  quanz  arvl 
tested  its  bandwidth,  radiation  pattern  and  heterodyne  sensitivity 
near  140  GHz.  Figure  2  shows  the  radiation  pattern  at  132  GHz 
plotted  on  a  linear  scale.  The  lack  of  nulls  or  significant  ripples  in 
the  main  beam  leads  us  to  believe  that  the  coupling  into  the  anten¬ 
nas  from  surface  waves  on  the  substrate  is  weak.  The  similar 
widths  of  the  E-  and  H-ptane  patterns  are  important  for  coiipling 
to  circular  lenses  and  reflecton. 

In  order  to  estimate  the  spectral  bandwidth  of  the  antenna,  video 
responsivity  measurements  were  made  at  several  different  frequen¬ 
cies.  A  p^  responsivity  of  170  V/W  was  measured  at  132  GHz 
and  at  least  85  V/W  was  measured  over  a  12%  bandwidth.  The 
theory  of  Pozar^  predicts  that  greater  bandwidth  could  be 
achieved  if  a  thicker  substrate  were  used;  however,  a  thicker  sub¬ 
strate  would  increase  the  surface  wave  coupling. 
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The  circuit  has  been  tested  as  a  planar  mixer  with  the  LO  and 
blackbody  RF  signal  combined  externally  in  a  quasi-opdcal  di- 
plexer.  Both  beams  were  focused  onto  the  substrate  by  an  offset 
ellipsoidal  mirror.  A  receiver  temperature  of  4900K  (DSB)  and  a 
total  conversion  loss  of  8.5  dB  were  measured  at  an  IF  frequency 
of  1.1  GHz.  The  mixer  required  4  mW  of  LO  power  from  the 
carcinotion  and  U  mA  bias  for  optimum  sensitivity.  This  LO 
power  level  suggests  that  high  power  sources  are  needed  as  local 
oscillators  if  large  arrays  of  this  type  are  to  be  driven  to  optimum 
sensitivity. 


Submillimeter  Measurements 

The  circuit  has  been  fabricated  on  70  and  33  pm  quaru  for  testing 
at  300  and  560  GHz,  respectively,  which  is  in  progress.  The 
diodes  used  typically  have  R,  -  3 J  Q  and  Cj(0)  -  9  fF,  yielding  a 
cutoff  frequency  of  5000  GHz.  In  order  to  examine  the  upper  fre¬ 
quency  limit  of  the  560  GHz  circuit,  we  have  tested  it  at  1760 
GHz  (171  pm)  using  a  methanol  laser  pumped  by  a  COi  laser. 
We  measur^  a  video  responsivity  of  5  mV/W.  At  171  pm  we 
believe  that  the  diode  is  acting  as  its  .own  antenna,  since  the  diode 
is  on  the  order  of  2X  in  length.  A  different  design  approach, 
such  as  monolithic  fabrication  on  a  thinner  substrate,  is  more  ap¬ 
propriate  at  this  wavelength  and  could  yield  improved  sensitivity. 


Conclusions 

We  have  developed  a  planar  millimeter  wave  mixer  fed  by  micros¬ 
trip  dipole  antennas  which  has  excellent  sensitivity  and  radiation 
patterns  at  140  GHz.  The  improvement  over  our  previous  results 
is  thought  to  be  due  to  better  low-pass  filter  and  transmission  line 
designs,  and  to  the  use  of  a  smaller,  lower  resistance  diode.  Scale 
modeling  at  3.2  GHz  has  been  used  to  estimate  the  circuit  im¬ 
pedance  at  the  diode  position  and  to  determine  the  coupling 
between  circuits  closely  spaced  in  an  array.  Measurements  at  171 
pm  suggest  that  a  different  consmicdon  approach  is  needed  at  this 
shon  wavelength. 
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Figure  1.  The  double-dipole  circuit  with  AIL  beam-lead  Schottky 
diode  bonded  midway  between  the  two  full-wavelength  dipoles. 
The  microstrip  circuit  shown  here  is  Au  on  150  pm  crystal  quartz 
for  operation  at  140  GHz. 


Figure  2.  The  radiation  pattern  at  132  GHz  of  the  antenna  pair 
shown  in  Figure  1. 
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Abstract 


A  compact,  solid  state  submill imeter-wave  heterodyne 
radiometer  has  been  developed  and  was  used  to  measure  spec¬ 
tral  characteristics  of  a  water  vapor  jet  in  a  space  simula¬ 
tion  chamber.  Features  of  the  557  GHz  water  vapor  line  pro¬ 
file  were  observed  in  significantly  greater  detail  than  in 
previous  experiments  through  an  increased  sensitivity  and 
improved  frequency  resolution  (600  kHz).  The  local  oscil¬ 
lator  of  the  radiometer  consisted  of  a  frequency  multiplica¬ 
tion  chain  which  was  driven  by  an  InP  Gunn  oscillator  at 
92.6  GHz,  and  which  contained  a  frequency  tripler  and  har¬ 
monic  mixer  in  cascade.  The  front  end  of  the  receiver  had  a 
noise  temperature  of  4500  K  (DSB)  at  555  GHz,  consumed  3  W 
and  weighed  3  kg.  This  advance  in  technology  is  particu¬ 
larly  relevant  to  submillimeter-wave  radiotiietry  from  a 
space-based  platform. 

Key  words:  water  vapor  spectroscopy,  submillimeter-wave 
radiometry,  heterodyne  radiometer,  submillimeter-wave  oscil¬ 
lator,  frequency  multiplier. 


Introduction 


In  the  last  few  years  at  Lincoln  Laboratory  we  have 
conducted  a  series  of  experiments  on  the  submillimeter  radi¬ 
ative  properties  of  a  jet  of  H2O  molecules  ejected  from  a 
sonic  nozzle  and  expanded  adiabatically  into  a  high-vacuum 
chamber  (l)-(3).  The  intensities  and  line  shapes  of  the  752 
and  557  GHz  rotational  transitions  of  the  cooled  H2O  mole¬ 
cules  were  measured  radiometrical ly  against  a  warm  back¬ 
ground  using  heterodyne  receivers  developed  here  and  at  the 
University  of  Massachusetts.  Although  the  radiometers  were 
operated  in  a  laboratory  environment,  the  experiments  were 
designed  to  simulate  heterodyne  spectroscopic  measurements 
that  could  be  carried  out  from  a  space-based  platform. 
Space  operation  would  allow  observations  in  areas  such  as 
astronomy  and  aeronomy  that  would  not  be  masked  by  atmo¬ 
spheric  water  vapor  absorption. 

The  752  GHz  experiment  used  a  far-IR  laser  pumped  by  a 
CO2  laser  as  the  local  oscillator  (LO)  (1),  while  the 
557  GHz  measurements  were  made  with  a  harmonic  mixer  and  a 
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282  GHz  carcinotron  as  the  LO  (3).  In  both  experiments  dif¬ 
ficulty  in  achieving  sufficient  LO  power  with  frequency 
stability  limited  the  quality  of  the  spectral  data.  More¬ 
over,  because  of  their  large  size,  weight  and  power  require¬ 
ments,  neither  LO  is  readily  adaptable  for  use  in  space. 

These  problems  have  prompted  us  to  develop  a  local 
oscillator  chain  for  the  submillimeter  band  which  uses  only 
solid  state  components.  With  this  new  oscillator  we  have 
demonstrated  the  operation  of  a  solid  state  submillimeter 
radiometer  near  557  GHz  with  best  DSB  noise  temperature  of 
4500  K,  front  end  (LO,  mixer  and  IF  amplifier)  weight  of 
3  kg  and  power  consumption  of  about  3  W.  The  radiometer  was 
used  to  repeat  the  water  jet  spectroscopy  experiment  in  the 
space  simulation  chamber,  yielding  results  substantially 
superior  to  the  two  previous  experiments. 


The  Radiometer 


A  block  diagram  of  the  present  radiometer  and  spectros¬ 
copy  experiment  is  shown  in  Figure  1.  The  local  oscillator 
of  the  radiometer  began  with  a  50  mW  InP  Gunn  oscillator  (4) 
at  92.58  GHz.  This  signal  was  fed  by  waveguide  to  a  fre¬ 
quency  tripler  containing  a  GaAs  Schottky  varactor 
diode  (5), (6).  The' varactor  diode  was  whisker-contacted  and 
mounted  in  a  crossed-guide  structure  with  each  waveguide 
terminated  by  a  movable  backshort.  The  tripler  produced 
approximately  3  mW  of  power  at  278  GHz,  which  was  fed  by 
waveguide  to  a  harmonic  mixer  (5).  A  diagram  of  the  har¬ 
monic  mixer,  which  contained  a  whisker-contacted  Schottky 
mixer  diode,  is  given  in  reference  3.  The  harmonic  mixer 
has  also  been  operated  as  a  frequency  doubler,  which  has 
permitted  the  chain  to  operate  as  a  sextupler.  When  driven 
by  an  IMPATT  oscillator,  this  chain  has  produced  50  wW  at 
556  GHz  (7).  The  two  multipliers  will  be  described  more 
fully  at  a  later  date  (8). 

Although  the  receiver  was  operated  at  557  GHz  in  the 
spectroscopy  experiment,  its  DSB  noise  temperature  improved 
from  7800  K  to  4500  K  when  the  signal  frequency  was  lowered 
to  555  GHz.  This  difference  was  due  to  the  fact  that  the 
mixer's  performance  was  dependent  on  the  LO  power  supplied 
to  it,  which  in  turn  varied  sharply  with  the  frequency- 
dependent  impedance  match  between  the  tripler  and  mixer. 
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In  the  spectroscopy  experiment  the  harmonic  mixer  com¬ 
bined  the  tripled  Gunn  signal  with  the  chopped  RF  signal 
containing  the  spectral  information  to  yield  an  intermediate 
frequency  (IF)  at  1.45  GHz.  This  IF  signal  was  amplified  by 
a  GaAs  MESFET  with  a  noise  temperature  of  50  K,  then  passed 
through  a  digitally-controlled  bandpass  filter  in  the  spec¬ 
trum  analyzer.  The  filter's  bandwidth  was  300  kHz  and  its 
center  frequency  was  stepped  through  the  spectral  line  in 
nominal  250  kHz  increments.  The  actual  step  size  included  a 
correction  for  the  frequency  drift  of  the  Gunn  oscillator 
with  respect  to  the  tenth  harmonic  of  an  X-band  cavity- 
stabilized  klystron  (stalo),  as  measured  by  the  frequency 
counter.  Although  this  tracking  filter  did  not  compensate 
for  rapid  changes  in  the  Gunn's  frequency  (jitter),  it  did 
correct  for  the  Gunn's  frequency  drift  due  to  thermal 
effects.  The  100  kHz  wide  jitter  in  the  Gunn  resulted  in  an 
overall  frequency  resolution  of  600  kHz.  This  could  be 
improved  significantly  by  further  frequency  stabilization. 
After  filtering,  the  signal  was  synchronously  detected, 
integrated  for  3  sec,  digitized  and  recorded  by  computer. 


Spectroscopy  Results 

As  in  the  previous  experiments  (l)-(3),  a  laboratory- 
scale  water  vapor  jet  from  a  1  cm  diameter  nozzle  was 
directed  along  the  axis  of  a  large  high- vacuum  chamber  with 
liquid  nitrogen  cryopanelled  walls.  Figure  2  shows  the 
highly  resolved  557  GHz  H2O  absorption  line  measured  against 
a  Hg  arc  continuum  at  a  point  50  cm  downstream  from  the  jet 
nozzle  for  a  nominal  flow  rate  of  0.5  g/sec.  To  illustrate 
the  substantial  improvement  in  sensitivity  and  resolution, 
the  line  profile  measured  previously  under  similar  condi¬ 
tions,  but  with  the  radiometer  using  a  carcinotron  LO  (3), 
is  sketched  for  comparison.  At  line  center,  the  signal  tem¬ 
perature  represents  the  adiabatically  cooled  rotational 
temperature  of  the  optically  thick  line.  The  minimum 
temperature  of  -75  K  agrees  with  the  value  reported 
earlier  (3^  The  linewidth  reflects  the  Doppler  broadening 
caused  by  the  range  of  flow  velocities  in  the  diverging 
molecular  beam,  and  its  nearly  flat  bottom  is  due  to  satu¬ 
rated  absorption  of  the  continuum  radiation  by  most  Doppler 
components.  The  steep  rise  and  fall  of  the  line  now  indi¬ 
cate  that  the  divergence  angle  of  the  vapor  jet  is  well 
defined. 
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Figure  2.  Optically-thick  H2O  absorption  line  profiles 
measured  at  50  cm  from  nozzle  against  Hg  arc  lamp  back¬ 
ground.  Comparison  is  made  with  data  from  earlier  experi¬ 
ment  using  a  carcinotron  LO. 


Figure  3  shows  the  unsaturated  H2O  absorption  line 
under  the  same  conditions  but  at  a  point  200  cm  from  the 
nozzle.  The  improved  resolution  has  contributed  at  least 
partly  to  the  observation  of  a  self-reversal  effect  more 
prominent  than  that  reported  originally  (2),  where  it  was 
attributed  to  Doppler-shifted  regions  of  varying  optical 
depths  and  temperatures  in  the  expanding  jet.  The  shape  and 
magnitude  of  the  self  reversal  were  again  sensitive  to  the 
particular  H2O  flow  rate  and  the  distance  downstream  from 
the  nozzle.  Without  additional  data,  it  is  not  possible  to 
separate  the  relative  contributions  of  jet  temperature  and 
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Figure  3.  Absorption  line  of  H2O  jet  illuminated  by  Hg  arc 
lamp  under  conditions  which  caused  a  self  reversal  around 
the  center  of  the  line.  Both  Figures  2  and  3  contain  cor¬ 
rections  for  the  constant  power  in  the  image  sideband  and 
for  2  dB  loss  in  the  TPX  chamber  window  (see,  e.g..  Ref.  3). 


opacity  to  the  self  reversal.  This  reflects  a  general  prob¬ 
lem  in  radiometry,  which  occurs  when  the  quantity  observed 
is  the  product  of  the  local  temperature  and  absorption, 
integrated  along  the  line  of  sight.  In  our  case  the 
detected  signal  is  the  sum  of  the  radiometric  temperature  of 
the  arc  lamp  attenuated  by  the  cold  water  vapor  and  the 
emission  from  the  vapor.  Additional  measurements,  other 
than  radiometric,  must  be  made  to  more  fully  understand  the 
dynamics  of  a  vapor  jet  expanding  adiabatically  into  a 
vacuum.  However,  the  present  radiometer  is  well  suited  for 
obtaining  the  necessary  spectroscopic  information. 
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Conclusion 


The  development  of  a  solid  state  submillimeter-wave 
oscillator  has  expanded  the  range  of  possible  experiments  in 
submillimeter  radiometry  and  spectroscopy.  Earlier  we  des¬ 
cribed  the  oscillator's  operation  as  a  "stand-alone"  source 
which  could  be  useful  for  direct-detection  spectroscopy  (7), 
and  here  we  have  discussed  its  application  as  the  LO  of  a 
compact  heterodyne  radiometer  for  high-resolution  spectros¬ 
copy.  The  sensitivity  and  stability  of  this  radiometer  are 
close  to  state-of-the-art  (9),  and  the  combination  of  these 
factors  with  the  receiver's  compactness  and  efficiency  make 
this  a  promising  system  for  use  in  space.  Space  operation 
may  require  minor  modifications  in  the  receiver  design, 
although  both  multipliers  used  in  this  experiment  have 
proved  to  be  rugged  and  reliable,  having  survived  three 
years  of  laboratory  tests  without  a  failure. 
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